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Preface 


This manual provides users of the VMS operating system with detailed usage 
and reference information on mathematics routines supplied in the MTH$ 
facility of the Run-Time Library. 

Run-Time Library routines can only be used in programs written in languages 
that produce native code for the VAX hardware. At present, these languages 
include VAX MACRO and the following compiled high-level languages: 


VAX Ada 
VAX BASIC 
VAX BLISS-32 
VAX C 
VAX COBOL 
VAX COBOL-74 
VAX CORAL 
VAX DIBOL 
VAX FORTRAN 
VAX Pascal 
VAX PL/I 
VAX RPG 
VAX SCAN 


Interpreted languages which can also access Run-Time Library routines 
include VAX DSM and DATATRIEVE. 


Intended Audience 

This manual is intended for system and application programmers who want 
to call Run-Time Library routines. 


Document Structure 

This manual is organized into two parts as follows: 

• The introductory chapters provide guidelines on using the MTH$ 
mathematics routines. 

• The MTH$ Reference Section provides detailed reference information on 
each mathematics routine contained in the MTH$ facility of the Run¬ 
Time Library. This information is presented using the documentation 
format described in the Introduction to the VMS Run-Time Library. Routine 
descriptions appear in alphabetical order by routine name. 










Preface 


Associated Documents 

The Run-Time Library routines are documented in a series of reference 
manuals. A general overview of the Run-Time Library and a description 
of how the Run-Time Library routines are accessed is presented in the 
Introduction to the VMS Run-Time Library. Descriptions of the other RTL 
facilities and their corresponding routines and usages are discussed in the 
following books: 

• The VMS RTL DECtalk (DTK$) Manual 

• The VMS RTL Library (LIB$) Manual 

• The VMS RTL General Purpose (OTS$) Manual 

• The VMS RTL Parallel Processing (PPL$) Manual 

• The VMS RTL Screen Management (SMG$) Manual 

• The VMS RTL String Manipulation (STR$) Manual 

The VAX Procedure Calling and Condition Handling Standard, which is 
documented in the Introduction to System Routines, contains useful information 
for anyone who wants to call Run-Time Library routines. 

Applications programmers of any language may refer to the Guide to Creating 
VMS Modular Procedures for the Modular Programming Standard and other* 
guidelines. 

High-level language programmers will find additional information on calling 
Run-Time Library routines in their language reference manual. Additional 
information may also be found in the language user's guide provided with 
your VAX language. 

The Guide to Using VMS Command Procedures may also be useful. 

For a complete list and description of the manuals in the VMS documentation 
set, see the Overview of VMS Documentation. 







Preface 


Conventions 


Convention 

Meaning 

[REjJ 

In examples, a key name (usually abbreviated) 
shown within a box indicates that you press 
a key on the keyboard; in text, a key name is 
not enclosed in a box. In this example, the key 
is the RETURN key. (Note that the RETURN 
key is not usually shown in syntax statements 
or in all examples; however, assume that you 
must press the RETURN key after entering a 
command or responding to a prompt.) 

CTRL/C 

A key combination, shown in uppercase with a 
slash separating two key names, indicates that 
you hold down the first key while you press the 
second key. For example, the key combination 
CTRL/C indicates that you hold down the key 
labeled CTRL while you press the key labeled C. 
In examples, a key combination is enclosed in a 
box. 

$ SHOW TIME 

05-JUN-1988 11:55:22 

In examples, system output (what the system 
displays) is shown in black. User input (what 
you enter) is shown in red. 

$ TYPE MYFILE.DAT 

In examples, a vertical series of periods, or 
ellipsis, means either that not all the data that 
the system would display in response to a 
command is shown or that not all the data a 
user would enter is shown. 

input-file, . . . 

In examples, a horizontal ellipsis indicates 
that additional parameters, values, or other 
information can be entered, that preceding 
items can be repeated one or more times, or 
that optional arguments in a statement have 
been omitted. 

[logical-name] 

Brackets indicate that the enclosed item is 
optional. (Brackets are not, however, optional 
in the syntax of a directory name in a file 
specification or in the syntax of a substring 
specification in an assignment statement.) 

quotation marks 
apostrophes 

The term quotation marks is used to refer 
to double quotation marks ("). The term 
apostrophe (') is used to refer to a single 
quotation mark. 


Other conventions used in the documentation of Run-Time Library routines 
are described in the Introduction to the VMS Run-Time Library. 





















Introduction to MTH$ 


The Run-Time Library mathematics routines may be called to perform a wide 
variety of computations including the following: 

• Complex exponentiation 

• Complex function evaluation 

• Exponentiation 

• Floating-point trigonometric function evaluation 

• Miscellaneous function evaluation 

The OTS$ facility provides additional language-independent arithmetic 
support routines. 

This introduction to Run-Time Library mathematics routines includes 
examples of how to call mathematics routines from BASIC, COBOL, 
FORTRAN, MACRO, PASCAL, and PL/I. 


Entry Point Names 

The names of the mathematics routines are formed by adding the MTH$ 
prefix to the function names. 

When function arguments and returned values are of the same data type, the 
first letter of the name indicates this data type. When function arguments and 
returned values are of different data types, the first letter indicates the data 
type of the returned value, and the second letter indicates the data type of the 
argument(s). 

The letters used as data type prefixes are listed below. 


Letter 

Data Type 

1 

Word 

J 

Longword 

D 

D_floating 

G 

G_floating 

H 

H_floating 

C 

F_floating complex 

CD 

D_floating complex 

CG 

G floating complex 


Generally, F-floating data types have no letter designation. For example, 
MTH$SIN returns an F-floating value of the sine of an F-floating argument 
and MTH$DSIN returns a D-floating value of the sine of a D-floating 
argument. However, in some of the miscellaneous functions, F-floating 
data types are referenced by the letter designation A. 
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Introduction to MTH$ 

1.2 Calling Conventions 


Calling Conventions 

All calls to mathematics routines, as described in the FORMAT section of 
each routine, accept arguments passed by reference. JSB entry points accept 
arguments passed by value. 

All mathematics routines return values in RO or R0/R1 except those routines 
for which the values cannot fit in 64 bits. D-floating complex, G-floating 
complex and H-floating values are data structures which are larger than 64 
bits. Routines that return values which cannot fit in registers R0/R1 return 
their function values into the first argument in the argument list. 

The notation JSB MTH$NAME—Rn, where n is the highest register number 
referenced, indicates that an equivalent JSB entry point is available. No 
registers are saved; only registers R0:Rn are changed. 

Routines with JSB entry points accept a single argument in R0:Rm, where m, 
which is defined below, is dependent on the data type. 


Data Type 

m 

F_floating 

0 

D_floating 

1 

G_floating 

1 

H_floating 

3 


A routine which returns one value returns it to registers R0:Rm. 

When a routine returns two values, for example MTH$SINCOS, the 
first value is returned in R0:Rm and the second value is returned in 
(R <m+l> :R <2*m+l> ). 

Note that for routines that return a single value, n> =m. For routines that 
return two values, n> =2*m + 1. 

All CALL entry points for mathematics routines do the following: 

• Disable floating-point underflow 

• Enable integer overflow 

• Cause no floating-point overflow or other arithmetic traps or faults 

• Preserve all other enabled operations across the CALL 

JSB entry points execute in the context of the caller with the enable operations 
as set by the caller. Since the routines do not cause arithmetic traps or faults, 
their operation is not affected by the setting of the arithmetic trap enables, 
except as noted. 

For more detailed information on CALL and JSB entry points, refer to the 
Introduction to the VMS Run-Time Library. 


Algorithms 

For those mathematics routines that have corresponding algorithms, the 
complete algorithm can be found in the Description section of the routine 
description appearing in the MTH$ Reference Section of this manual. 
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1.4 Condition Handling 


1.4 Condition Handling 

Error conditions are indicated by using the VAX signaling mechanism. The 
VAX signaling mechanism signals all conditions in mathematics routines 
as SEVERE by calling LIB$SIGNAL. When a SEVERE error is signaled, the 
image is caused to exit after printing an error message. A user-established 
condition handler can be written to cause execution to continue at the point 
of the error by returning SS$_CONTINUE. A mathematics routine returns to 
its caller after the contents of R0/R1 have been restored from the mechanism 

argument vector CHF$L_MCH_SAVR0/R1. Thus, the user-established 

handler should correct CHF$L_MCH_SAVR0/R1 to the desired function 
value to be returned to the caller of the mathematics routine. 

D-floating complex, G-floating complex, and H-floating values cannot be 
corrected with a user-established condition handler, because R2/R3 are not 
available in the mechanism argument vector. 

Note that it is more reliable to correct RO and R1 to resemble RO and R1 of 
a double-precision floating-point value. A double-precision floating-point 
value correction works for both single- and double-precision values. If 
the correction is not performed, the floating-point reserved operand -0.0 
is returned. A floating-point reserved operand is a floating-point datum 
with a sign bit of 1 and a biased exponent of zero. Accessing the floating¬ 
point reserved operand will cause a reserved operand fault. See the VMS 
RTL Library (LIB$) Manual for a complete description of how to write user 
condition handlers for SEVERE errors. 

A few mathematics routines signal floating underflow if the calling program 
(JSB or CALL) has enabled floating underflow faults or traps. 

All mathematics routines access input arguments and the real and imaginary 
parts of complex numbers using floating-point instructions. Therefore, a 
reserved operand fault can occur in any mathematics routine. 


1.5 Complex Numbers 

A complex number y is defined as an ordered pair of real numbers r and i, 
where r is the real part and i is the imaginary part of the complex number. 

Y=(r,i) 

VMS supports three floating-point complex types: F-floating complex, D- 
floating complex, and G-floating complex. There is no H-floating complex 
data type. 

Run-Time Library mathematics routines that use complex arguments require 
two x-floating values to be passed by reference for each argument. The 
first x-floating value contains r, the real part of the complex number. The 
second x-floating value contains i, the imaginary part of the complex number. 
Similarly, Run-Time Library mathematics routines that return complex 
function values return two x-floating values. Some Language Independent 
Support (OTS$) routines also calculate complex functions. 

Note that complex functions have no JSB entry points. 
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Introduction to MTH$ 

1.6 Routines Not Documented in the MTH$ Reference Section 


1.6 


Routines Not Documented in the MTH$ Reference Section 


The mathematics routines in Table 1-1 are not found in the reference section 
of this manual. Instead, their entry points and argument information are 
listed in Appendix A of this manual. 

A reserved operand fault can occur for any floating-point input argument 
in any mathematics routine. Other condition values signaled by each 
mathematics routine are indicated in the footnotes. 


Table 1—1 Additional Mathematics Routines 


Entry Point 

Function 

Absolute Value Routines 

MTH$ABS 

F-floating absolute value 

MTH$DABS 

D-floating absolute value 

MTH$GABS 

G-floating absolute value 

MTH$HABS 

H-floating absolute value 1 

MTHSIIABS 

Word absolute value 2 

MTH$JIABS 

Longword absolute value 2 

Bitwise AND Operator Routines 

MTH$IIAND 

Bitwise AND of two word arguments 

MTH$JIAND 

Bitwise AND of two longword arguments 

F-floating Conversion Routines 

MTH$DBLE 

Convert F-floating to D-floating (exact) 

MTH$GDBLE 

Convert F-floating to G-floating (exact) 

MTH$IIFIX 

Convert F-floating to word (truncated) 2 

MTHSJIFIX 

Convert F-floating to longword (truncated) 2 


Returns value to the first argument; value exceeds 64 bits. 
2 Integer overflow exceptions can occur. 
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1.6 Routines Not Documented in the MTH$ Reference Section 


Table 1-1 (Cont.) Additional Mathematics Routines 

Entry Point Function 


Floating-Point Positive Difference Routines 


MTH$DIM 

MTHSDDIM 

MTH$GDIM 

MTH$HDIM 

MTHSIIDIM 

MTH$JIDIM 


Positive difference of two F-floating arguments 3,4 
Positive difference of two D-floating arguments 3,4 
Positive difference of two G-floating arguments 3,4 
Positive difference of two H-floating arguments 1,3,4 
Positive difference of two word arguments 2 
Positive difference of two longword arguments 2 


Bitwise Exclusive OR Operator Routines 


MTHSIIEOR Bitwise exclusive OR of two word arguments 

MTH$JIEOR Bitwise exclusive OR of two longword arguments 


Integer to Floating-point Conversion Routines 


MTH$FLOATI 

MTH$DFLOTI 

MTH$GFLOTI 

MTH$FLOATJ 

MTH$DFLOTJ 

MTH$GFLOTJ 


Convert word to F-floating (exact) 
Convert word to D-floating (exact) 
Convert word to G-floating (exact) 
Convert longword to F-floating (exact) 
Convert word to D-floating (exact) 
Convert longword to G-floating (exact) 


Conversion to Greatest Floating-point Integer Routines 


MTH$FLOOR 

MTHSDFLOOR 

MTH$GFLOOR 

MTH$HFLOOR 


Convert F-floating to greatest F-floating integer 
Convert D-floating to greatest D-floating integer 
Convert G-floating to greatest G-floating integer 
Convert H-floating to greatest H-floating integer 1 


Returns value to the first argument; value exceeds 64 bits. 
2 Integer overflow exceptions can occur. 

3 Floating-point overflow exceptions can occur. 
4 Floating-point underflow exceptions can occur. 
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1.6 Routines Not Documented in the MTH$ Reference Section 


Table 1-1 (Cont.) 

Additional Mathematics Routines 

Entry Point 

Function 

Floating-point Truncation Routines 

MTH$AINT 

Convert F-floating to truncated F-floating 3 

MTH$DINT 

Convert D-floating to truncated D-floating 

MTH$IIDINT 

Convert D-floating to truncated word 2 

MTHSJIDINT 

Convert D-floating to truncated longword 2 

MTH$GINT 

Convert G-floating to truncated G-floating 

MTH$IIGINT 

Convert G-floating to truncated word 2 

MTH$JIGINT 

Convert G-floating to truncated longword 2 

MTH$HINT 

Convert H-floating to truncated H-floating 1,3 

MTH$IIHINT 

Convert H-floating to truncated word 2 

MTH$JIHINT 

Convert H-floating to truncated longword 2 

MTH$IINT 

Convert F-floating to truncated word 2 

MTHSJINT 

Convert F-floating to truncated longword 2 

Bitwise Inclusive OR Operator Routines 

MTH$IIOR 

Bitwise inclusive OR of two word arguments 

MTH$JIOR 

Bitwise inclusive OR of two longword arguments 

Maximum Value Routines 

MTH$AIMAXO 

F-floating maximum of n word arguments 

MTHSAJMAXO 

F-floating maximum of n longword arguments 

MTHSIMAXO 

Word maximum of n word arguments 

MTH$JMAXO 

Longword maximum of n longword arguments 

MTH$AMAX1 

F-floating maximum of n F-floating arguments 2 

MTH$DMAX1 

D-floating maximum of n D-floating arguments 

MTH$GMAX1 

G-floating maximum of n G-floating arguments 

MTH$HMAX1 

H-floating maximum of n H-floating arguments 1 

MTH$IMAX1 

Word maximum of n F-floating arguments 2 

MTH$JMAX1 

Longword maximum of n F-floating arguments 2 


Returns value to the first argument; value exceeds 64 bits, 
integer overflow exceptions can occur. 

3 Floating-point overflow exceptions can occur. 
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1.6 Routines Not Documented in the MTH$ Reference Section 


Table 1 -1 (Cont.) Additional Mathematics Routines 
Entry Point Function 

Minimum Value Routines 


MTH$AIMINO 

MTH$AJMINO 

MTH$IMINO 

MTHSJMINO 

MTH$AMIN1 

MTH$DMIN1 

MTH$GMIN1 

MTH$HMIN1 

MTH$IMIN1 

MTH$JMIN1 


F-floating minimum of n word arguments 
F-floating minimum of n longword arguments 
Word minimum of n word arguments 
Longword minimum of n longword arguments 
F-floating minimum of n F-floating arguments 2 
D-floating minimum of n D-floating arguments 
G-floating minimum of n G-floating arguments 
H-floating minimum of n H-floating arguments' 
Word minimum of n F-floating arguments 2 
Longword minimum of n F-floating arguments 2 


Remainder Routines 


MTHSAMOD 

MTH$DMOD 

MTHSGMOD 

MTHSHMOD 

MTH$IMOD 

MTHSJMOD 


Remainder of two F-floating arguments, arg 1 /arg2 3 
Remainder of two D-floating arguments, arg1/arg2 3 
Remainder of two G-floating arguments, arg1/arg2 3 
Remainder of two H-floating arguments, arg1/arg2 1,3 
Remainder of two word arguments, arg 1 /arg2 5 
Remainder of two longword arguments, arg1/arg2 5 


Floating-point Conversion to Nearest Value Routines 


MTHSANINT 

MTH$DNINT 

MTHSIIDNNT 

MTH$JIDNNT 

MTHSGNINT 

MTHSIIGNNT 

MTHSJIGNNT 

MTHSHNINT 


Convert F-floating to nearest F-floating integer 
Convert D-floating to nearest D-floating integer 3 
Convert D-floating to nearest word integer 
Convert D-floating to nearest longword integer 
Convert G-floating to nearest G-floating integer 3 
Convert G-floating to nearest word integer 2 
Convert G-floating to nearest longword integer 2 
Convert H-floating to nearest H-floating integer 1 


1 Returns value to the first argument; value exceeds 64 bits, 
integer overflow exceptions can occur. 

3 Floating-point overflow exceptions can occur. 
5 Divide-by-zero exceptions can occur. 
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1.6 Routines Not Documented in the MTH$ Reference Section 


Table 1—1 (Cont.) Additional Mathematics Routines 


Entry Point 

Function 

MTH$IIHNNT 

MTH$JIHNNT 

MTH$ININT 

MTH$JNINT 

Convert H-floating to nearest word integer 2 

Convert H-floating to nearest longword integer 2 

Convert F-floating to nearest word integer 2 

Convert F-floating to nearest longword integer 3,6 

Bitwise Complement Operator Routines 

MTH$INOT 

MTH$JNOT 

Bitwise complement of word argument 

Bitwise complement of longword argument 

Floating-point Multiplication Routines 

MTH$DPROD 

MTH$GPROD 

D-floating product of two F-floating arguments 3 

G-floating product of two F-floating arguments 3 

Bitwise Shift Operator Routines 

MTHSIISHFT 

MTHSJISHFT 

Bitwise shift of word 

Bitwise shift of longword 

Floating-point Sign Function Routines 

MTHSSGN 

MTHSSIGN 

MTHSDSIGN 

MTHSGSIGN 

MTH$HSIGN 

F- or D-floating sign function 

F-floating transfer of sign of y to sign of x 

D-floating transfer of sign of y to sign of x 

G-floating transfer of sign of y to sign of x 

H-floating transfer of sign of y to sign of x 1 


Returns value to the first argument; value exceeds 64 bits, 
integer overflow exceptions can occur. 

3 Floating-point overflow exceptions can occur. 

6 Returns contents of RO if a negative argument is input. 
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1.6 Routines Not Documented in the MTH$ Reference Section 


Table 1 -1 (Cont.) Additional Mathematics Routines 

Entry Point Function 

MTH$IISIGN Word transfer of sign of y to sign of x 

MTH$JISIGN Longword transfer of sign of y to sign of x 

Conversion of Double to Single Floating-point Routines 

MTHSSNGL Convert D-floating to F-floating (rounded) 3 

MTHSSNGLG Convert G-floating to F-floating (rounded) 3,4 


3 Floating-point overflow exceptions can occur. 
4 Floating-point underflow exceptions can occur. 


1.7 Examples of Calls to Run-Time Library Mathematics Routines 
1.7.1 BASIC Example 

The following BASIC program uses the H-floating data type. BASIC also 
supports the D-floating, F-floating and G-floating data types, but does not 
support the complex data types. 


10 ! + 

! Sample program to demonstrate a call to MTH$HEXP from BASIC. 

I - 

EXTERNAL SUB MTH$HEXP ( HFLOAT, HFLOAT ) 

DECLARE HFLOAT X.Y ! X and Y are H-floating 

DIGITS$ = '###.#################################' 

X = ’1.2345678901234567891234567892’H 
CALL MTH$HEXP (Y.X) 

A$ = ’MTHSHEXP of ’ + DIGITSS + ’ is ’ + DIGITSS 

PRINT USING A$, X, Y 

END 


The output from this program is as follows: 

MTH$HEXP of 1.234567890123456789123456789200000 
is 3.436893084346008004973301321342110 


1.7.2 COBOL Example 

The following COBOL program uses the F-floating and D-floating data types. 
COBOL does not support the G-floating and H-floating data types or the 
complex data types. 

This COBOL program calls MTH$EXP and MTH$DEXP. 
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IDENTIFICATION DIVISION. 

PROGRAM-ID. FLOATING-POINT. 

* 

* Calls MTH$EXP using a Floating Point data type. 

* Calls MTH$DEXP using a Double Floating Point data type. 

* 

ENVIRONMENT DIVISION. 

DATA DIVISION. 

WORKING-STORAGE SECTION. 

01 FLOAT.PT COMP-1. 

01 ANSWER.F COMP-1. 

01 DOUBLE.PT COMP-2. 

01 ANSWER.D COMP-2. 

PROCEDURE DIVISION. 

PO. 

MOVE 12.34 TO FLOAT.PT. 

MOVE 3.456 TO DOUBLE.PT. 

CALL "MTHSEXP" USING BY REFERENCE FLOAT.PT GIVING ANSWER.F. 
DISPLAY " MTHSEXP of ", FLOAT.PT CONVERSION, " is ", 

ANSWER.F CONVERSION. 


CALL "MTHSDEXP" USING BY REFERENCE DOUBLE.PT GIVING ANSWER.D. 
DISPLAY " MTHSDEXP of ", DOUBLE.PT CONVERSION, " is ", 

ANSWER.D CONVERSION . 


STOP RUN. 


The output from this example program is as follows: 

MTHSEXP of 1.234000E+01 is 2.286620E+05 
MTHSDEXP of 3.456000000000000E+00 is 
3.168996280537917E+01 


1.7.3 FORTRAN Examples 

The first two FORTRAN programs below use the D-floating and H-floating 
data types. The third FORTRAN program below uses the F-floating complex 
data type. FORTRAN supports the four floating data types and the three 
complex data types. 


Q c+ 

C This FORTRAN program computes exp(x) in 
C double precision by using the RTL routine " MTHSDEXP x ". 

C 

C Declare X,Y and MTHSDEXP as double precision values. 

C MTHSDEXP(X) will return a double precision value to variable Y. 

C- 

REAL*8 X,Y,MTHSDEXP 
X = 3.456 
Y = MTHSDEXP(X) 

WRITE(6,1) X,Y 

1 FORMAT( * \’MTHSDEXP(\F20.15,’) IS \F20.15) 

END 

The output generated by this FORTRAN example is as follows: 

MTHSDEXP(3.456000000000000) IS 
31.689962805379165 
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g c+ 

C This FORTRAN program computes exp(x) using 
C the RTL routine MTH$HEXP. MTH$HEXP is CALLed by 
C MTH$HEXP(return_value , argument) 

C 

C Declare X,Y as H-floating point values. 

C Given X MTH$HEXP will return the value of exp(X) in Y by the call 
C CALL MTH$HEXP(Y,X). 

C- 

REAL+16 X,Y 

X = 1.2345678901234567891234567892 
CALL MTH$HEXP(Y,X) 

WRITEC6.1) X,Y 

1 FORMAT(* *, ’MTH$HEXP of \E35.30,* is \E35.30) 

END 


This FORTRAN program generates the following output: 

MTH$HEXP of .123456789012345678912345678920E+01 
is .343689308434600800497330132134E+01 

S C+ 

C This FORTRAN program computes the complex log 

C of x using the RTL routine MTH$CL0G. This program also demonstrates 
C two ways the user can create a complex number. 

C 

C Declare Z,Z_L0G,MTH$CMPLX, and MTHSCLOG as complex values and R and I 
C as real values. MTH$CMPLX takes two real arguments and returns one 
C complex number: Z = MTH$CMPLX(R,I) is a complex number with "real" 

C part R and "imaginary" part I. 

C 

C Given a complex number Z, MTH$CL0G(Z) returns the complex natural 
C logarithm of Z. 

C- 

COMPLEX Z,Z.LOG,MTH$CMPLX,MTH$CL0G 
REAL*4 R,I 
R = 3.142563 
I = 7.4367846 
Z = MTH$CMPLX(R,I) 

C+ 

C Z is a complex number with real part R and imaginary part I. 

C- 

TYPE *, ’ The complex number z is’,z 
C+ 

C Compute the natural logarithm of Z = (2,1). 

C Directly define the complex number Z. 

C- 


Z = (2.0,1.0) 

Z_L0G = MTHSCLOG(Z) 

TYPE *, * The complex log of (2,1) is \Z_L0G 
END 


The output generated by this program is as follows: 


The complex number z is (3.142563,7.436785) 
The complex log of (2,1) is 
(0.8047190,0.4636476) 
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1.7.4 MACRO Examples 

MACRO and BLISS support JSB entry points as well as CALLS and CALLG 
entry points. Both MACRO and BLISS support the four floating data types 
and the three complex data types. 

The MACRO programs below illustrate the use of the CALLS and CALLG 
instructions, as well as JSB entry points. 

Q .TITLE EXAMPLE.JSB 

; + 

; This example calls MTH$DEXP by using a Macro JSB command. 

; The JSB command expects R0/R1 to contain the quadword input value X. 

; The result of the JSB will be located in R0/R1. 

.EXTRN MTH$DEXP_R6 ;MTH$DEXP is an external routine. 

.PSECT DATA, PIC, EXE, NOWRT 
X: .DOUBLE 2.0 ; X is 2.0 

.ENTRY EXAMPLE.JSB, ~M<> 

MOVQ X, RO ; X is in registers RO and R1 

JSB G~MTH$DEXP_R6 ; The result is returned in R0/R1. 

RET 

.END EXAMPLE.JSB 

This MACRO program generates the following output: 

RO <~ 732541EC 
R1 <— ED6EC6A6 

That is, MTH$DEXP(2) is 7.3890560989306502 


2 .TITLE EXAMPLE.CALLG 

; + 

; This example calls MTH$HEXP by using a Macro CALLG command. 

; The CALLG command expects that the address of the return value 
; Y, the address of the input value X, and the argument count 2 be 
; stored in memory: this program stores this information in ARGUMENTS. 
; The result of the CALLG will be located in R0/R1. 


.EXTRN MTH$HEXP 
.PSECT DATA, PIC, 
ARGUMENTS: 

.LONG 2 
.ADDRESS Y, X 


.H.FLOATING 2 
.H.FLOATING 0 
.ENTRY EXAMPLE.G, 
CALLG ARGUMENTS, 
RET 

.END EXAMPLE.G 


; MTH$HEXP is an external routine. 

EXE, WRT 

; The CALLG will use two arguments. 

; The first argument must be the address 
; receiving the computed value, while 
; the second argument is used to 
; compute exp(X). 

; X = 2.0 

; Y is the result, initially set to 0. 

~M<> 

G~MTH$HEXP ; CALLG returns the value to Y. 
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The output generated by this MACRO program is as follows: 


address of Y <-- D8E64003 
<— 4DDA4B8D 
<— 3A3BDCC3 
<— B68BA206 

That is, MTH$HEXP of 2.0 returns 
7.38905609893065022723042746057501 


.TITLE EXAMPLE.CALLS 


This example calls MTHSHEXP by using the Macro CALLS command. 

The CALLS command expects the SP to contain the H-floating address of 
the return value, the address of the input argument X and the argument 
count 2. The result of the CALLS will be located in registers R0-R3. 


Y: 

X: 


.EXTRN MTHSHEXP 

.PSECT DATA, PIC, EXE, 

.H_FL0ATING 0 
.H.FL0ATING 2 
.ENTRY EXAMPLES, ~M<> 
M0VAL X, -(SP) 

MOVAL Y, -(SP) 

CALLS Y, G~MTHSHEXP 

RET 

.END EXAMPLES 


; MTHSHEXP is an external routine. 

WRT 

; Y is the result, initially set to 0. 

; X = 2 

; The address of X is in the SP. 

; The address of Y is in the SP 
; The value is returned to the address of Y. 


The output generated by this program is as follows: 


address of Y <— D8E64003 
<-- 4DDA4B8D 
<— 3A3BDCC3 
<-- B68BA206 

That is, MTHSHEXP of 2.0 returns 
7.38905609893065022723042746057501 


.TITLE C0MPLEX.EX1 

This example calls MTHSCLOG by using a MACRO CALLG command. 

To compute the complex natural logarithm of Z = (2.0,1.0) register 

RO is loaded with 2.0, the real part of Z, and register R1 is loaded 

with 1.0, the imaginary part of Z. The CALLG to MTHSCLOG 

returns the value of the natural logarithm of Z in 

registers RO and Rl. RO gets the real part of Z and R1 

gets the imaginary part. 


ARGS: 


REAL: 

IMAG: 


.EXTRN MTHSCLOG 

.PSECT DATA, PIC, EXE, NOWRT 

.LONG 1 ; The CALLG will use one argument. 

.ADDRESS REAL ; The one argument that the CALLG 

; uses is the address of the argument 


.FLOAT 
.FLOAT 
.ENTRY 
CALLG 


2 ; 

1 ; 

C0MPLEX.EX1, ~M<> 
ARGS, G~MTH$CL0G; 


of MTHSCLOG. 
real part of Z is 2.0 
imaginary part Z is 1.0 

MTHSCLOG return the real part of the 
complex natural logarithm in RO and 


the imaginary part in Rl. 
RET 

.END C0MPLEX_EX1 
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This program generates the following output: 

RO <— 0210404E 

R1 <— 63383FED 

That is. MTH$CL0G(2.0,1.0) is 
(0.8047190,0.4636476) 

g .TITLE C0MPLEX.EX2 

; + 

; This example calls MTH$CL0G by using a MACRO CALLS command. 

; To compute the complex natural logarithm of Z = (2.0,1.0) register 
; R0 is loaded with 2.0, the real part of Z, and register R1 is loaded 
; with 1.0, the imaginary part of Z. The CALLS to MTH$CL0G 
; returns the value of the natural logarithm of Z in registers R0 

; and R1. R0 gets the real part of Z and R1 gets the imaginary 

; part. 



.EXTRN 

MTHSCLOG 



.PSECT 

DATA, PIC, EXE, NOWRT 

REAL: 

.FLOAT 

2 ; 

; real part of Z is 2.0 

IMAG: 

.FLOAT 

1 ; 

; imaginary part Z is 1.0 


.ENTRY 

C0MPLEX.EX2, ~M<> 


MOVAL 

REAL, -(SP) 

SP <-- address of Z. Real part of Z : 




in @(SP) and imaginary part is in 


CALLS 

#1, G"MTH$CL0G 

@(SP)+4. 

MTH$CL0G return the real part of the 




complex natural logarithm in RO and 
the imaginary part in R1. 


RET 




.END 

C0MPLEX_EX2 



This MACRO example program generates the following output: 

R0 <— 0210404E 

R1 <— 63383FED 

That is, MTH$CL0G(2.0,1.0) is 
(0.8047190,0.4636476) 


1.7.5 PASCAL Examples 

The following PASCAL programs use the D-floating and H-floating data 
types. PASCAL also supports the F-floating and G-floating data types. 
PASCAL does not support the complex data types, however. 

Q {+> 

{ Sample program to demonstrate a call to MTHSDEXP from PASCAL. 

{-} 

PROGRAM CALL_MTH$DEXP (OUTPUT); 

OJ 

{ Declare variables used by this program. 

{-} 

VAR 

X : DOUBLE := 3.456; { X,Y are D-floating unless overridden } 

Y : DOUBLE; { with /DOUBLE qualifier on compilation } 
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{+> 

{ Declare the RTL routine used by this program. 

{-> 

[EXTERNAL.ASYNCHRONOUS] FUNCTION MTH$DEXP (VAR value : DOUBLE) : DOUBLE; EXTERN; 


BEGIN 

Y := MTH$DEXP (x); 

WRITELN (’MTH$DEXP of X:5:3. ’ is ’. Y:20:16); 

END. 


The output generated by this PASCAL program is as follows: 

MTH$DEXP of 3.456 is 31.6899656462382318 



B 




{+> 

{ Sample program to demonstrate a call to MTH$HEXP from PASCAL. 

{-> 

PROGRAM CALL_MTH$HEXP (OUTPUT); 

W 

{ Declare variables used by this program. 

{-> 

VAR 

X : QUADRUPLE := 1.2345678901234567891234567892; { X is H-floating > 
Y : QUADRUPLE; { Y is H-floating > 

{+} 

{ Declare the RTL routine used by this program. 

{-} 

[EXTERNAL,ASYNCHRONOUS] PROCEDURE MTHSHEXP (VAR h.exp : QUADRUPLE; 
value : QUADRUPLE); EXTERN; 

BEGIN 

MTH$HEXP (Y.X); 

WRITELN (’MTH$HEXP of X:30:28, ' is Y:35:33); 

END. 


This PASCAL program generates the following output: 

MTHSDEXP of 3.456 is 31.6899656462382318 


1.7. 
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PL/I Examples 

The following PL/I programs use the D-floating and H-floating data types to 
test entry points. PL/I also supports the F-floating and G-floating data types. 
PL/I does not support the complex data types, however. 
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* 

This program tests a MTH$D entry point * 

* 

PROC OPTIONS (MAIN) ; 

DCL (MTH$DEXP) 

ENTRY (FLOAT(53)) RETURNS (FLOAT(53)); 

DCL OPERAND FLOAT(53); 

DCL RESULT FLOAT(53); 

/*** Begin test ***/ 

OPERAND = 3.456; 

RESULT = MTH$DEXP(OPERAND); 

PUT EDIT (’MTH$DEXP of OPERAND, ’ is \ RESULT)(A(12),F(5,3),A(4),F(20,15)); 
END TEST; 


/* 

* 

* 

* 

*/ 

TEST: 


The output generated by this PL/I program is as follows: 

MTHSDEXP of 3.456 is 31.689962805379165 


2 /* 

* * 

* This program tests a MTH$H entry point. * 

* Note that in the PL/I statement below, the /G-float switch * 

* is needed to compile both G- and H-floating point MTH$ routines. */ 


TEST: PROC OPTIONS (MAIN) ; 

DCL (MTH$HEXP) 

ENTRY (FLOAT (113), FLOAT (113)) ; 

DCL OPERAND FLOAT (113); 

DCL RESULT FLOAT (113); 

/*** Begin test ***/ 

OPERAND = 1.234578901234567891234567892; 

CALL MTH$HEXP(RESULT,OPERAND); 

PUT EDIT ('MTH$HEXP of ’, OPERAND, ’ is ’, RESULT) (A(12),F(29,27),A(4),F(29,27)); 
end test; 


To run this program, you must use the following DCL commands: 

$ PLI/G.FLOAT EXAMPLE 
$ LINK EXAMPLE 
$ RUN EXAMPLE 

This program generates the following output: 

MTHSHEXP of 1.234578901234567891234567892 is 
3.436930928565989790506225633 
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$UPDSEC 


SUPDSEC 

Update Section File on Disk 

The Update Section File on Disk service writes all modified pages in an 
active private or global section back into the section file on disk. One or 
more I/O requests are queued, based on the number of pages that have 
been modified. 

FORMAT 

SYS$UPDSEC inadr ,[retadr] ,[acmode] ,[updflg] 

,[efn] ,[iosb] ,[astadr] ,[astprm] 

RETURNS 

VMS usage: cond_value 
type: longword (unsigned) 

access: write only 

mechanism: by value 

Longword condition value. All system services (except $EXIT) return by 
immediate value a condition value in RO. Condition values that this service 
returns are listed under CONDITION VALUES RETURNED. 

ARGUMENTS 

inadr 

VMS usage: address_range 
type: longword (unsigned) 

access: read only 

mechanism: by reference-array reference or descriptor 

Starting and ending virtual addresses of the pages that are to be written to the 
section file if they have been modified. The inadr argument is the address 
of a 2-longword array containing, in order, the starting and ending process 
virtual addresses. Only the virtual page number portion of each virtual 
address is used; the low-order 9 bits are ignored. 

The $UPDSEC service scans pages starting at the address contained in the 
first longword specified by inadr and ending at the address contained in the 
second longword. Within this range, $UPDSEC locates read/write pages that 
have been modified and writes them (contiguously, if possible) to the section 
file on disk. Unmodified pages are also written to disk if they share the same 
cluster with modified pages. 

If the starting and ending virtual addresses are the same, a single page is 
written to the section file if the page has been modified. 

The address specified by the second longword may be smaller than the 
address specified by the first longword. 

retadr 

VMS usage: address_range 
type: longword (unsigned) 

access: write only 

mechanism: by reference-array reference or descriptor 

Addresses of the first and last pages that were actually queued for writing, in 
the first $QIO request, back to the section file on disk. The retadr argument 
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DESCRIPTION 


CONDITION 

VALUES 

RETURNED 


SYSTEM SERVICE DESCRIPTIONS 

$UNWIND 


If you do not specify newpc, execution resumes at the location specified by 
the PC in the call frame of the procedure that receives control when the 
unwinding operation is complete. 


The actual unwind is not performed immediately. Rather, the return 

addresses in the call stack are modified so that when the condition handler 
returns, the unwind procedure is called from each frame being unwound. 

During the actual unwinding of the call stack, $UNWIND examines each 
frame in the call stack to see if a condition handler has been declared. If a 
handler has been declared, $UNWIND calls the handler with the condition 
value SS$_UNWIND (indicating that the call stack is being unwound) in 
the condition name argument of the signal array. When you call a condition 
handler with this condition value, that handler can perform any procedure- 
specific clean-up operations that may be required. After the condition handler 
returns, the call frame is removed from the stack. 


SS$_NORMAL 

SS$_ACCVIO 

SS$_INSFRAME 

SS$_NOSIGNAL 

SS$_UNWINDING 


The service completed successfully. 

The call stack is not accessible to the caller. 

This condition is detected when the call stack is 
scanned to modify the return address. 

There are insufficient call frames to unwind to the 
specified depth. 

No signal is currently active for an exception 
condition. 

An unwind operation is already in progress. 
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$UNWIND 

Unwind Call Stack 

The Unwind Call Stack service unwinds the procedure call stack; that is, 
it removes a specified number of call frames from the stack. Optionally, 
it may return control to a new program counter (PC) unwinding the stack 
The SUNWIND service is intended to be called from within a condition¬ 
handling routine. 

FORMAT 

SYS$UNWIND [depadr],[newpc] 

RETURNS 

VMS usage: cond—value 
type: longword (unsigned) 

access: write only 

mechanism: by value 

Longword condition value. All system services (except $EXIT) return by 
immediate value a condition value in RO. Condition values that this service 
returns are listed under CONDITION VALUES RETURNED. 

ARGUMENTS 

depadr 

VMS usage: longword_unsigned 
type: longword (unsigned) 

access: read only 

mechanism: by reference 

Depth to which the procedure call stack is to be unwound. The depadr 
argument is the address of a longword value. The value 0 specifies the call 
frame of the procedure that was executing when the condition occurred (that 
is, no call frames are unwound), 1 specifies the caller of that frame, 2 specifies 
the caller of the caller of that frame, and so on. 

If depadr specifies the value 0, no unwind occurs and $UNWIND returns a 
successful condition value in RO. 

If you do not specify depadr, $UNWIND unwinds the stack to the call frame 
of the procedure that called the procedure that established the condition 
handler that is calling the SUNWIND service. This is the default and the 
normal method of unwinding the procedure call stack. 

newpc 

VMS usage: address 
type: longword (unsigned) 

access: read only 

mechanism: by reference 

New value for the program counter (PC); this value replaces the current value 
of the PC in the call frame of the procedure that receives control when the 
unwinding operation is complete. The newpc argument is a longword value 
containing the address at which execution is to resume. 

Execution resumes at this address when the unwinding operation is complete. 
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CONDITION 

VALUES 

RETURNED 


acmode 

VMS usage: access_mode 
type: longword (unsigned) 

access: read only 

mechanism: by value 

Access mode on behalf of which the request is being made. The acmode 
argument is a longword containing the access mode. The $PSLDEF macro 
defines the symbols for the four access modes. 

The most privileged access mode used is the access mode of the caller. To 
unlock any specified page, the resultant access mode must be equal to or 
more privileged than the access mode of the owner of that page. 


SS$_WASCLR 
SS$_W ASSET 

SSS—ACCVIO 


SS$_NOPR!V 


The service completed successfully. At least one 
of the specified pages was previously unlocked. 

The service completed successfully. All of the 
specified pages were previously locked in the 
working set. 

The inadr argument cannot be read by the 
caller; the retadr argument cannot be written 
by the caller; or a page in the specified range is 
inaccessible or does not exist. 

A page in the specified range is in the system 
address space. 
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This BASIC program prompts for a cosine value and determines the angle 
. . ., Ti_ p~,,p rrortmrAhpH bv this Droeram is as follows. 


$ RUN ACOS 

Cosine value betwen -1 and 
The angle with that cosine 


+1 ? .5 

is 1.0472 radians 


g PROGRAM GETANGLECINPUT,OUTPUT); 

{ This PASCAL program uses MTHSACOS to determine 
{ the angle which has the cosine given as input. 

{-> 

VAR 

COS : REAL; 

FUNCTION MTH$AC0S(COS : REAL) : REAL; 

EXTERN; 


BEGIN 


END. 


WRITE(’Cosine value between -1 and +1: ). 

WRITELNCThe angle with that cosine is MTH$AC0S(C0S). 
’ radians’); 


This PASCAL program prompts for a cosine value and determines the a n gl e 
drat has that cosine. The output generated by this program ts as follows. 


$ RUN ACOS 

Cosine value between -1 and +1. 
The angle with that cosine is 1 


.5 

.04720E+00 radians 
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MTH$xACOSD 


MTHSxACOSD Arc Cosine of Angle Expressed in 

Degrees 

Given the cosine of an angle, the Arc Cosine of Angle Expressed in 
Degrees routine returns that angle (in degrees). 9 tXpressed ,n 


FORMAT 

MTHSACOSD cosine 

MTHSDACOSD cosine 

MTHSGACOSD cosine 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$AC0SD_R4 


MTH$DACOSD_R7 


MTH$GACOSD_R7 

RETURNS 

Each of the above JSB entries accepts as input one of the floating-point types. 

\/MQ 1 iconn- XI__ 


F_floating, D_floating. G_floatinq 
write only 


type: 
access: 

mechanism: by value 
Angle in degrees. The angle returned will have a value in the 

0 < angle < 180 


range 


n H $ ACOSD returns an F-floating number. MTH$DACOSD returns a 
D-floating number. MTH$GACOSD returns a G-floating number 


cosine 

VMS usage: 
type: 
access: 
mechanism: 


floating-point 

F_floating, G_floating. D_floating 
read only 
by reference 


Cosine of the angle whose value (in degrees) is to be return^H • 

argument is the address of a floating-point number that is this cosine. The 
so u e value of cosine must be less than or. equal to 1. For MTHSACOSD 
a°n fl C S P eClhes a " F - floatin S number. For MTH$DACOSD, cosine Specifies' 
num£ 8 " Umber - FO ' M ™* GA COSD. cosine specifies a G-floaSg 
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DESCRIPTION The angle in degrees whose cosine is X is computed as: 


Value of 
Cosine 

Angle Returned 

0 

90 

1 

0 

-1 

180 

0<X< 1 

zATAND(zSQRT(l - X 2 )/X ), where zATAND and zSQRT 
are the Math Library arc tangent and square root routines, 
respectively, of the appropriate data type 

-1 <X <0 

zATAND(zSQRT(l - X 2 )/X) + 180 

1<\X\ 

The error MTH$INVARGMAT is signaled 


The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HACOSD. 


SS$_ROPRAND Reserved operand. The MTH$xACOSD routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit of 
one and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_INVARGMAT Invalid argument. The absolute value of cosine 

is greater than 1. LIBSSIGNAL copies the 
floating-point reserved operand to the mechanism 
argument vector CHF$I—MCH_SAVR0/R1. The 
result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L_MCH_SAVR0/R1. 


EXAMPLE 

PROGRAM ACOSD(INPUT,OUTPUT); 

{+> 

{ This PASCAL program demonstrates the use of 
{ MTH$AC0SD. 

{-> 

FUNCTION MTH$AC0SD(COS : REAL): REAL; EXTERN; 

VAR 

COSINE : REAL; 

RET.STATUS : REAL; 

BEGIN 

COSINE := 0.5; 

RET.STATUS := MTH$AC0SD(COSINE); 

WRITELN(’The angle, in degrees, is: RET.STATUS); 

END. 


CONDITION 

VALUES 

SIGNALED 
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MTH$xACOSD 


The output generated by this PASCAL example program is as follows 
The angle, expressed in degrees, is: 6.00000E+01 
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MTH$xASIN 


MTH$xASIN Arc Sine in Radians 



Given the sine of an angle, the Arc Sine in Radians routine returns that 
angle (in radians). 

FORMAT 

MTHSASIN sine 

MTHSDASIN sine 

MTHSGASIN sine 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$ASIN_R4 

MTH$DASIN_R7 

MTH$GASIN_R7 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 

type: F_floating, D—floating, G_floating 

access: write only 

mechanism: by value 

Angle in radians. The angle returned will have a value in the range 

—7t/2 < angle < 7r/2 

MTH$ASIN returns an F-floating number. MTH$DASIN returns a D-floating 
number. MTH$GASIN returns a G-floating number. 

ARGUMENTS 

sine 

VMS usage: floating-point 

type: F_floating, D—floating, G_floating 

access: read only 

mechanism: by reference 

The sine of the angle whose value (in radians) is to be returned. The sine 
argument is the address of a floating-point number that is this sine. The 
absolute value of sine must be less than or equal to 1 . For MTH$ASIN, sine 
specifies an F-floating number. For MTH$DASIN, sine specifies a D-floating 
number. For MTH$GASIN, sine specifies a G-floating number. 













MTH$xASIN 


DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


The angle in radians whose sine is X is computed as: 


Value of Sine 

Angle Returned 

0 

0 

1 

7t/2 

-1 

—tt/2 

0 < |X| < 1 

zATAN(X/zSQRT(l - X 2 )), where zATAN and zSQRT 
are the Math Library arc tangent and square root routines, 
respectively, of the appropriate data type 

1<|X| 

The error MTH$INVARGMAT is signaled 


The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HASIN. 


SS$_ROPRAND Reserved operand. The MTH$xASIN routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_INVARGMAT Invalid argument. The absolute value of sine 

is greater than 1. LIB$SIGNAL copies the 
floating-point reserved operand to the mechanism 
argument vector CHF$L_MCH_SAVR0/R1. The 
result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L_MCH_SAVR0/R1. 
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MTHSxASIND 


MTH$xASIND Arc Sine in Degrees 



Given the sine of an angle, the Arc Sine in Degrees routine returns that 
angle (in degrees). 

FORMAT 

MTHSASIND sine 

MTHSDASIND sine 

MTHSGASIND sine 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$ASIND_R4 

MTH$DASIND_R7 

MTH$GASIND_R7 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

Angle in degrees. The angle returned will have a value in the range 

—90 < angle < 90 

MTH$ASIND returns an F-floating number. MTH$DASIND returns a D- 
floating number. MTH$GASIND returns a G-floating number. 

ARGUMENTS 

sine 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

Sine of the angle whose value (in degrees) is to be returned. The sine 
argument is the address of a floating-point number that is this sine. The 
absolute value of sine must be less than or equal to 1. For MTH$ASIND, 
sine specifies an F-floating number. For MTH$DASIND, sine specifies a 
D-floating number. For MTH$GASIND, sine specifies a G-floating number. 
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MTH$xASIND 


DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


The angle in degrees whose sine is X is computed as: 


Value of Sine 

Value Returned 

0 

0 

1 

90 

-1 

-90 

0 < \X\ < 1 

zATAND(X/zSQRT( 1 - X 2 )), where zATAND and zSQRT 
are the Math Library arc tangent and square root routines, 
respectively, of the appropriate data type 

1<|X| 

The error MTH$INVARGMAT is signaled 


The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HASIND. 


SS$_ROPRAND Reserved operand. The MTH$xASIND routine 

encountered a floating point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit of 
one and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_INVARGMAT Invalid argument. The absolute value of sine 

is greater than 1. LIB$SIGNAL copies the 
floating-point reserved operand to the mechanism 

argument vector CHF$I_MCH_SAVR0/R1. The 

result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L_MCH_SAVR0/R1. 
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MTH$xATAN 


MTH$xATAN Arc Tangent in Radians 



Given the tangent of an angle, the Arc Tangent in Radians routine returns 
that angle (in radians). 

FORMAT 

MTH$ATAN tangent 

MTH$DATAN tangent 

MTH$GATAN tangent 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$ATAN_R4 

MTH$DATAN_R7 

MTH$GATAN_R7 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F—floating, D_floating, G—floating 

access: write only 

mechanism: by value 

Angle in radians. The angle returned will have a value in the range 

—7t/2 < angle < 7r/2 

MTH$ATAN returns an F-floating number. MTH$DATAN returns a D- 
floating number. MTH$GATAN returns a G-floating number. 

ARGUMENTS 

tangent 

VMS usage: floating-point 
type: F—floating, D_floating, G—floating 

access: read only 

mechanism: by reference 

The tangent of the angle whose value (in radians) is to be returned. The 
tangent argument is the address of a floating-point number that is this 
tangent. For MTH$ATAN, tangent specifies an F-floating number. For 
MTH$DATAN, tangent specifies a D-floating number. For MTH$GATAN, 
tangent specifies a G-floating number. 
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MTH$xATAN 


DESCRIPTION 


CONDITION 

VALUE 

SIGNALED 


In radians, the computation of the arc tangent function is based on the 
following identities: 

arctanpO = X - X 3 /3 + X 5 /5 - X 7 /7 + ... 

arctan(X) = X + X * Q(X 2 ), 

where Q(y) = -Y /3 + Y 2 /5 - Y 3 /7 + ... 

arctan(X) = X * P{X 2 ), 

where P(Y) = 1 - Y/3 + Y 2 /b - Y 3 /7 + ... 

arctan(A’) = ir/2 — arct£in(l/X) 

arctan(A') = arctan(A) + circtan((A’ — A)/( 1 + A* X)) 

for any real A 

The angle in radians whose tangent is X is computed as: 


Value of X 

Angle Returned 

0<X<3/32 

X + X*Q(X 2 ) 

3/32 < X<11 

AT AN (A) + V * (P(V 2 )), where A and ATAN(A) are 
chosen by table lookup and V = (X - A)/(l + A * X) 

11 < X 

tt/2 - W * (P(W 2 )) where W = 1/X 

X<0 

—zAT AN(\X\) 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HATAN. 

SS$_ROPRAND 

Reserved operand. The MTH$xATAN routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$xATAND 


MTH$xATAND Arc Tangent in Degrees 



Given the tangent of an angle, the Arc Tangent in Degrees routine returns 
that angle (in degrees). 

FORMAT 

MTHSATAND tangent 

MTHSDATAND tangent 

MTHSGATAND tangent 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$ATAN D_R4 

MTH$DATAND_R7 

MTH$GATAND_R7 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 

type: F_floating, D_ floating, G_floating 

access: write only 

mechanism: by value 

Angle in degrees. The angle returned will have a value in the range 

—90 < angle < 90 

MTH$ATAND returns an F-floating number. MTH$DATAND returns a 
D-floating number. MTH$GATAND returns a G-floating number. 

ARGUMENTS 

tangent 

VMS usage: floating-point 

type: F_floating, D_ floating, G_floating 

access: read only 

mechanism: by reference 

The tangent of the angle whose value (in degrees) is to be returned. The 
tangent argument is the address of a floating-point number that is this 
tangent. For MTH$ATAND, tangent specifies an F-floating number. For 
MTH$DATAND, tangent specifies a D-floating number. For MTH$GATAND, 
tangent specifies a G-floating number. 
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MTH$xATAND 


DESCRIPTION 


CONDITION 

VALUE 

SIGNALED 


The computation of the arc tangent function is based on the following 
identities: 

arctan(X) = (180/tt) * ( X - X 3 /3 + X 5 /b - X 1 /7 + ...) 
arctan(X) = 64 * X + X * Q(X 2 ), 

where Q{Y) = 180/tt * [(1 - 64 * tt/180)] - Y/ 3 + Y 2 /b - Y 3 / 7 + y 4 /9 
arctan(X) = X * P(X 2 ), 

where P(Y) = 180/tt * [1 - y/3 + Y 2 /5 - Y 3 / 7 + y 4 /9...] 
arctan(X) = 90 — arctan(l/X) 

arctan(X) = arctan(A) + axctan((X - A)/( 1 + A * X)) 

The angle in degrees whose tangent is X is computed as: 


Tangent 

Angle Returned 

X<3/32 

64 * X + X * Q(X 2 ) 

3/32 < X<11 

ATAND(A) + V * P{V 2 ), where A and ATAND(A) are 
chosen by table lookup and V = (X - A)/( 1 + A * X) 

11 < X 

90 — W * ( P(W 2 )), where W = 1/X 

X<0 

-zATAND(IXI) 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HATAND. 

SS$_ROPRAND 

Reserved operand. The MTH$xATAND routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 














MTH$xATAN2 


MTH$xATAN2 Arc Tangent in Radians with Two 



Arguments 

Given sine and cosine, the Arc Tangent in Radians with Two Arguments 
routine returns the angle (in radians) whose tangent is given by the 
quotient of sine and cosine, (sine/cosine). 

FORMAT 

MTH$ATAN2 sine,cosine 

MTH$DATAN2 sine,cosine 

MTH$GATAN2 sine,cosine 

Each of the above formats accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

Angle in radians. MTH$ATAN2 returns an F-floating number. 
MTH$DATAN2 returns a D-floating number. MTH$GATAN2 returns a 
G-floating number. 

ARGUMENTS 

sine 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

Dividend. The sine argument is the address of a floating-point number that 
is this dividend. For MTH$ATAN2, sine specifies an F-floating number. For 
MTH$DATAN2, sine specifies a D-floating number. For MTH$GATAN2, sine 
specifies a G-floating number. 

cosine 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

Divisor. The cosine argument is the address of a floating-point number that 
is this divisor. For MTH$ATAN2, cosine specifies an F-floating number. For 
MTH$DATAN2, cosine specifies a D-floating number. For MTH$GATAN2, 
cosine specifies a G-floating number. 


MTH-17 











MTH$xATAN2 


DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


The angle in radians whose tangent is Y/X is computed as follows, where / is 
defined in the description of MTH$zCOSH. 


Value of Input Arguments 

Angle Returned 


X = 0 or Y/X > 2 (/+1) 

7r/2 * (signY) 


X > 0 and Y/X< 2 (/+1) 

z AT AN(Y/X) 


X < 0 and Y/X< 2 </+1) 

7T * (signY) + zATAN(Y/X) 



The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HATAN2. 


SS$_ROPRAND Reserved operand. The MTH$xATAN2 routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_INVARGMAT Invalid argument. Both cosine and sine are zero. 

LIB$SIGNAL copies the floating-point reserved 
operand to the mechanism argument vector 
CHF$L_MCH_SAVR0/R1. The result is the 
floating-point reserved operand unless you have 

written a condition handler to change CHF$I_ 

MCH_SAVR0/R1. 
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MTH$xAT AN D2 


MTH$xATAND2 Arc Tangent in Degrees with 



Two Arguments 

Given sine and cosine, the Arc Tangent in Degrees with Two Arguments 
routine returns the angle (in degrees) whose tangent is given by the 
quotient of sine and cosine, (sine/cosine). 

FORMAT 

MTH$ATAND2 sine,cosine 

MTH$DATAND2 sine,cosine 

MTH$GATAND2 sine,cosine 

Each of the above formats accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

Angle (in degrees). MTH$ATAND2 returns an F-floating number. 
MTH$DATAND2 returns a D-floating number. MTH$GATAND2 returns 
a G-floating number. 

ARGUMENTS 

sine 

VMS usage: floating-point 
type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

Dividend. The sine argument is the address of a floating-point number that 
is this dividend. For MTH$ATAND2, sine specifies an F-floating number. For 
MTH$DATAND2, sine specifies a D-floating number. For MTH$GATAND2, 
sine specifies a G-floating number. 

cosine 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

Divisor. The cosine argument is the address of a floating-point number 
that is this divisor. For MTH$ATAND2, cosine specifies an F-floating 
number. For MTH$DATAND2, cosine specifies a D-floating number. For 
MTH$GATAND2, cosine specifies a G-floating number. 
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MTH$xATAND2 


DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


The angle in degrees whose tangent is Y/X is computed below and where / is 
defined in the description of MTH$zCOSH. 


Value of Input Arguments 

Angle Returned 

X = 0 or Y/X > 2 (/ * ,) 

90 * ( signY ) 

X>0 and Y/X < 2 ,/+,) 

zATAND(Y/X) 

X < 0 and Y/X < 2 (/+1 > 

180 * (aijny) + zATAND(Y/X) 


The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HATAND2. 


SS$_ROPRAND Reserved operand. The MTH$xATAND2 routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_INVARGMAT Invalid argument. Both cosine and sine are zero. 

LIB$SIGNAL copies the floating-point reserved 
operand to the mechanism argument vector 
CHF$L_MCH_SAVRO/R 1. The result is the 
floating-point reserved operand unless you have 

written a condition handler to change CHF$I_ 

MCH_SAVRO/R 1. 
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MTH$xATANH 


MTH$xATANH Hyperbolic Arc Tangent 



Given the hyperbolic tangent of an angle, the Hyperbolic Arc Tangent 
routine returns the hyperbolic arc tangent of that angle. 

FORMAT 

MTH$ATANH hyperbolic-tangent 

MTH$DATANH hyperbolic-tangent 

MTH$GATANH hyperbolic-tangent 

Each of the above formats accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F—floating, D_floating, G_floating 

access: write only 

mechanism: by value 

The hyperbolic arc tangent of hyperbolic-tangent. MTH$ATANH returns 
an F-floating number. MTH$DATANH returns a D-floating number. 
MTH$GATANH returns a G-floating number. 

ARGUMENTS 

hyperbolic-tangent 

VMS usage: floating-point 
type: F—floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

Hyperbolic tangent of an angle. The hyperbolic-tangent argument is the 
address of a floating-point number that is this hyperbolic tangent. For 
MTH$ATANH, hyperbolic-tangent specifies an F-floating number. For 
MTH$DATANH, hyperbolic-tangent specifies a D-floating number. For 
MTH$GATANH, hyperbolic-tangent specifies a G-floating number. 

DESCRIPTION 

The hyperbolic arc tangent function is computed as follows: 

Value of x Value Returned 

\X\ < 1 zATANH(X) = zLOG((X + l)/(X - l))/2 

|X|>1 An invalid argument is signaled 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HATANH. 
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MTH$xATANH 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_INVARGMAT 


Reserved operand. The MTH$xATANH routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Invalid argument: |X|>1. LIB$SIGNAL copies the 
floating-point reserved operand to the mechanism 

argument vector CHF$I_MCH_SAVR0/R1. The 

result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L_MCH_SAVR0/R1. 
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MTHSCxABS 


MTH$CxABS Complex Absolute Value 



The Complex Absolute Value routine returns the absolute value of a 
complex number (r,i). 

FORMAT 

MTH$CABS complex-number 

MTH$CDABS complex-number 

MTH$CGABS complex-number 

Each of the above three formats accepts as input one of the three floating¬ 
point complex types. 

RETURNS 

VMS usage: floating-point 

type: F_ floating, D_floating, G_ floating 

access: write only 

mechanism: by value 

The absolute value of a complex number. MTH$CABS returns an F-floating 
number. MTH$CDABS returns a D-floating number. MTH$CGABS returns a 
G-floating number. 

ARGUMENT 

complex-number 

VMS usage: complex—number 

type: F_ floating complex, D_floating complex, G_ floating 

complex 
access: read only 

mechanism: by reference 

A complex number (r,i), where r and i are both floating-point complex values. 
The complex-number argument is the address of this complex number. For 
MTH$CABS, complex-number specifies an F-floating complex number. For 
MTH$CDABS, complex-number specifies a D-floating complex number. For 
MTH$CGABS, complex-number specifies a G-floating complex number. 

DESCRIPTION 

The complex absolute value is computed as follows, where MAX is the larger 
of Irl and lil, and MIN is the smaller of Irl and lil. 

result = MAX * SQRT((M IN/M AX) 2 + 1) 
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MTH$CxABS 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_FLOOVEMAT 


Reserved operand. The MTH$CxABS routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point overflow in Math Library when both 
r and i are large. 


EXAMPLES 

Q c+ 

C This FORTRAN example forms the absolute value of an 
C F-floating complex number using MTH$CABS and the 
C FORTRAN random number generator RAN. 

C 

C Declare Z as a complex value and MTH$CABS as a REAL*4 value. 

C MTH$CABS will return the absolute value of Z: Z_NEW = MTH$CABS(Z) 

C- 

COMPLEX Z 

COMPLEX CMPLX 

REAL*4 Z_NEW,MTH$CABS 

INTEGER M 

M = 1234567 

C+ 

C Generate a random complex number with the FORTRAN generic CMPLX. 

C- 

Z = CMPLX(RAN(M),RAN(M)) 

C+ 

C Z is a complex number (r,i) with real part "r" and 
C imaginary part "i". 

C- 

TYPE *, * The complex number z is’.z 

TYPE *, ’It has real part’.REAL(Z),’and imaginary part’,AIMAG(Z) 
TYPE *, ’ ’ 

C+ 

C Compute the complex absolute value of Z. 

C- 


Z.NEW = MTH$CABS(Z) 

TYPE *, * The complex absolute value of’.z,* is’,Z_NEW 
END 


This example uses an F-floating complex number for complex-number. The 
output of this FORTRAN example is as follows: 

The complex number z is (0.8535407,0.2043402) 

It has real part 0.8535407 and imaginary part 0.2043402 

The complex absolute value of (0.8535407,0.2043402) is 0.8776597 
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MTHSCxABS 


2 c+ 

C This FORTRAN example forms the absolute 
C value of a G-floating complex number using 
C MTHSCGABS and the FORTRAN random number 

C generator RAN. 

C 

C Declare Z as a complex value and MTH$CGABS as a 
C REAL*8 value. MTHSCGABS will return the absolute 
C value of Z: Z.NEW = MTHSCGABS(Z). 

C- 

C0MPLEX*16 Z 

REAL*8 Z.NEW.MTHSCGABS 

C+ 

C Generate a random complex number with the FORTRAN 
C generic CMPLX. 

C- 

Z = (12.34567890123,45.536376385345) 

TYPE *, * The complex number z is’.z 
TYPE *, * ’ 

C+ 

C Compute the complex absolute value of Z. 

C- 

Z.NEW = MTHSCGABS(Z) 

TYPE *, ’ The complex absolute value of’.z,* is’,Z.NEW 
END 

This FORTRAN example uses a G-floating complex number for complex- 
number. Because this example uses a G-floating number, it must be compiled 
as follows: 

$ FORTRAN/G MTHEX.FOR 

Notice the difference in the precision of the output generated: 

The complex number z is (12.3456789012300.45.5363763853450) 

The complex absolute value of (12.3456789012300,45.5363763853450) is 
47.1802645376230 
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MTH$CCOS 


MTH$CCOS Cosine of a Complex Number 



(F-floating Value) 

The Cosine of a Complex Number (F-floating Value) routine returns the 
cosine of a complex number as an F-floating value. 

FORMAT 

MTH$CCOS complex-number 

RETURNS 

VMS usage: complex-number 
type: F_floating complex 

access: write only 

mechanism: by value 

The complex cosine of the complex input number. MTH$CCOS returns an 
F-floating complex number. 

ARGUMENTS 

complex-number 

VMS usage: complex—number 
type: F_floating complex 

access: read only 

mechanism: by reference 

A complex number (r,i) where r and i are floating-point numbers. The 
complex-number argument is the address of this complex number. For 
MTH$CCOS, complex-number specifies an F-floating complex number. 

DESCRIPTION 

The complex cosine is calculated as follows: 

result = ( COS(r) * COSH(i), -SIN(r) * SINH(i)) 

The routine descriptions for the D- and G-floating point versions of this 
routine are listed alphabetically under MTH$CxCOS. 

CONDITION 

VALUES 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$CCOS routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_FLOOVEMAT Floating-point overflow in Math Library: the 

absolute value of i is greater than about 88.029 


for F-floating values. 


MTH-26 














MTHSCCOS 


EXAMPLE 

o 

C This FORTRAN example forms the complex 
C cosine of an F-floating complex number using 
C MTH$CC0S and the FORTRAN random number 
C generator RAN. 

C 

C Declare Z and MTHSCCOS as complex values. 

C MTHSCCOS will return the cosine value of 
C Z: Z.NEW = MTHSCCOS(Z) 

C- 

COMPLEX Z,Z.NEW,MTHSCCOS 
COMPLEX CMPLX 
INTEGER M 
M = 1234567 

C+ 

C Generate a random complex number with the 
C FORTRAN generic CMPLX. 

C- 

Z = CMPLX (RAN (M) .RAN (M)) 

C+ 

C Z is a complex number (r,i) with real part "r" and 
C imaginary part "i". 

C- 

TYPE *, * The complex number z is’.z 

TYPE *, * It has real part’.REAL(Z),’and imaginary part’,AIMAG(Z) 
TYPE *, ’ * 

C+ 

C Compute the complex cosine value of Z. 

C- 


Z.NEW = MTHSCCOS(Z) 

TYPE *, ’ The complex cosine value of’.z,’ is’,Z.NEW 
END 

This FORTRAN example demonstrates the use of MTH$CCOS, using the 
MTHSCCOS entry point. The output of this program is as follows: 

The complex number z is (0.8535407,0.2043402) 

It has real part 0.8535407 and imaginary part 0.2043402 

The complex cosine value of (0.8535407,0.2043402) is (0.6710899.-0.1550672) 


MTH-27 








MTH$CxCOS 


MTH$CxCOS Cosine of a Complex Number 



The Cosine of a Complex Number routine returns the cosine of a complex 
number. 

FORMAT 

MTH$CDCOS complex-cosine, complex-number 
MTH$CGCOS complex-cosine ,complex-number 

Each of the above formats accepts as input one of the floating-point complex 
types. 

RETURNS 

None. 

ARGUMENTS 

complex-cosine 

VMS usage: cornplex_number 

type: D_floating complex, G_ floating complex 

access: write only 

mechanism: by reference 

Complex cosine of the complex-number. The complex cosine routines that 
have D-floating and G-floating complex input values write the address of 
the complex cosine into the complex-cosine argument. For MTH$CDCOS, 
the complex-cosine argument specifies a D-floating complex number. For 
MTH$CGCOS, the complex-number argument specifies a G-floating complex 
number. 

complex-number 

VMS usage: complex_number 

type: D_floating complex, G_ floating complex 

access: read only 

mechanism: by reference 

A complex number (r,i) where r and i are floating-point numbers. The 
complex-number argument is the address of this complex number. For 
MTH$CDCOS, complex-number specifies a D-floating complex number. For 
MTH$CGCOS, complex-number specifies a G-floating complex number. 

DESCRIPTION 

The complex cosine is calculated as follows: 

result = ( COS(r) * COSH{i), -SIN{r) * SINH(i)) 
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MTH$CxCOS 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_FLOOVEMAT 


Reserved operand. The MTH$CxCOS routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point overflow in Math Library: the 
absolute value of i is greater than about 88.029 
for F-floating and D-floating values or greater than 
709.089 for G-floating values. 


EXAMPLE 

c+ 

c This FORTRAN example forms the complex 
C cosine of a D-floating complex number using 
C MTH$CDC0S and the FORTRAN random number 

C generator RAN. 

C 

C Declare Z and MTH$CDC0S as complex values. 

C MTH$CDC0S will return the cosine value of 
C Z: Z.NEW = MTH$CDCOS(Z) 

C- 

COMPLEX*16 Z,Z.NEW,MTH$CDC0S 
COMPLEX*16 DCMPLX 
INTEGER M 
M = 1234567 

C+ 

C Generate a random complex number with the 
C FORTRAN generic DCMPLX. 

C- 

Z = DCMPLX(RAN(M),RAN(M)) 

C+ 

C Z is a complex number (r,i) with real part "r" and 
C imaginary part "i". 

C- 

TYPE *, * The complex number z is’.z 
TYPE *, * * 

C+ 

C Compute the complex cosine value of Z. 

C- 


Z.NEW = MTH$CDCOS(Z) 

TYPE *, * The complex cosine value of’.z,* is’,Z.NEW 
END 


MTH-29 









MTHSCxCOS 


This FORTRAN example program demonstrates the use of MTH$CxCOS, 
using the MTH$CDCOS entry point. Notice the high precision of the output 
generated: 

The complex number z is (0.8535407185554504,0.2043401598930359) 

The complex cosine value of (0.8535407185554504,0.2043401598930359) is 
(0.6710899028500762.-0.1550672019621661) 
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MTH$CEXP 


MTH$CEXP Complex Exponential (F-floating 



Value) 

The Complex Exponential (F-floating Value) routine returns the complex 
exponential of a complex number as an F-floating value. 

FORMAT 

MTH$CEXP complex-number 

RETURNS 

VMS usage: complex-number 
type: F_floating complex 

access: write only 

mechanism: by value 

Complex exponential of the complex input number. MTH$CEXP returns an 
F-floating complex number. 

ARGUMENTS 

complex-number 

VMS usage: complex-number 
type: F-floating complex 

access: read only 

mechanism: by reference 

Complex number whose complex exponential is to be returned. This complex 
number has the form (r,i), where r is the real part and i is the imaginary part. 
The complex-number argument is the address of this complex number. For 
MTH$CEXP, complex-number specifies an F-floating number. 

DESCRIPTION 

The complex exponential is computed as follows: 

complex — exponent = ( EXP(r ) * COS(i),EXP{r) * SIN(i)) 

The routine descriptions for the D- and G-floating point versions of this 
routine are listed alphabetically under MTH$CxEXP. 
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MTHSCEXP 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_FLOOVEMAT 


Reserved operand. The MTHSCEXP routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point overflow in Math Library: the 
absolute value of r is greater than about 88.029 
for F-floating values. 


EXAMPLE 

c+ 

C This FORTRAN example forms the complex exponential 
C of an F-floating complex number using MTH$CEXP 
C and the FORTRAN random number generator RAN 
C 

C Declare Z and MTH$CEXP as complex values. MTH$CEXP 
C will return the exponential value of Z: Z_NEW = MTH$CEXP(Z) 

C- 

COMPLEX Z,Z_NEW,MTHSCEXP 
COMPLEX CMPLX 
INTEGER M 
M = 1234567 

C+ 

C Generate a random complex number with the 
C FORTRAN generic CMPLX. 

C- 

Z = CMPLX(RAN(M),RAN(M)) 

C+ 

C Z is a complex number (r,i) with real part "r" 

C and imaginary part "i". 

C- 

TYPE *, * The complex number z is’.z 

TYPE *, ’ It has real part’.REAL(Z),’and imaginary part’.AIMAG(Z) 
TYPE *, ’ ’ 

C+ 

C Compute the complex exponential value of Z 
C- 

Z.NEW = MTHSCEXP(Z) 

TYPE *, ’ The complex exponential value of’.z,’ is’.Z.NEW 
END 


This FORTRAN program demonstrates the use of MTH$CEXP as a function 
call. The output generated by this example is as follows: 

The complex number z is (0.8535407.0.2043402) 

It has real part 0.8535407 and imaginary part 0.2043402 
The complex exponential value of (0.8535407,0.2043402) is 
(2.299097,0.4764476) 
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MTHSCxEXP 


MTH$CxEXP Complex Exponential 



The Complex Exponential routine returns the complex exponential of a 
complex number. 

FORMAT 

MTH$CDEXP complex-exponent ,complex-number 
MTH$CGEXP complex-exponent ,complex-number 

Each of the above formats accepts as input one of the floating-point complex 
types. 

RETURNS 

None. 

ARGUMENTS 

complex-exponent 

VMS usage: complex-number 

type: D—floating complex, G_ floating complex 

access: write only 

mechanism: by reference 

Complex exponential of complex-number. The complex exponential routines 
that have D-floating complex and G-floating complex input values write the 
complex-exponent into this argument. For MTH$CDEXP, complex-exponent 
argument specifies a D-floating complex number. For MTH$CGEXP, 
complex-exponent specifies a G-floating complex number. 

complex-number 

VMS usage: complex-number 

type: D—floating complex, G_ floating complex 

access: read only 

mechanism: by reference 

Complex number whose complex exponential is to be returned. This complex 
number has the form (r,i), where r is the real part and i is the imaginary 
part. The complex-number argument is the address of this complex number. 
For MTH$CDEXP, complex-number specifies a D-floating number. For 
MTH$CGEXP, complex-number specifies a G-floating number. 

DESCRIPTION 

The complex exponential is computed as follows: 

complex — exponent = ( EXP(r ) * COS(i), EXP(r) * SIN(i)) 
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MTH$CxEXP 


CONDITION 

VALUES 

SIGNALED 


SSS—ROPRAND 


MTH$_FLOOVEMAT 


Reserved operand. The MTH$CxEXP routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point overflow in Math Library: the 
absolute value of r is greater than about 88.029 
for D-floating values or greater than about 709.089 
for G-floating values. 


EXAMPLE 

c+ 

C This FORTRAN example forms the complex exponential 
C of a G-floating complex number using MTH$CGEXP 
C and the FORTRAN random number generator RAN. 

C 

C Declare Z and MTH$CGEXP as complex values. 

C MTH$CGEXP will return the exponential value 
C of Z: CALL MTH$CGEXP(Z_NEW,Z) 

C- 

C0MPLEX*16 Z.Z.NEW 
COMPLEX*16 MTH$GCMPLX 
REAL*8 R,I 
INTEGER M 
M = 1234567 

C+ 

C Generate a random complex number with the FORTRAN 
C- generic CMPLX. 

C- 


R = RAN(M) 

I = RAN(M) 

Z = MTH$GCMPLX(R,I) 

TYPE *, * The complex number z is’.z 
TYPE *. * * 


Compute the complex exponential value of Z. 


CALL MTH$CGEXP(Z_NEW,Z) 

TYPE *, 1 The complex exponential value of’,z,’ is’,Z_NEW 
END 


This FORTRAN example demonstrates how to access MTH$CGEXP as a 
procedure call. Because G-floating numbers are used, this program must be 
compiled using the command "FORTRAN/G filename". 
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MTH$CxEXP 


Notice the high precision of the output generated: 

The complex number z is (0.853540718555450,0.204340159893036) 

The complex exponential value of (0.853540718555450,0.204340159893036) is 
(2.29909677719458,0.476447678044977) 
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MTHSCLOG 


MTH$CLOG Complex Natural Logarithm 



(F-floating Value) 

The Complex Natural Logarithm (F-floating Value) routine returns the 
complex natural logarithm of a complex number as an F-floating value. 

FORMAT 

MTH$CLOG complex-number 

RETURNS 

VMS usage: complex-number 
type: F_floating complex 

access: write only 

mechanism: by value 

The complex natural logarithm of a complex number. MTHSCLOG returns 
an F-floating complex number. 

ARGUMENTS 

complex-number 

VMS usage: complex—number 
type: F_floating complex 

access: read only 

mechanism: by reference 

Complex number whose complex natural logarithm is to be returned. This 
complex number has the form (r,i), where r is the real part and i is the 
imaginary part. The complex-number argument is the address of this 
complex number. For MTHSCLOG, complex-number specifies an F-floating 
number. 

DESCRIPTION 

The complex natural logarithm is computed as follows: 

CLOG(x) = ( LOG{CABS(x)),ATAN2(i,r )) 

The routine descriptions for the D- and G-floating point versions of this 
routine are listed alphabetically under MTHSCxLOG. 

CONDITION 

VALUE 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$CLOG routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 


of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTHSCLOG 


EXAMPLE 

Examples of using MTH$CLOG from VAX MACRO (using both the CALLS 
and the CALLG instructions) appear in the introductory section of this 
manual. 
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MTH$CxLOG 


MTH$CxLOG Complex Natural Logarithm 



The Complex Natural Logarithm routine returns the complex natural 
logarithm of a complex number. 

FORMAT 

MTH$CDLOG complex-natural-log f complex-number 
MTH$CGLOG complex-natural-log ,complex-number 

Each of the above formats accepts as input one of the floating-point complex 
types. 

RETURNS 

None. 

ARGUMENTS 

complex-natural-log 

VMS usage: complex-number 

type: D_floating complex, G_floating complex 

access: write only 

mechanism: by reference 

Natural logarithm of the complex number specified by complex-number. 

The complex natural logarithm routines that have D-floating complex 
and G-floating complex input values write the address of the complex 
natural logarithm into complex-natural-log. For MTH$CDLOG, the 
complex-natural-log argument specifies a D-floating complex number. 

For MTH$CGLOG, the complex-natural-log argument specifies a G-floating 
complex number. 

complex-number 

VMS usage: complex-number 

type: D-floating complex, G-floating complex 

access: read only 

mechanism: by reference 

Complex number whose complex natural logarithm is to be returned. This 
complex number has the form (r,i), where r is the real part and i is the 
imaginary part. The complex-number argument is the address of this 
complex number. For MTH$CDLOG, complex-number specifies a D-floating 
number. For MTH$CGLOG, complex-number specifies a G-floating number. 

DESCRIPTION 

The complex natural logarithm is computed as follows: 

CLOG(x ) = {LOG{CABS(x)),ATAN2{i,r)) 
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MTH$CxLOG 


SS$_ROPRAND Reserved operand. The MTH$CxLOG routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 


EXAMPLE 

C+ 

C This FORTRAN example forms the complex logarithm 
C of a D-floating complex number by using MTH$CDL0G 

C and the FORTRAN random number generator RAN. 

C 

C Declare Z and MTH$CDL0G as complex values. Then MTH$CDL0G 
c will return the logarithm of Z: CALL MTH$CDLOG(Z_NEW,Z). 

C 

C Declare Z,Z_L0G, and MTH$DCMPLX as complex values, 

C and R and I as real values. MTH$DCMPLX takes two real 
C arguments and returns one complex number. 

C 

C Given a complex number Z, MTH$CDLOG(Z) returns the 
C complex natural logarithm of Z. 

C- 

COMPLEX*16 Z,Z_NEW,MTH$DCMPLX 
REAL*8 R,I 

R = 3.1425637846746565 FORTRAN/G MTHEX.FOR 
I = 7.43678469887 
Z = MTH$DCMPLX(R,I) 

C+ 

C Z is a complex number (r,i) with real part "r" and imaginary 
C part "i". 

C- 

TYPE *, * The complex number z is’.z 
TYPE *. ’ > 

CALL MTH$CDL0G(Z_NEW,Z) 

TYPE The complex logarithm of’,z,’ is’,Z_NEW 
END 

This FORTRAN example program uses MTH$CDLOG by calling it as a 
procedure. The output generated by this program is as follows: 

The complex number z is (3.142563784674657,7.436784698870000) 

The complex logarithm of (3.142563784674657,7.436784698870000) is 
(2.088587642177504,1.170985519274141) 


CONDITION 

VALUE 

SIGNALED 
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MTHSCMPLX 


MTH$CMPLX Complex Number Made from 



F-floating-Point 

The Complex Number Made from F-floating-Point routine returns a 
complex number from two floating-point input values. 

FORMAT 

MTH$CMPLX real-part, imaginary-part 

RETURNS 

VMS usage: complex_number 
type: F_floating complex 

access: write only 

mechanism: by value 

A complex number. MTH$CMPLX returns an F-floating complex number. 

ARGUMENTS 

real-part 

VMS usage: floating-point 
type: F_floating 

access: read only 

mechanism: by reference 

Real part of a complex number. The real-part argument is the address 
of a floating-point number that contains this real part, r, of (r,i). For 
MTH$CMPLX, real-part specifies an F-floating number. 

imaginary-part 

VMS usage: floating-point 
type: F_floating 

access: read only 

mechanism: by reference 

Imaginary part of a complex number. The imag-parg argument is the address 
of a floating-point number that contains this imaginary part, i, of (r,i). For 
MTH$CMPLX, imaginary-part specifies an F-floating number. 

DESCRIPTION 

The MTH$CMPLX routines return a complex number from two F-floating 
input values. The routine descriptions for the D- and G-floating point 
versions of this routine are listed alphabetically under MTH$xCMPLX. 
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MTH$CMPLX 


SS$_ROPRAND Reserved operand. The MTH$CMPLX routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 


EXAMPLE 

c+ 

C This FORTRAN example forms two F-floating 
C point complex numbers using MTH$CMPLX 
C and the FORTRAN random number generator RAN. 

C 

C Declare Z and MTH$CMPLX as complex values, and R 

C and I as real values. MTH$CMPLX takes two real 

C F-floating point values and returns one COMPLEX+8 number. 

C 

C Note, since CMPLX is a generic name in FORTRAN, it would be 
C sufficient to use CMPLX. 

C CMPLX must be declare to be of type C0MPLEX*8. 

C 

C Z = CMPLX(R,I) 

C- 

COMPLEX Z,MTH$CMPLX,CMPLX 
REAL*4 R,I 
INTEGER M 
M = 1234567 
R = RAN(M) 

I = RAN(M) 

Z = MTH$CMPLX(R,I) 

C+ 

C Z is a complex number (r,i) with real part "r" and 
C imaginary part "i". 

C- 

TYPE *, * The two input values are:',R,I 
TYPE *, * The complex number z is’.z 
z = CMPLX(RAN(M),RAN(M)) 

TYPE *, ’ * 

TYPE *, * Using the FORTRAN generic CMPLX with random R and I:* 

TYPE *, * The complex number z is’.z 
END 

This FORTRAN example program demonstrates the use of MTH$CMPLX. 
The output generated by this program is as follows: 

The two input values are: 0.8535407 0.2043402 

The complex number z is (0.8535407,0.2043402) 

Using the FORTRAN generic CMPLX with random R and I: 

The complex number z is (0.5722565,0.1857677) 


CONDITION 

VALUE 

SIGNALED 
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MTH$xCMPLX 


MTH$xCMPLX Complex Number Made from D- 



or G-floating-Point 

The Complex Number Made from D- or G-floating-Point routine returns a 
complex number from two D- or G-floating input values. 

FORMAT 

MTH$DCMPLX complx,ireal-part ,imaginary-part 
MTH$GCMPLX complx,real-part,imaginary-part 

Each of the above formats accepts as input one of floating-point complex 
types. 

RETURNS 

None. 

ARGUMENTS 

complx 

VMS usage: complex_number 

type: D_floating complex, G-floating complex 

access: write only 

mechanism: by reference 

The floating-point complex value of a complex number. The complex 
exponential functions that have D-floating complex and G-floating complex 
input values write the address of this floating-point complex value into 
complx. For MTH$DCMPLX, complx specifies a D-floating complex number. 
For MTH$GCMPLX, complx specifies a G-floating complex number. For 
MTH$CMPLX, complx is not used. 

real-part 

VMS usage: floating-point 
type: D_floating, G_floating 

access: read only 

mechanism: by reference 

Real part of a complex number. The real-part argument is the address 
of a floating-point number that contains this real part, r, of (r,i). 

For MTH$DCMPLX, real-part specifies a D-floating number. For 
MTH$GCMPLX, real-part specifies a G-floating number. 

imaginary-part 

VMS usage: floating-point 
type: D_floating, G_floating 

access: read only 

mechanism: by reference 

Imaginary part of a complex number. The imag-parg argument is the 
address of a floating-point number that contains this imaginary part, i, of 
(r,i). For MTH$DCMPLX, imaginary-part specifies a D-floating number. For 
MTH$GCMPLX, imaginary-part specifies a G-floating number. 
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MTH$xCMPLX 


CONDITION 

VALUE 

SIGNALED 


SS$_ROPRAND 


Reserved operand. The MTH$xCMPLX routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 


EXAMPLE 

c+ 

C This FORTRAN example forms two D-floating 
C point complex numbers using MTH$CMPLX 
C and the FORTRAN random number generator RAN. 

C 

C Declare Z and MTH$DCMPLX as complex values, and R 

C and I as real values. MTH$DCMPLX takes two real 

C D-floating point values and returns one 

C COMPLEX*16 number. 

C 

C- 

C0MPLEX*16 Z 
REAL*8 R,I 
INTEGER M 
M = 1234567 
R = RAN(M) 

I = RAN (M) 

CALL MTH$DCMPLX(Z,R,I) 


C+ 

C Z is a complex number (r.i) with real part "r" and imaginary 
C part "i". 

C- 

TYPE *, * The two input values are:’,R,I 
TYPE *, ’ The complex number z is’.Z 
END 


This FORTRAN example demonstrates how to make a procedure call 
to MTH$DCMPLX. Notice the difference in the precision of the output 
generated. 

The two input values are: 0.8535407185554504 0.2043401598930359 

The complex number z is (0.8535407185554504,0.2043401598930359) 
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MTH$CONJG 


MTH$CONJG Conjugate of a Complex Number 



(F-floating Value) 

The Conjugate of a Complex Number (F-floating Value) routine returns the 
complex conjugate (r,-i) of a complex number (r,i) as an F-floating value. 

FORMAT 

MTH$CONJG complex-number 

RETURNS 

VMS usage: complex_number 
type: F_floating complex 

access: write only 

mechanism: by value 

Complex conjugate of a complex number. MTH$CONJG returns an F-floating 
complex number. 

ARGUMENTS 

complex-number 

VMS usage: complex-number 
type: F—floating complex 

access: read only 

mechanism: by reference 

A complex number (r,i), where r and i are floating-point numbers. The 
complex-number argument is the address of this floating-point complex 
number. For MTH$CONJG, complex-number specifies an F-floating number. 

DESCRIPTION 

The MTH$CONJG routine return the complex conjugate (r,-i) of a complex 
number (r,i) as an F-floating value. The routine descriptions for the D- 
and G-floating point versions of this routine are listed alphabetically under 
MTH$xCONJG. 

CONDITION 

VALUE 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$CONJG routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 


of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$xCONJG 


MTH$xCONJG Conjugate of a Complex Number 



The Conjugate of a Complex Number routine returns the complex 
conjugate (r,-i) of a complex number (r,i). 

FORMAT 

MTH$DCONJG complex-conjugate ,complex-number 
MTH$GCONJG complex-conjugate, complex-number 

Each of the above formats accepts as input one of the floating-point complex 
types. 

RETURNS 

None. 

ARGUMENTS 

complex-conjugate 

VMS usage: complex.number 

type: D.floating complex, G_ floating complex 

access: write only 

mechanism: by reference 

The complex conjugate (r,-i) of the complex number specified by complex- 
number. MTH$DCONJG and MTH$GCONJG write the address of this 
complex conjugate into complex-conjugate. For MTH$DCONJG, the 
complex-conjugate argument specifies the address of a D-floating complex 
number. For MTH$GCONJG, the complex-conjugate argument specifies the 
address of a G-floating complex number. 

complex-number 

VMS usage: complex-number 

type: D.floating complex, G. floating complex 

access: read only 

mechanism: by reference 

A complex number (r,i), where r and i are floating-point numbers. The 
complex-number argument is the address of this floating-point complex 
number. For MTH$DCONJG, complex-number specifies a D-floating 
number. For MTH$GCONJG, complex-number specifies a G-floating 
number. 

CONDITION 

VALUE 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$xCONJG routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 


of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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EXAMPLE 

c+ 

C This FORTRAN example forms the complex conjugate 
C of a G-floating complex number using MTH$GCONJG 
C and the FORTRAN random number generator RAN. 

C 

C Declare Z, Z_NEW, and MTH$GCONJG as a complex values. 

C MTH$GCONJG will return the complex conjugate 
C value of Z: Z_NEW = MTH$GCONJG(Z). 

C- 

C0MPLEX*16 Z,Z_NEW,MTH$GC0NJG 
COMPLEX*16 MTH$GCMPLX 
REAL*8 R,I,MTH$GREAL,MTH$GIMAG 
INTEGER M 
M = 1234567 

C+ 

C Generate a random complex number with the 
C FORTRAN generic CMPLX. 

C- 

R = RAN (M) 

I = RAN (M) 

Z = MTH$GCMPLX(R,I) 

TYPE *, ’ The complex number z is’.z 
TYPE 1,MTH$GREAL(Z),MTH$GIMAG(Z) 

1 FORMAT(’ with real part ’.F20.16,’ and imaginary part’,F20.16) 
TYPE *, * ’ 

C+ 

C Compute the complex absolute value of Z. 

C- 


Z_NEW = MTH$GCONJG(Z) 

TYPE *, ’ The complex conjugate value of’.z,’ is’,Z_NEW 
TYPE 1,MTH$GREAL(Z_NEW),MTH$GIMAG(Z_NEW) 

END 


This FORTRAN example demonstrates how to make a function call to 
MTH$GCONJG. Because G-floating numbers are used, the examples must be 
compiled with the statement "FORTRAN/G filename". 

The output generated by this program is as follows: 

The complex number z is (0.853540718555450,0.204340159893036) 
with real part 0.8535407185554504 
and imaginary part 0.2043401598930359 

The complex conjugate value of 

(0.853540718555450,0.204340159893036) is 
(0.853540718555450,-0.204340159893036) 
with real part 0.8535407185554504 
and imaginary part -0.2043401598930359 
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MTH$xCOS Cosine of Angle Expressed in 



Radians 

The Cosine of Angle Expressed in Radians routine returns the cosine of a 
given angle (in radians). 

FORMAT 

MTH$COS angle-in-radians 

MTH$DCOS angle-in-radians 

MTH$GCOS angle-in-radians 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$COS_R4 

MTH$DCOS_R7 

MTH$GCOS_R7 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F—floating, D_floating, G—floating 

access: write only 

mechanism: by value 

Cosine of the angle. MTH$COS returns an F-floating number. MTH$DCOS 
returns a D-floating number. MTH$GCOS returns a G-floating number. 

ARGUMENTS 

angle-in-radians 

VMS usage: floating-point 
type: F—floating, D_floating, G—floating 

access: read only 

mechanism: by reference 

The angle in radians. The angle-in-radians argument is the address of a 
floating-point number. For MTH$COS, angle-in-radians is an F-floating 
number. For MTH$DCOS, angle-in-radians specifies a D-floating number. 
For MTHSGCOS, angle-in-radians specifies a G-floating number. 

DESCRIPTION 

See the MTH$xSINCOS routine for the algorithm used to compute the cosine. 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HCOS. 
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CONDITION 

VALUE 

SIGNALED 


SS$_ROPRAND 


Reserved operand. The MTH$xCOS procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$xCOSD Cosine of Angle Expressed in 



Degrees 

The Cosine of Angle Expressed in Degrees routine returns the cosine of a 
given angle (in degrees). 

FORMAT 

MTH$COSD angle-in-degrees 

MTH$DCOSD angle-in-degrees 

MTH$GCOSD angle-in-degrees 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$COSD_R4 

MTH$DCOSD_R7 

MTH$GCOSD_R7 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage, floating-point 
type: F—floating, D_floating, G—floating 

access: write only 

mechanism: by value 

Cosine of the angle. MTH$COSD returns an F-floating number. 
MTH$DCOSD returns a D-floating number. MTH$GCOSD returns a G- 
floating number. 

ARGUMENTS 

angle-in-degrees 

VMS usage: floating-point 
type: F—floating, D_floating, G—floating 

access: read only 

mechanism: by reference 

Angle (in degrees). The angle-in-degrees argument is the address of a 
floating-point number. For MTH$COSD, angle-in-degrees specifies an F- 
floating number. For MTH$DCOSD, angle-in-degrees specifies a D-floating 
number. For MTH$GCOSD, angle-in-degrees specifies a G-floating number. 

DESCRIPTION 

See the MTH$SINCOSD routine for the algorithm used to compute the 
cosine. 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HCOSD. 
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CONDITION 

VALUE SS$_ROPRAND 

SIGNALED 


Reserved operand. The MTH$xCOSD procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$xCOSH Hyperbolic Cosine 



The Hyperbolic Cosine routine returns the hyperbolic cosine of the input 
value. 

FORMAT 

MTHSCOSH floating-point-input-value 

MTH$DCOSH floating-point-input-value 

MTH$GCOSH floating-point-input-value 

Each of the above formats accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F—floating, D—floating, G—floating 

access: write only 

mechanism: by value 

The hyperbolic cosine of the input value floating-point-input-value. 
MTH$COSH returns an F-floating number. MTH$DCOSH returns a D- 
floating number. MTH$GCOSH returns a G-floating number. 

ARGUMENTS 

floating-point-input- value 

VMS usage: floating-point 
type: F—floating, D—floating, G—floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address of 
this input value. For MTH$COSH, floating-point-input-value specifies an 
F-floating number. For MTH$DCOSH, floating-point-input-value specifies a 
D-floating number. For MTH$GCOSH, floating-point-input-value specifies 
a G-floating number. 

DESCRIPTION 

Computation of the hyperbolic cosine depends on the magnitude of the 
input argument. The range of the function is partitioned using four data¬ 
type-dependent constants: a(z), b(z), and c(z). The subscript z indicates the 
data type. The constants depend on the number of exponent bits (e) and the 
number of fraction bits (/) associated with the data type (z). 

The values of e and / are: 

z e f 

F 8 24 

D 8 56 

G 11 53 
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The values of the constants in terms of e and / are: 


Variable 

Value 

a(z) 

2<-//2) 

b(z) 

CEILING[ (/ + l)/2 * ln(2)] 

c(z) 

(2 e 1 ) * ln(2) 

Based on the above definitions, zCOSH(X) is computed as follows: 

Value of X 

Value Returned 


IXI < a(z) 

1 

a(z) < IXI < .25 

Computed using a power series expansion in \X\ 2 

.25 < IXI < b(z) 

(zEXP(\X\) + l/zEXP(\X\))/2 

b(z) <1X1 < c(z) 

zEXP(\X\)/2 

c(z) < Ixl 

Overflow occurs 


This routine description for the H-floating point value is listed alphabetically 
under MTH$HCOSH. 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_FLOOVEMAT 


Reserved operand. The MTH$xCOSH procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point overflow in Math Library: the 
absolute value of floating-point-input-value is 
greater than about yyy ; LIBSSIGNAL copies the 
reserved operand to the signal mechanism vector. 
The result is the reserved operand -0.0 unless a 
condition handler changes the signal mechanism 
vector. 

The values of yyy are: 

MTH$C0SH—88.722 
MTH$DCOSH—88.722 
MTH$GCOSH—709.782 


MTH-52 












MTHSCSIN 


MTHSCSIN 

Sine of a Complex Number 
(F-floating Value) 

The Sine of a Complex Number (F-floating Value) routine returns the sine 
of a complex number (r,i) as an F-floating value. 

FORMAT 

MTH$CSIN complex-number 

RETURNS 

VMS usage: complex_nurnber 
type: F_floating complex 

access: write only 

mechanism: by value 

Complex sine of the complex number. MTH$CSIN returns an F-floating 
complex number. 

\ 

ARGUMENTS 

complex-number 

VMS usage: complex-number 
type: F_floating complex 

access: read only 

mechanism: by reference 

A complex number (r,i), where r and i are floating-point numbers. The 
complex-number argument is the address of this complex number. For 
MTH$CSIN, complex-number specifies an F-floating complex number. 

DESCRIPTION 

The complex sine is computed as follows: 

complex — sine = ( SIN(r ) * COSH(i),COS(r) * SINH{i)) 

The routine descriptions for the D- and G-floating point versions of this 
routine are listed alphabetically under MTH$CxSIN. 

CONDITION 

VALUES 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$CSIN procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_FLOOVEMAT Floating-point overflow in Math Library: the 

absolute value of i is greater than about 88.029 


for F-floating values. 
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MTH$CxSIN Sine of a Complex Number 



The Sine of a Complex Number routine returns the sine of a complex 
number (r,i). 

FORMAT 

MTH$CDSIN complex-sine, complex-number 
MTH$CGSIN complex-sine,complex-number 

Each of the above formats accepts as input one of the floating-point complex 
types. 

RETURNS 

None. 

ARGUMENTS 

complex-sine 

VMS usage: complex—number 
type: D_floating complex, G—floating complex 

access: write only 

mechanism: by reference 

Complex sine of the complex number. The complex sine routines with D- 
floating complex and G-floating complex input values write the complex sine 
into this complex-sine argument. For MTH$CDSIN, complex-sine specifies 
a D-floating complex number. For MTH$CGSIN, complex-sine specifies a 
G-floating complex number. 

complex-number 

VMS usage: complex—number 
type: D_floating complex, G-floating complex 

access: read only 

mechanism: by reference 

A complex number (r,i), where r and i are floating-point numbers. The 
complex-number argument is the address of this complex number. For 
MTH$CDSIN, complex-number specifies a D-floating complex number. For 
MTH$CGSIN, complex-number specifies a G-floating complex number. 

DESCRIPTION 

The complex sine is computed as follows: 

complex - sine = ( SIN(r) * COSH(i),COS{r) * SINH(i)) 
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CONDITION 

VALUES 

SIGNALED 


SSS—ROPRAND 


MTH$_FLOOVEMAT 


Reserved operand. The MTH$CxSIN procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point overflow in Math Library: the 
absolute value of i is greater than about 88.029 
for D-floating values or greater than about 709.089 
for G-floating values. 


EXAMPLE 

c+ 

C This FORTRAN example forms the complex 

C sine of a G-floating complex number using 

C MTHSCGSIN and the FORTRAN random number 

C generator RAN. 

C 

C Declare Z and MTHSCGSIN as complex values. 

C MTHSCGSIN will return the sine value 
C of Z: CALL MTHSCGSIN(Z.NEW,Z) 

C- 

C0MPLEX*16 Z,Z.NEW 
COMPLEX*16 DCMPLX 
REAL*8 R,I 
INTEGER M 
M = 1234567 

C+ 

C Generate a random complex number with the 
C FORTRAN generic DCMPLX. 

C- 

R = RAN(M) 

I = RAN(M) 

Z = DCMPLX(R,I) 

C+ 

C Z is a complex number (r,i) with real part "r" and 
C imaginary part "i". 

C- 

TYPE *, * The complex number z is’.z 
TYPE *, ’ * 

C+ 

C Compute the complex sine value of Z. 

C- 

CALL MTHSCGSIN(Z.NEW,Z) 

TYPE *, * The complex sine value of*,z,* is’,Z.NEW 
END 
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This FORTRAN example demonstrates a procedure call to MTH$CGSIN. 
Because this program uses G-floating numbers, it must be compiled with the 
statement "FORTRAN/G filename". 

The output generated by this program is as follows: 

The complex number z is (0.853540718555450,0.204340159893036) 

The complex sine value of (0.853540718555450,0.204340159893036) is 
(0.769400835484975,0.135253340912255) 
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MTH$CSQRT Complex Square Root (F-floating 

Value) 

The Complex Square Root (F-floating Value) routine returns the complex 
square root of a complex number (r,i). 

FORMAT 

MTH$CSQRT complex-number 

RETURNS 

VMS usage: complex-number 
type: F_floating complex 

access: write only 

mechanism: by value 

The complex square root of complex-number. MTH$CSQRT returns an 
F-floating number. 

ARGUMENTS 

complex-number 

VMS usage: complex-number 
type: F_floating complex 

access: read only 

mechanism: by reference 

Complex number (r,i). The complex-number argument contains the address 
of this complex number. For MTH$CSQRT, complex-number specifies an 
F-floating number. 

DESCRIPTION 

The complex square root is computed as follows. 


First, calculate ROOT and Q using the following equations: 


ROOT = SQRT((ABS{r) + (CABS(r,i))/2) 
Q = *7(2 * ROOT) 

Then, the complex result is given as follows: 


r 

i 

CSQRT((r,i» 

>0 

Any 

(ROOT,Q) 

<0 

>0 

(Q,ROOT) 

<0 

<0 

(-Q.-ROOT) 


The routine descriptions for the D- and G-floating point versions of this 
routine are listed alphabetically under MTH$CxSQRT. 
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CONDITION 

VALUE 

SIGNALED 


SS$_ROPRAND 


Reserved operand. The MTH$CSQRT procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$CxSQRT Complex Square Root 



The Complex Square Root routine returns the complex square root of a 
complex number (r,i). 

FORMAT 

MTH$CDSQRT complex-square-root,complex-number 
MTH$CGSQRT complex-square-root, complex-number 

Each of the above formats accepts as input one of the floating-point complex 
types. 

RETURNS 

None. 

ARGUMENTS 

complex-square-root 

VMS usage: complex—number 

type: D_floating complex, G_ floating complex 

access: write only 

mechanism: by reference 

Complex square root of the complex number specified by complex-number. 
The complex square root routines that have D-floating complex and G- 
floating complex input values write the complex square root into complex- 
square-root. For MTH$CDSQRT, complex-square-root specifies a D-floating 
complex number. For MTH$CGSQRT, complex-square-root specifies a 
G-floating complex number. 

complex-number 

VMS usage: complex-number 

type: D_floating complex, G_floating complex 

access: read only 

mechanism: by reference 

Complex number (r,i). The complex-number argument contains the address 
of this complex number. For MTH$CDSQRT, complex-number specifies a D- 
floating number. For MTH$CGSQRT, complex-number specifies a G-floating 
number. 

DESCRIPTION 

The complex square root is computed as follows. 

First, calculate ROOT and Q using the following equations: 

ROOT = SQRT((ABS{r) + {CABS(r,i))/ 2) 


Q = *7(2 * ROOT) 
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Then, the complex result is given as follows: 


r 

i 

CSQRT((r,i)) 

>0 

any 

(ROOT,Q) 

<0 

>0 

(Q,ROOT) 

<0 

<0 

(-Q,-ROOT) 


CONDITION 

VALUE 

SIGNALED 


SS$_ROPRAND 


Reserved operand. The MTH$CxSQRT procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 


EXAMPLE 

c+ 

C This FORTRAN example forms the complex square 
C root of a D-floating complex number using 
C MTH$CDSQRT and the FORTRAN random number 

C generator RAN. 

C 

C Declare Z and Z_NEW as complex values. MTH$CDSQRT 
C will return the complex square root of 
C Z: CALL MTH$CDSQRT(Z_NEW,Z). 

C- 

C0MPLEX*16 Z,Z_NEW 
COMPLEX*16 DCMPLX 
INTEGER M 
M = 1234567 

C+ 

C Generate a random complex number with the 
C FORTRAN generic CMPLX. 

C- 

Z = DCMPLX(RAN(M),RAN(M)) 

C+ 

C Z is a complex number (r,i) with real part "r" and imaginary 
C part "i". 

C- 

TYPE *, * The complex number z is’,z 
TYPE *, * ’ 

C+ 

C Compute the complex complex square root of Z. 

C- 

CALL MTH$CDSQRT(Z_NEW,Z) 

TYPE *, * The complex square root of’.z,* is*,Z_NEW 
END 
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This FORTRAN example program demonstrates a procedure call to 
MTH$CDSQRT. The output generated by this program is as follows: 

The complex number z is (0.8535407185554504,0.2043401598930359) 

The complex square root of (0.8535407185554504,0.2043401598930359) is 
(0.9303763973040062,0.1098158554350485) 
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MTH$CVT_x_x Convert One Double-Precision 

Value 

The Convert One Double-Precision Value routines convert one double¬ 
precision value to the destination data type and return the result as a 
function value. MTH$CVT_D_G converts a D-floating value to G-floating 
and MTH$CVT_G_D converts a G-floating value to a D-floating value. 

FORMAT 

MTH$CVT_D_G floating-point-input-val 

MTH$CVT_G _D floating-point-input-val 

RETURNS 

VMS usage: floating-point 
type: G_floating, D—floating 

access: write only 

mechanism: by value 

The converted value. MTH$CVT_D_G returns a G-floating value. 
MTH$CVT_G_D returns a D-floating value. 

ARGUMENT 

floating-point-input-val 

VMS usage: floating-point 
type: D—floating, G_floating 

access: read only 

mechanism: by reference 

The input value to be converted. The floating-point-input-val argument is 
the address of this input value. For MTH$CVT_D_G, the floating-point- 
input-val argument specifies a D-floating number. For MTH$CVT_G_D, the 
floating-point-input-val argument specifies a G-floating number. 

DESCRIPTION 

These procedures are designed to function as hardware conversion 
instructions. They fault on reserved operands. If floating-point overflow 
is detected, an error is signaled. If floating-point underflow is detected and 
floating-point underflow is enabled, an error is signaled. 
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CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_FLOOVEMAT 

MTH$_FLOUNDMAT 


Reserved operand. The MTH$CVT_x_x procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point overflow in Math Library. 
Floating-point underflow in Math Library. 
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MTH$CVT_xA_xA 

Convert an Array of Double-Precision Values 



The Convert an Array of Double-Precision Values routines convert a 
contiguous array of double-precision values to the destination data type 
and return the results as an array. MTH$CVT_DA_GA converts D-floating 
values to G-floating and MTH$CVT_GA_DA converts G-floating values to 
D-floating. 

FORMAT 

MTH$CVT_DA_GA floating-point-input-array 

, floating-point-dest-array 
[, array-size] 

MTH$CVT_GA_DA floating-point-input-array 

, floating-point-dest-array 
[, array-size] 

RETURNS 

MTH$CVT_DA_GA and MTH$CVT_GA_DA return the address of the 
output array to the floating-point-dest-array argument. 

ARGUMENTS 

floating-point-input-array 

VMS usage: floating-point 

type: D_floating, G_floating 

access: read only 

mechanism: by reference, array reference 

Input array of values to be converted. The floating-point-input-array 
argument is the address of an array of floating-point numbers. For 
MTH$CVT_DA_GA, floating-point-input-array specifies an array of D- 
floating numbers. For MTH$CVT_GA_DA, floating-point-input-array 
specifies an array of a G-floating numbers. 

floating-point-dest-array 

VMS usage: floating-point 

type: G_floating, D_floating 

access: write only 

mechanism: by reference, array reference 

Output array of converted values. The floating-point-dest-array argument 
is the address of an array of floating-point numbers. For MTH$CVT_DA_ 
GA, floating-point-dest-array specifies an array of G-floating numbers. 

For MTH$CVT_GA_DA, floating-point-dest-array specifies an array of 
D-floating numbers. 
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array-size 

VMS usage: longword—signed 
type: longword (signed) 

access: read only 

mechanism: by reference 

Number of array elements to be converted. The default value is 1. The 
array-size argument is the address of a longword containing this number of 
elements. 

DESCRIPTION 

These procedures are designed to function as hardware conversion 
instructions. They fault on reserved operands. If floating-point overflow 
is detected, an error is signaled. If floating-point underflow is detected and 
floating-point underflow is enabled, an error is signaled. 

CONDITION 

VALUES 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$CVT_xA_xA 

procedure encountered a floating-point reserved 
operand due to incorrect user input. A floating¬ 
point reserved operand is a floating-point datum 
with a sign bit of 1 and a biased exponent of zero. 
Floating-point reserved operands are reserved for 
future use by DIGITAL. 

MTH$_FLOOVEMAT Floating-point overflow in Math Library. 

MTH$_FLOUNDMAT Floating-point underflow in Math Library. 
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MTH$xEXP Exponential 



The Exponential routine returns the exponential of the input value. 

FORMAT 

MTH$EXP floating-point-input-value 

MTH$DEXP floating-point-input-value 

MTH$GEXP floating-point-input-value 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$EXP_R4 

MTH$DEXP_R6 

MTH$GEXP_R6 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

The exponential of floating-point-input-value. MTH$EXP returns an F- 
floating number. MTH$DEXP returns a D-floating number. MTH$GEXP 
returns a G-floating number. 

ARGUMENTS 

floating-point-input- value 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address of 
a floating-point number. For MTH$EXP, floating-point-input-value specifies 
an F-floating number. For MTH$DEXP, floating-point-input-value specifies 
a D-floating number. For MTH$GEXP, floating-point-input-value specifies a 
G-floating number. 
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MTH$xEXP 


DESCRIPTION The exponential of x is computed as: 


Value of x 

Value Returned 

X > c(z) 

Overflow occurs 

X< - c(z) 

0 

|X| < 2" (/+1) 

1 

Otherwise 

2 Y *2 U * 2 W 


where: 

Y = INTEGERS * ln2(E)) 

V = FRAC(x * ln2(E)) * 16 
U = INTEGER^)/ 16 

W = FRAC{V)/16 

2 W = polynomial approximation of degree 4,8, or 8 for z = F, D, or G. 

See also the section on the hyperbolic cosine for definitions of f and c(z). 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HEXP. 


CONDITION 

VALUES 

SIGNALED 


SSS—ROPRAND 


MTH$_FLOOVEMAT 


Reserved operand. The MTH$xEXP routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point overflow in Math Library: floating- 
point-input-value is greater than yyy; LIB$SIGNAL 
copies the reserved operand to the signal 
mechanism vector. The result is the reserved 
operand -0.0 unless a condition handler changes 
the signal mechanism vector. 

The values of yyy are approximately: 

MTHSEXP—88.029 

MTH$DEXP—88.029 

MTH$GEXP—709.089 
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MTHSxEXP 


MTH$_FLOUNDMAT Floating-point underflow in Math Library: floating- 

point-input-value is less than or equal to yyy 
and the caller (CALL or JSB) has set hardware 
floating-point underflow enable. The result is set 
to 0.0. If the caller has not enabled floating-point 
underflow (the default), a result of 0.0 is returned 
but no error is signaled. 

The values of yyy are approximately: 

MTH$EXP-88.722 

MTH$DEXP-88.722 

MTH$GEXP-709.774 


EXAMPLE 


IDENTIFICATION DIVISION. 

PROGRAM-ID. FLOATING-POINT. 

* 

* Calls MTH$EXP using a Floating Point data type. 

* Calls MTH$DEXP using a Double Floating Point data type. 

* 

ENVIRONMENT DIVISION. 

DATA DIVISION. 

WORKING-STORAGE SECTION. 

01 FLOAT.PT COMP-1. 

01 ANSWER.F COMP-1. 

01 DOUBLE.PT COMP-2. 

01 ANSWER.D COMP-2. 

PROCEDURE DIVISION. 

PO. 

MOVE 12.34 TO FLOAT.PT. 

MOVE 3.456 TO DOUBLE.PT. 

CALL "MTH$EXP" USING BY REFERENCE FLOAT.PT GIVING ANSWER.F. 
DISPLAY " MTH$EXP of ", FLOAT.PT CONVERSION, " is ", 

ANSWER.F CONVERSION. 


CALL "MTH$DEXP" USING BY REFERENCE DOUBLE.PT GIVING ANSWER.D. 
DISPLAY " MTH$DEXP of ", DOUBLE.PT CONVERSION, " is ", 

ANSWER.D CONVERSION . 


STOP RUN. 


This sample program demonstrates calls to MTH$EXP and MTH$DEXP from 
COBOL. 

The output generated by this program is as follows: 

MTH$EXP of 1.234000E+01 is 2.286620E+05 
MTH$DEXP of 3.456000000000000E+00 is 
3.168996280537917E+01 
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MTH$HACOS 


MTH$HACOS Arc Cosine of Angle Expressed in 



Radians (H-floating Value) 

Given the cosine of an angle, the Arc Cosine of Angle Expressed in 
Radians (H-floating Value) routine returns that angle (in radians) in H- 
floating-point precision. 

FORMAT 

MTH$HACOS h-radians ,cosine 

jsb entries 

MTH$HACOS_R8 

RETURNS 

None. 

ARGUMENTS 

h-radians 

VMS usage: floating-point 
type: H—floating 

access: write only 

mechanism: by reference 

Angle (in radians) whose cosine is specified by cosine. The h-radians 
argument is the address of an H-floating number that is this angle. 
MTH$HACOS writes the address of the angle into h-radians. 

cosine 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The cosine of the angle whose value (in radians) is to be returned. The cosine 
argument is the address of a floating-point number that is this cosine. The 
absolute value of cosine must be less than or equal to 1. For MTH$HACOS, 
cosine specifies an H-floating number. 
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MTH$HACOS 


DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


The angle in radians whose cosine is X is computed as: 


Value of 
Cosine 

Value Returned 

0 

7t/2 

1 

0 

-1 

7T 

0<X< 1 

zATAN(zSQRT ( 1 - X 2 )/X), where zATAN and zSQRT are the 
Math Library arc tangent and square root routines, respectively, 
of the appropriate data type 

-1 <X<0 

zATAN(zSQRT{ 1 - X 2 )/X) + tt 

1<|X| 

The error MTH$INVARGMAT is signaled 


SS$_ROPRAND Reserved operand. The MTH$xACOS routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit of 
one and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_INVARGMAT Invalid argument. The absolute value of cosine 

is greater than 1. LIBSSIGNAL copies the 
floating-point reserved operand to the mechanism 

argument vector CHF$I_MCH_SAVR0/R1. The 

result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L_MCH_SAVR0/R1. 


MTH-70 












MTH$HACOSD 


MTH$HACOSD Arc Cosine of Angle Expressed in 

Degrees (H-Floating Value) 


Given the cosine of an angle, the Arc Cosine of Angle Expressed in 
Degrees (H-Floating Value) routine returns that angle (in degrees) as an 
H-floating value. 

FORMAT 

MTH$HACOSD h-degrees,cosine 

jsb entries 

MTH$HACOSD_R8 

RETURNS 

None. 


ARGUMENTS 


h-degrees 

VMS usage: floating-point 
type: H—floating 

access: write only 

mechanism: by reference 

Angle (in degrees) whose cosine is specified by cosine. The h-degrees 
argument is the address of an H-floating number that is this angle. 
MTH$HACOSD writes the address of the angle into h-degrees. 

cosine 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

Cosine of the angle whose value (in degrees) is to be returned. The cosine 
argument is the address of a floating-point number that is this cosine. The 
absolute value of cosine must be less than or equal to 1. For MTH$HACOSD, 
cosine specifies an H-floating number. 
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MTH$HACOSD 


DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


The angle in degrees whose cosine is X is computed as: 


Value of 
Cosine 

Angle Returned 

0 

90 

1 

0 

-1 

180 

0<X< 1 

zATAND(zSQRT(l - X 2 )/X), where zATAND and zSQRT 
are the Math Library arc tangent and square root routines, 
respectively, of the appropriate data type 

-1 <X<0 

zATAND(zSQRT(l - X 2 )/X) + 180 

1<\X\ 

The error MTH$INVARGMAT is signaled 


SSS—ROPRAND Reserved operand. The MTH$xACOSD routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit of 
one and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_INVARGMAT Invalid argument. The absolute value of cosine 

is greater than 1. LIB$SIGNAL copies the 
floating-point reserved operand to the mechanism 

argument vector CHF$I_MCH_SAVR0/R1. The 

result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L _MCH_S A VRO/R1. 
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MTH$HASIN 


MTH$HASIN Arc Sine in Radians (H-floating 

Value) 


Given the sine of an angle, the Arc Sine in Radians (H-floating Value) 
routine returns that angle (in radians) as an H-floating value. 

FORMAT 

MTH$HASIN h-radians,sine 

jsb entries 

MTH$HASIN_R8 

RETURNS 

None. 


ARGUMENTS 


h-radians 

VMS usage: 
type: 
access: 
mechanism: 


floating-point 
H-floating 
write only 
by reference 


Angle (in radians) whose sine is specified by sine. The h-radians argument is 
the address of an H-floating number that is this angle. MTH$HASIN writes 
the address of the angle into h-radians. 


sine 

VMS usage: 
type: 
access: 
mechanism: 


floating-point 
H_floating 
read only 
by reference 


The sine of the angle whose value (in radians) is to be returned. The sine 
argument is the address of a floating-point number that is this sine. The 
absolute value of sine must be less than or equal to 1. For MTH$HASIN, 
sine specifies an H-floating number. 


MTH-73 











MTH$HASIN 


DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


The angle in radians whose sine is X is computed as: 


Value of 

Sine 

Angle Returned 

0 

0 

1 

7t/2 

-1 

— 7t/2 

0 < |X| < 1 

zATAN(X/zSQRT( 1 - X 2 )), where zATAN and zSQRT are the 
Math Library arc tangent and square root routines, respectively, 
of the appropriate data type 

1<|X| 

The error MTH$INVARGMAT is signaled 


SS$_ROPRAND 


MTH$_INVARGMAT 


Reserved operand. The MTH$xASIN routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Invalid argument. The absolute value of sine 
is greater than 1. LIBSSIGNAL copies the 
floating-point reserved operand to the mechanism 

argument vector CHF$I_MCH_SAVR0/R1. The 

result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L_MCH_SAVR0/R1. 


I 
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MTHSHASIND 


MTHSHASIND Arc Sine in Degrees (H-Floating 

Value) 

Given the sine of an angle, the Arc Sine in Degrees (H-Floating Value) 
routine returns that angle (in degrees) as an H-floating value. 

FORMAT 

MTHSHASIND h-degrees, sine 

jsb entries 

MTH$HASIND_R8 

RETURNS 

None. 

ARGUMENTS 

h-degrees 

VMS usage: floating-point 
type: H—floating 

access: write only 

mechanism: by reference 

Angle (in degrees) whose sine is specified by sine. The h-degrees argument 
is the address of an H-floating number that is this angle. MTH$HASIND 
writes the address of the angle into h-degrees. 

sine 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

Sine of the angle whose value (in degrees) is to be returned. The sine 
argument is the address of a floating-point number that is this sine. The 
absolute value of sine must be less than or equal to 1. For MTHSHASIND, 
sine specifies an H-floating number. 
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MTH$HASIND 


DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


The angle in degrees whose sine is X is computed as: 


Value of 

Sine 

Value Returned 

0 

0 

1 

90 

-1 

-90 

0 < |X| < 1 

zATAND(X/zSQRT( 1 - X 2 )), where zATAND and zSQRT 
are the Math Library arc tangent and square root routines, 
respectively, of the appropriate data type 

1< \X\ 

The error MTH$INVARGMAT is signaled 


SS$_ROPRAND Reserved operand. The MTH$xASIND routine 

encountered a floating point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit of 
one and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_INVARGMAT Invalid argument. The absolute value of sine 

is greater than 1. LIB$SIGNAL copies the 
floating-point reserved operand to the mechanism 
argument vector CHF$L_MCH_SAVR0/R1. The 
result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L_MCH_SAVR0/R1. 
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MTH$HATAN 


MTH$HATAN Arc Tangent in Radians (H-floating 

Value) 

Given the tangent of an angle, the Arc Tangent in Radians (H-floating 

Value) routine returns that angle (in radians) as an H-floating value. 

FORMAT 

MTHSHATAN h-radians , tangent 

jsb entries 

MTH$HATAN_R8 

RETURNS 

None. 

ARGUMENTS 

h-radians 

VMS usage: floating-point 
type: H_floating 

access: write only 

mechanism: by reference 

Angle (in radians) whose tangent is specified by tangent. The h-radians 
argument is the address of an H-floating number that is this angle. 
MTH$HATAN writes the address of the angle into h-radians. 

tangent 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The tangent of the angle whose value (in radians) is to be returned. The 
tangent argument is the address of a floating-point number that is this 
tangent. For MTH$HATAN, tangent specifies an H-floating number. 

DESCRIPTION 

In radians, the computation of the arc tangent function is based on the 

following identities: 

arctan(X) — X — X 3 /3 - 1 - X 5 /5 — X 1 /7 + ... 

arctan(X) = X + X * Q(X 2 ), 

where Q{Y) = -Y/ 3 + Y 2 /5 - Y 3 /7 + ... 

arctan(X) = X * P(X 2 ), 

where P{Y) = 1 -Y/3 + Y 2 /5 -Y /7 + ... 

arctan(X) = tt /2 - arctan(l/X) 

arctan(X) = arctan(>l) -1- arctan((X - A)/{ 1 + A * X)) 

for any real A 
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MTH$HATAN 


CONDITION 

VALUE 

SIGNALED 


The angle in radians whose tangent is X is computed as: 


Value of X 

Angle Returned 

0<X<3/32 

X + X*Q(X 2 ) 

3/32 < X<11 

AT AN (A) + V * (P(V 2 )), where A and ATAN(A) are 
chosen by table lookup and V = (X - A)/{I + A * X) 

11 <X 

n/2 -W * (P(W 2 )) where W = 1/X 

X<0 

-zATAN(\X\) 


SS$_ROPRAND Reserved operand. The MTH$xATAN routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$HATAND 


MTH$HATAND Arc Tangent in Degrees 



(H-floating Value) 

Given the tangent of an angle, the Arc Tangent in Degrees (H-floating 
Value) routine returns that angle (in degrees) as an H-floating point value. 

FORMAT 

MTH$HATAND h-degrees,tangent 

jsb entries 

MTH$HATAND_R8 

RETURNS 

None. 

ARGUMENTS 

h-degrees 

VMS usage: floating-point 
type: H-floating 

access: write only 

mechanism: by reference 

Angle (in degrees) whose tangent is specified by tangent. The h-degrees 
argument is the address of an H-floating number that is this angle. 
MTH$HATAND writes the address of the angle into h-degrees. 

tangent 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The tangent of the angle whose value (in degrees) is to be returned. The 
tangent argument is the address of a floating-point number that is this 
tangent. For MTH$HATAND, tangent specifies an H-floating number. 

DESCRIPTION 

The computation of the arc tangent function is based on the following 
identities: 

arctan(X) = 180/7T * ( X — X 3 /3 + X 5 /5 — X 7 /7 + ...) 

arctan(X) = 64 * X -1- X * Q(X 2 ), 

where Q{Y) = 180/7T * [(1 - 64 * 7t/180) - Y/ 3+ 

y 2 / 5- y 3 /7 + r 4 /M 

arctan(X) = X * P(X 2 ), 

where P{Y) = 180/tt * [1 - Y/ 3 + Y 2 /5 - Y 3 / 7+ 
y 4 /9...] 

arctan(X) = 90 — arctan(l/X) 

arctan(X) = arctan(A) + arctan((X — A)/( 1 + A* X)) 
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MTH$HATAND 


The angle in degrees whose tangent is X is computed as: 


Tangent 

Angle Returned 

X<3/32 

64*X + X*Q(X 2 ) 

3/32 < X<11 

ATAND(A) + V * P(V 2 ), where A and ATAND(A) are 
chosen by table lookup and V = (X - A)/{\ + A * X) 

11 <X 

90 - W * (P(W 2 )), where W = 1/X 

X<0 

—zAT AND(\X\) 


CONDITION 

VALUE 

SIGNALED 


SS$_ROPRAND 


Reserved operand. The MTH$xATAND routine 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$HATAN2 


MTH$HATAN2 Arc Tangent in Radians 

(H-floating Value) with Two 
Arguments 


Given sine and cosine, the Arc Tangent in Radians (H-floating Value) with 
Two Arguments routine returns the angle (in radians) as an H-floating value 
whose tangent is given by the quotient of sine and cosine, (sine/cosine). 

FORMAT 

MTH$HATAN2 h-radians ,sine,cosine 

RETURNS 

None. 

ARGUMENTS 

h-radians 

VMS usage: floating-point 
type. H-floating 

access: write only 

mechanism: by reference 

Angle (in radians) whose tangent is specified by (sine/cosine). The h- 
radians argument is the address of an H-floating number that is this angle. 
MTH$HATAN2 writes the address of the angle into h-radians. 


sine 

VMS usage: floating-point 
type. H-floating 

access: read only 

mechanism: by reference 

Dividend. The sine argument is the address of a floating-point number that 
is this dividend. For MTH$HATAN2, sine specifies an H-floating number. 


cosine 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

Divisor. The cosine argument is the address of a floating-point number that 
is this divisor. For MTH$HATAN2, cosine specifies an H-floating number. 
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MTH$HATAN2 


DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


The angle in radians whose tangent is Y/X is computed as follows, where /is 
defined in the description of MTH$zCOSH. 

Value of Input Arguments 

Angle Returned 

X = 0 or Y/X > 2 ,/+1) 

X > 0 and Y/X< 2 (/+1) 

X <0 and Y/X< 2 </+ ” 

rr/2 * ( signY ) 
zATAN(Y/X) 

7T * (signY) + zATAN(Y/X) 

. 



SS$_ROPRAND Reserved operand. The MTH$HATAN2 routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_INVARGMAT Invalid argument. Both cosine and sine are zero. 

LIBSSIGNAL copies the floating-point reserved 
operand to the mechanism argument vector 
CHF$L_MCH_SAVR0/R1. The result is the 
floating-point reserved operand unless you have 

written a condition handler to change CHF$I_ 

MCH_SAVR0/R1. 
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MTH$HATAND2 


MTH$HATAND2 Arc Tangent in Degrees 

(H-floating Value) with Two 
Arguments 



Given sine and cosine, MTH$xHTAND2 returns the angle (in degrees) 
whose tangent is given by the quotient of sine and cosine, (sine/cosine). 

FORMAT 

MTH$HATAND2 h-degrees ,sine ,cosine 

RETURNS 

None. 

ARGUMENTS 

h-degrees 

VMS usage: floating-point 
type: H_floating 

access: write only 

mechanism: by reference 

Angle (in degrees) whose tangent is specified by (sine/cosine). The h- 
degrees argument is the address of an H-floating number that is this angle. 
MTH$HATAND2 writes the address of the angle into h-degrees. 

sine 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

Dividend. The sine argument is the address of a floating-point number that 
is this dividend. For MTH$HATAND2, sine specifies an H-floating number. 


cosine 

VMS usage: 
type: 
access: 
mechanism: 


floating-point 
H_floating 
read only 
by reference 


Divisor. The cosine argument is the address of a floating-point number that 
is this divisor. For MTH$HATAND2, cosine specifies an H-floating number. 
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MTH$HATAND2 


DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


MTH-84 


The angle in degrees whose tangent is Y/X is computed below. The value of 
/ is defined in the description of MTH$zCOSH. 


Value of Input Arguments Angle Returned 


X = 0 or Y/X > 2 (/t " 

X > 0 and Y/X< 2 (/+1) 

X < 0 and Y/X< 2 ,/+1) 

90 * ( signY ) 
zATANDiY/X) 

180 * (signY) + zATAND(Y/X) 




SS$_ROPRAND Reserved operand. The MTH$HATAND2 routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_JNVARGMAT Invalid argument. Both cosine and sine are zero. 

LIB$SIGNAL copies the floating-point reserved 
operand to the mechanism argument vector 
CHF$L _MCH _SAVRO/R1. The result is the 
floating-point reserved operand unless you have 

written a condition handler to change CHF$I_ 

MCH—SAVR0/R1. 












MTH$HATANH 


MTH$HATANH Hyperbolic Arc Tangent 



(H-floating Value) 

Given the hyperbolic tangent of an angle, the Hyperbolic Arc Tangent (H- 
floating Value) routine returns the hyperbolic arc tangent (as an H-floating 
value) of that angle. 

FORMAT 

MTH$HATANH h-atanh ,hyperbolic-tangent 

RETURNS 

None. 

ARGUMENTS 

h-atanh 

VMS usage: floating-point 
type: H —floating 

access: write only 

mechanism: by reference 

Hyperbolic arc tangent of the hyperbolic tangent specified by hyperbolic- 
tangent. The h-atanh argument is the address of an H-floating number that 
is this hyperbolic arc tangent. MTH$HATANH writes the address of the 
hyperbolic arc tangent into h-atanh. 

hyperbolic-tangent 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

Hyperbolic tangent of an angle. The hyperbolic-tangent argument is the 
address of a floating-point number that is this hyperbolic tangent. For 
MTH$HATANH, hyperbolic-tangent specifies an H-floating number. 

DESCRIPTION 

The hyperbolic arc tangent function is computed as follows: 

Value of x Value Returned 

\X\ < 1 zATANH(X) = zLOG((X + 1)/(X - l))/2 

|X|>1 An invalid argument is signaled 


MTH-85 













MTH$HATANH 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND Reserved operand. The MTH$xATANH routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_INVARGMAT Invalid argument: IXI > 1. LIB$SIGNAL copies the 

floating-point reserved operand to the mechanism 

argument vector CHF$I_MCH_SAVR0/R1. The 

result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L _MCH_S A VRO/R1. 
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MTH$HCOS 


MTH$HCOS Cosine of Angle Expressed in 



Radians (H-floating Value) 

The Cosine of Angle Expressed in Radians (H-floating Value) routine 
returns the cosine of a given angle (in radians) as an H-floating value. 

FORMAT 

MTH$HCOS h-cosine ,angle-in-radians 

jsb entries 

MTH$HCOS_R5 

RETURNS 

None. 

ARGUMENTS 

h-cosine 

VMS usage: floating-point 
type: H_floating 

access: write only 

mechanism: by reference 

Cosine of the angle specified by angle-in-radians. The h-cosine argument is 
the address of an H-floating number that is this cosine. MTH$HCOS writes 
the address of the cosine into h-cosine. 

angle-in-radians 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The angle in radians. The angle-in-radians argument is the address of 
a floating-point number. For MTH$HCOS ( angle-in-radians specifies an 
H-floating number. 

DESCRIPTION 

See the MTH$xSINCOS routine for the algorithm used to compute the cosine. 

CONDITION 

VALUE 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$HCOS procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 


reserved operands are reserved for future use by 
DIGITAL. 
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MTH$HCOSD Cosine of Angle Expressed in 



Degrees (H-floating Value) 

The Cosine of Angle Expressed in Degrees (H-floating Value) routine 
returns the cosine of a given angle (in degrees) as an H-floating value. 

FORMAT 

MTH$HCOSD h-cosine ,angle-in-degrees 

jsb entries 

MTH$HCOSD_R5 

RETURNS 

None. 

ARGUMENTS 

h-cosine 

VMS usage: floating-point 
type: H—floating 

access: write only 

mechanism: by reference 

Cosine of the angle specified by angle-in-degrees. The h-cosine argument is 
the address of an H-floating number that is this cosine. MTH$HCOSD writes 
this cosine into h-cosine. 

angle-in-degrees 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

Angle (in degrees). The angle-in-degrees argument is the address of a 
floating-point number. For MTH$HCOSD, angle-in-degrees specifies an 
H-floating number. 

DESCRIPTION 

See the MTH$SINCOSD routine for the algorithm used to compute the 
cosine. 

CONDITION 

VALUE 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$HCOSD procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 


of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$HCOSH Hyperbolic Cosine (H-floating 



Value) 

The Hyperbolic Cosine routine returns the hyperbolic cosine of the input 
value as an H-floating value. 

FORMAT 

MTH$HCOSH h-cosh ,floating-point-input-value 

RETURNS 

None. 

ARGUMENTS 

h-cosh 

VMS usage: floating-point 
type: H_floating 

access: write only 

mechanism: by reference 

Hyperbolic cosine of the input value specified by floating-point-input-value. 
The h-cosh argument is the address of an H-floating number that is this 
hyperbolic cosine. MTH$HCOSH writes the address of the hyperbolic cosine 
into h-cosh. 

floating-point-input- value 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address of 
this input value. For MTH$HCOSH, floating-point-input-value specifies an 
H-floating number. 

DESCRIPTION 

Computation of the hyperbolic cosine depends on the magnitude of the 
input argument. The range of the function is partitioned using four data¬ 
type-dependent constants: a(z), b(z), and c(z). The subscript z indicates the 
data type. The constants depend on the number of exponent bits (e) and the 
number of fraction bits (/) associated with the data type (z). 

The values of e and / are as follows: 

e = 15 


/ = 113 
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The values of the constants in terms of e and / are: 


Variable 

Value 

a(z) 

2-// 2 

b(z) 

(/ + l)/2 * ln(2) 

c(z) 

2'- 1 * ln(2) 

Based on the above definitions, zCOSH(X) is computed as follows: 

Value of X 

Value Returned 


i*i < <*(*) 

1 

a(*)<|X| < .25 

Computed using a power series expansion in \X\ 2 

.25<|X| < b(z) 

(zEXP(\X\) + l/zEXP(\X\))/2 

b(z)<\X\ < c(z) 

zEXP{\X\)/2 

c(z)<\X\ 

Overflow occurs 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_FLOOVEMAT 


Reserved operand. The MTH$HCOSH procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point overflow in Math Library: the 
absolute value of floating-point-input-value is 
greater than about yyy ; LIB$SIGNAL copies the 
reserved operand to the signal mechanism vector. 
The result is the reserved operand -0.0 unless a 
condition handler changes the signal mechanism 
vector. The value of yyy is 11356.523. 
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MTH$HEXP Exponential (H-floating Value) 



The Exponential routine returns the exponential of the input value as an 
H-floating value. 

FORMAT 

MTHSHEXP h-exp,floating-point-input-value 

jsb entries 

MTH$HEXP_R6 

RETURNS 

None. 

ARGUMENTS 

h-exp 

VMS usage: floating-point 
type: H—floating 

access: write only 

mechanism: by reference 

Exponential of the input value specified by floating-point-input-value. 

The h-exp argument is the address of an H-floating number that is this 
exponential. MTH$HEXP writes the address of the exponential into h-exp. 

floating-point-input- value 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address 
of a floating-point number. For MTH$HEXP, floating-point-input-value 
specifies an H-floating number. 

DESCRIPTION 

The exponential of x is computed as: 

Value of x Value Returned 

x > c(z) Overflow occurs 

x< — c(z) 0 

|z| < 2- (/+11 1 

Otherwise 2 Y * 2 U * 2 W 

where: 

Y = INTEGERS * ln2{E)) 

V = FRAC(x * ln2(E)) * 16 

U = INTEGER(V)/16 
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W = FRAO(V)/ 16 

2 W = polynomial approximation of degree 14 for z = H. 

See also the section on the hyperbolic cosine for definitions of f and c(z). 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND Reserved operand. The MTH$xEXP routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_FLOOVEMAT Floating-point overflow in Math Library: floating¬ 

point-input-value is greater than yyy; LIB$SIGNAL 
copies the reserved operand to the signal 
mechanism vector. The result is the reserved 
operand -0.0 unless a condition handler changes 
the signal mechanism vector. The value of yyy is 
approximately 11355.830 for MTH$HEXP. 

MTH$_FLOUNDMAT Floating-point underflow in Math Library: floating- 

point-input-value is less than or equal to yyy and 
the caller (CALL or JSB) has set hardware floating¬ 
point underflow enable. The result is set to 0.0. If 
the caller has not enabled floating-point underflow 
(the default), a result of 0.0 is returned but no 
error is signaled. The value of yyy is approximately 
-11356.523 for MTH$HEXP. 
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MTH$HLOG Natural Logarithm (H-floating Value) 



The Natural Logarithm (H-floating Value) routine returns the natural (base 
e) logarithm of the input argument as an H-floating value. 

FORMAT 

MTH$HLOG h-natlog , floating-point-input-value 

jsb entries 

MTH$HLOG_R8 

RETURNS 

None. 

ARGUMENTS 

h-natlog 

VMS usage: floating-point 
type: H—floating 

access: write only 

mechanism: by reference 

Natural logarithm of floating-point-input-value. The h-natlog argument 
is the address of an H-floating number that is this natural logarithm. 
MTH$HLOG writes the address of this natural logarithm into h-natlog. 

floating-point-input- value 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address of 
a floating-point number that is this value. For MTH$HLOG, floating-point- 
input-value specifies an H-floating number. 

DESCRIPTION 

Computation of the natural logarithm routine is based on the following: 

1 ln(X * Y) = ln(X) + ln(Y) 

2 ln(l + X) = X-X 2 /2 + X 3 /3-X 4 /4... 

for IXI < 1 

3 lnpf) = ln(A) + 2 * (Y + V 3 /3 + V 5 /5 + V 7 /7...) 

where V = {X - A)/(X + A), A> 0, 
and p(y) = 2 * {1 + y/3 + y 2 / 5...) 

For x = 2 n * /, where n is an integer and f is in the interval of 0.5 to 1, 
define the following quantities: 

If n> 1, then N = n — 1 and F — 2f 

If n< 0, then N — n and F = / 

* 
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From (1) above it follows that: 

4 ln(X) = N * ln(2) + ln(F) 

Based on the above relationships, zLOG is computed as follows: 

1 If |F — 1| < 2 -5 , 

zLOG(X) = N * zLOG{2) + W+ W * p{W), 
where W = F-l. 

2 Otherwise, 

zLOG{X) = N * zLOG{2) + zLOG(A) + V * p(F 2 ), 
where V = (F - A)/(F + A) and A and zLOG(A) 
are obtained by table look up. 


SS$_ROPRAND Reserved operand. The MTH$HLOG procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_LOGZERNEG Logarithm of zero or negative value. Argument 

floating-point-input-value is less than or equal 
to 0.0. LIBSSIGNAL copies the floating-point 
reserved operand to the mechanism argument 

vector CHF$I_MCH_SAVR0/R1. The result is the 

floating-point reserved operand unless you have 

written a condition handler to change CHF$I_ 

MCH_SAVR0/R1. 


CONDITION 

VALUES 

SIGNALED 
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MTH$HL0G2 Base 2 Logarithm (H-floating 



Value) 

The Base 2 Logarithm (H-floating Value) routine returns the base 2 
logarithm of the input value specified by floating-point-input-value as 
an H-floating value. 

FORMAT 

MTH$HLOG2 h-log2 ,floating-point-input-value 

RETURNS 

None. 

ARGUMENTS 

h-log2 

VMS usage: floating-point 
type: H_floating 

access: write only 

mechanism: by reference 

Base 2 logarithm of floating-point-input-value. The h-log2 argument is the 
address of an H-floating number that is this base 2 logarithm. MTH$HLOG2 
writes the address of this logarithm into h-log2. 

floating-point-input- value 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address 
of a floating-point number that is this input value. For MTH$HLOG2, 
floating-point-input-value specifies an H-floating number. 

DESCRIPTION 

The base 2 logarithm function is computed as follows: 

zLOG2(X) = zLOG2(E) * zLOG(X) 
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CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_LOGZERNEG 


Reserved operand. The MTH$HLOG2 procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Logarithm of zero or negative value. Argument 
floating-point-input-value is less than or equal 
to 0.0. LIB$SIGNAL copies the floating-point 
reserved operand to the mechanism argument 
vector CHF$L_MCH_SAVR0/R1. The result is the 
floating-point reserved operand unless you have 

written a condition handler to change CHF$I_ 

MCH_SAVR0/R 1. 
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MTH$HLOG10 Common Logarithm (H-floating 



Value) 

The Common Logarithm (H-floating Value) routine returns the common 
(base 10) logarithm of the input argument as an H-floating value. 

FORMAT 

MTH$HLOG 10 h-log 10,floating-point-input-value 

jsb entries 

MTH$HLOG10_R8 

RETURNS 

None. 

ARGUMENTS 

h-log 10 

VMS usage: floating-point 
type: H_floating 

access: write only 

mechanism: by reference 

Common logarithm of the input value specified by floating-point-input- 
value. The h-loglO argument is the address of an H-floating number that is 
this common logarithm. MTHSHLOGIO writes the address of the common 
logarithm into h-loglO. 

floating-point-input- value 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address 
of a floating-point number. For MTH$HLOG10, floating-point-input-value 
specifies an H-floating number. 

DESCRIPTION 

The common logarithm function is computed as follows: 

zLOG10(X) = zLOG10{E) * zLOG{X) 
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CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND Reserved operand. The MTH$HLOG10 procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_LOGZERNEG Logarithm of zero or negative value. Argument 

floating-point-input-value is less than or equal 
to 0.0. LIB$SIGNAL copies the floating-point 
reserved operand to the mechanism argument 
vector CHF$L_MCH_SAVR0/R1. The result is the 
floating-point reserved operand unless you have 

written a condition handler to change CHF$I_ 

MCH_SAVR0/R1. 
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MTHSHSIN 

Sine of Angle Expressed in Radians 
(H-floating Value) 

The Sine of Angle Expressed in Radians (H-floating Value) routine returns 
the sine of a given angle (in radians) as an H-fioating value. 

FORMAT 

MTH$HSIN h-sine ,angle-in-radians 

jsb entries 

MTH$HSIN_R5 

RETURNS 

None. 

ARGUMENTS 

h-sine 

VMS usage: floating-point 
type: H—floating 

access: write only 

mechanism: by reference 

The sine of the angle specified by angle-in-radians. The h-sine argument is 
the address of an H-floating number that is this sine. MTH$HSIN writes the 
address of the sine into h-sine. 

angle-in-radians 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

Angle (in radians). The angle-in-radians argument is the address of a 
floating-point number that is this angle. For MTHSHSIN, angle-in-radians 
specifies an H-floating number. 

DESCRIPTION 

See the MTH$SINCOS routine for the algorithm used to compute this sine. 

CONDITION 

VALUE 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$HSIN procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 


of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$HSIND Sine of Angle Expressed in Degrees 



(H-floating Value) 

The Sine of Angle Expressed in Degrees (H-floating Value) routine returns 
the sine of a given angle (in degrees) as an H-floating value. 

FORMAT 

MTH$HSIND h-sine,angle-in-degrees 

jsb entries 

MTH$HSIND_R5 

RETURNS 

None. 

ARGUMENTS 

h-sine 

VMS usage: floating-point 
type: H—floating 

access: write only 

mechanism: by reference 

Sine of the angle specified by angle-in-degrees. The h-sine argument is the 
address of an H-floating number that is this sine. MTH$HSIND writes the 
address of the angle into h-sine. 

angle-in-degrees 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

Angle (in degrees). The angle-in-degrees argument is the address of a 
floating-point number that is this angle. For MTH$HSIND, angle-in-degrees 
specifies an H-floating number. 

DESCRIPTION 

See MTH$SINCOSD for the algorithm used to compute the sine. 
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CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND Reserved operand. The MTH$HSIND procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased ecponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_FLOUNDMAT Floating-point underflow in Math Library. The 

absolute value of the input angle is less than 
180/7n*2 -m (where m = 16,384 for H-floating). 
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MTH$HSINH Hyperbolic Sine (H-floating Value) 



The Hyperbolic Sine (H-floating Value) routine returns the hyperbolic sine 
of the input value specified by floating-point-input-value as an H-floating 
value. 

FORMAT 

MTH$HSINH h-sinh,floating-point-input-value 

RETURNS 

None. 

ARGUMENTS 

h-sinh 

VMS usage: floating-point 
type: H_floating 

access: write only 

mechanism: by reference 

Hyperbolic sine of the input value specified by floating-point-input-value. 
The h-sinh argument is the address of an H-floating number that is this 
hyperbolic sine. MTH$HSINH writes the address of the hyperbolic sine into 

h-sinh. 

floating-point-input- value 

VMS usage: floating-point 
type: H_f loafing 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address of 
a floating-point number that is this value. For MTH$HSINH, floating-point- 
input-value specifies an H-floating number. 

DESCRIPTION 

Computation of the hyperbolic sine function depends on the magnitude of 
the input argument. The range of the function is partitioned using four data 
type dependent constants: a(z), b(z), and c(z). The subscript z indicates the 
data type. The constants depend on the number of exponent bits ( e ) and the 
number of fraction bits (/) associated with the data type (z). 

The values of e and / are as follows: 


e = 15 
/ = 113 
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The values of the constants in terms of e and / are: 


Variable 

Value 

a(z) 

2<-//2> 

b(z) 

(/ + l)/2 * ln(2) 

c(z) 

2 e 1 *ln(2) 

Based on the above definitions, zSINH(X) is computed as follows: 

Value of X 

Value Returned 


1*1 < «(*) 

X 

a(z)<|X| < 1.0 

zSINH(X) is computed using a power series 
expansion in \X\ 2 

1.0<|X| < b{z) 

(zEXP(X) - zEXP(-X))/ 2 

6(z)<|X| < c(z) 

SIGN(X)*zEXP(\X\)/2 

c(z)<\X\ 

Overflow occurs 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND Reserved operand. The MTH$HSINH procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_FLOOVEMAT Floating-point overflow in Math Library: the 

absolute value of floating-point-input-value 
is greater than yyy. LIB$SIGNAL copies the 
floating-point reserved operand to the mechanism 
argument vector CHF$L_MCH_SAVR0/R1. The 
result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L_MCH_SAVR0/R1. The value of yyy is 
approximately 11356.523. 
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MTH$HSQRT Square Root (H-floating Value) 



The Square Root (H-floating Value) routine returns the square root of the 
input value floating-point-input-value as an H-floating value. 

FORMAT 

MTH$HSQRT h-sqrt , floating-point-input-value 

jsb entries 

MTH$HSQRT_R8 

RETURNS 

None. 

ARGUMENTS 

h-sqrt 

VMS usage: floating-point 
type: H_floating 

access: write only 

mechanism: by reference 

Square root of the input value specified by floating-point-input-value. The 
h-sqrt argument is the address of an H-floating number that is this square 
root. MTH$HSQRT writes the address of the square root into h-sqrt. 

floating-point-input- value 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

Input value. The floating-point-input-value argument is the address of 
a floating-point number that contains this input value. For MTH$HSQRT, 
floating-point-input-value specifies an H-floating number. 

DESCRIPTION 

The square root of X is computed as follows: 

If X < 0, an error is signaled. 

Let X = 2 K *F 

where: 

K is the exponential part of the floating-point data 

F is the fractional part of the floating-point data 

If K is even: 

X = 2 (2 * p) * f . 
zSQRT(X) = 2 P * zSQRT(F), 

1/2 <F < 1, where P = K/2 
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If K is odd: 

X = 2 (2 * Pt1) * F = 2 ,2 * p+2) * (F/ 2), 
zSQRT(X) = 2 (Pt1) * zSQRT{F/ 2), 
l/4<F/2 < 1/2, where p = (K-1 )/2 

Let F' = A * F + B, when K is even: 

A = 0.95F6198 (hex) 

B = 0.6BA5918 (hex) 

Let F' = A* (F/2) + B, when K is odd: 

A = 0.D413CCC (hex) 

B = 0.4C1E248 (hex) 

Let K' = P, when K is even 
Let K' = P+1, when K is odd 

Let V[0] = 2 k> * F' be a straight line approximation within the given interval 
using coefficients A and B which minimize the absolute error at the midpoint 
and endpoint. 

Starting with Y[0], n Newton-Raphson iterations are performed: 

Y[n + 1] = 1/2 * (Y[n] + X/Y\n\) 

where n = 5 for H-floating. 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_SQUROONEG 


Reserved operand. The MTH$HSQRT procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Square root of negative number. Argument 
floating-point-input-value is less than 0.0. 
LIBSSIGNAL copies the floating-point reserved 
operand to the mechanism argument vector 
CHF$L_MCH_SAVR0/R1. The result is the 
floating-point reserved operand unless you have 
written a condition handler to change CHF$I— 
MCH_SAVR0/R1. 
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MTH$HTAN Tangent of Angle Expressed in 



Radians (H-floating Value) 

The Tangent of Angle Expressed in Radians (H-floating Value) routine 
returns the tangent of a given angle (in radians) as an H-floating value. 

FORMAT 

MTH$HTAN h-tan,angle-in-radians 

jsb entries 

MTH$HTAN_R5 

RETURNS 

None. 

ARGUMENTS 

h-tan 

VMS usage: floating-point 
type: H_floating 

access: write only 

mechanism: by reference 

Tangent of the angle specified by angle-in-radians. The h-tan argument is 
the address of an H-floating number that is this tangent. MTH$HTAN writes 
the address of the tangent into h-tan. 

angle-in-radians 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The input angle (in radians). The angle-in-radians argument is the address of 
a floating-point number that is this angle. For MTH$HTAN, angle-in-radians 
specifies an H-floating number. 

DESCRIPTION 

When the input argument is expressed in radians, the tangent function is 
computed as follows: 

1 If |X| < 2then zTAN{X) = X (see the section on MTH$zCOSH 
for the definition of f) 

2 Otherwise, call MTH$zSINCOS to obtain zSIN(X) and zCOS(X); then 

a. If zCOS(X) = 0, signal overflow 

b. Otherwise, zTAN(X) = zSIN{X)/zCOS(X) 
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CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND Reserved operand. The MTH$HTAN procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_FLOOVEMAT Floating-point overflow in math library. 
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MTH$HTAND Tangent of Angle Expressed in 



Degrees (H-floating Value) 

The Tangent of Angle Expressed in Degrees (H-floating Value) routine 
returns the tangent of a given angle (in degrees) as an H-floating value. 

FORMAT 

MTH$HTAND h-tan ,angle-in-degrees 

jsb entries 

MTH$HTAND_R5 

RETURNS 

None. 

ARGUMENTS 

h-tan 

VMS usage: floating-point 
type: H_floating 

access: write only 

mechanism: by reference 

Tangent of the angle specified by angle-in-degrees. The h-tan argument 
is the address of an H-floating number that is this tangent. MTH$HTAND 
writes the address of the tangent into h-tan. 

angle-in-degrees 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The input angle (in degrees). The angle-in-degrees argument is the address 
of a floating-point number which is this angle. For MTH$HTAND, angle-in- 
degrees specifies an H-floating number. 

DESCRIPTION 

When the input argument is expressed in degrees, the tangent function is 
computed as follows: 

« 

1 If \X\ < (180/7t) 2 /and underflow signaling is enabled, underflow 
is signaled (see the section on MTH$zCOSH for the definition of e). 

2 Otherwise, if |X| < (180 /tt) * 2then zTAND(X) = (tt/180) * X. See 
the description of MTH$zCOSH for the definition of /. 

3 Otherwise, call MTH$zSINCOSD to obtain zSIND(X) and zCOSD(X). 

a. Then, if zCOSD(X) = 0, signal overflow 

b. Else, zTAND(X) = zSIND(X)/zCOSD{X) 
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MTH$HTAND 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND Reserved operand. The MTH$HTAND procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_FLOOVEMAT Floating-point overflow in math library. 
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MTH$HTANH Compute the Hyperbolic Tangent 



(H-floating Value) 

! 

The Compute the Hyperbolic Tangent (H-floating Value) routine returns the 
hyperbolic tangent of the input value as an H-floating value. 

FORMAT 

MTH$HTANH h-tanh, floating-point-input-value 

RETURNS 

None. 

ARGUMENTS 

h-tanh 

VMS usage: floating-point 
type: H—floating 

access: write only 

mechanism: by reference 

Hyperbolic tangent of the value specified by floating-point-input-value. The 
h-tanh argument is the address of a H-floating number that is this hyperbolic 
tangent. MTH$HTANH writes the address of the hyperbolic tangent into 

h-tanh. 

floating-point-input- value 

VMS usage: floating-point 
type: H_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address of 
a floating-point number that contains this input value. For MTH$HTANH, 
floating-point-input-value specifies an H-floating number. 

DESCRIPTION 

For MTH$HTANH, the hyperbolic tangent of X is computed using a value of 
56 for g and a value of 40 for h. The hyperbolic tangent of X is computed as 
follows: 

Value of x Hyperbolic Tangent Returned 

|X|<2~ 9 X 

2~° < |X|<0.25 zSINH(X)/zCOSH(X) 

0.25 < |X| < h (zEXP( 2 * X) - l)/(zEXP(2 * X) + 1) 

/i<|X| sign(X) * 1 
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MTH$HTANH 


CONDITION 

VALUE 

SIGNALED 


SS$_ROPRAND Reserved operand. The MTHSHTANH procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$xlMAG Imaginary Part of a Complex 



Number 

The Imaginary Part of a Complex Number routine returns the imaginary 
part of a complex number. 

FORMAT 

MTH$AIMAG complex-number 

MTH$DIMAG complex-number 

MTH$GIMAG complex-number 

Each of the above three formats corresponds to one of the three floating-point 
complex types. 

RETURNS 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

Imaginary part of the input complex-number. MTH$AIMAG returns an F- 
floating number. MTH$DIMAG returns a D-floating number. MTH$GIMAG 
returns a G-floating number. 

ARGUMENT 

complex-number 

VMS usage: complex-number 

type: F_floating complex, D_floating complex, 

G_ floating complex 
access: read only 

mechanism: by reference 

The input complex number. The complex-number argument is the address 
of this floating-point complex number. For MTH$AIMAG, complex-number 
specifies an F-floating number. For MTH$DIMAG, complex-number specifies 
a D-floating number. For MTH$GIMAG, complex-number specifies a G- 
floating number. 

CONDITION 

VALUE 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$xlMAG routine 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 


of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$xlMAG 


EXAMPLE 

c+ 

C This FORTRAN example forms the imaginary part of 
C a G-floating complex number using MTH$GIMAG 
C and the FORTRAN random number generator 
C RAN. 

C 

C Declare Z as a complex value and MTH$GIMAG as 
C a REAL*8 value. MTH$GIMAG will return the imaginary 
C part of Z: Z_NEW = MTH$GIMAG(Z). 

C- 

C0MPLEX*16 Z 
COMPLEX*16 DCMPLX 
REAL*8 R,I,MTH$GIMAG 
INTEGER M 
M = 1234567 

C+ 

C Generate a random complex number with the 
C FORTRAN generic CMPLX. 

C- 

R = RAN (M) 

I = RAN(M) 

Z = DCMPLX(R,I) 

C+ 

C Z is a complex number (r,i) with real part "r" and 
C imaginary part "i". 

C- 

TYPE *, 1 The complex number z is’.z 

TYPE *, * It has imaginary part’,MTH$GIMAG(Z) 

END 


This FORTRAN example demonstrates a procedure call to MTH$GIMAG. 
Because this example uses G-floating numbers, it must be compiled with the 
statement "FORTRAN/G filename". 

The output generated by this program is as follows: 

The complex number z is (0.8535407185554504,0.2043401598930359) 

It has imaginary part 0.2043401598930359 
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MTH$xLOG 


MTH$xLOG Natural Logarithm 



The Natural Logarithm routine returns the natural (base e) logarithm of the 
input argument. 

FORMAT 

MTH$ALOG floating-point-input-value 

MTH$DLOG floating-point-input-value 

MTH$GLOG floating-point-input-value 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$ALOG_R5 

MTH$DLOG_R8 

MTH$GLOG_R8 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F_floating, D—floating, G_floating 

access: write only 

mechanism: by value 

The natural logarithm of floating-point-input-value. MTH$ALOG returns an 
F-floating number. MTH$DLOG returns a D-floating number. MTH$GLOG 
returns a G-floating number. 

ARGUMENTS 

floating-point-input- value 

VMS usage: floating-point 

type: F_floating, D—floating, G_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address 
of a floating-point number that is this value. For MTH$ALOG, floating¬ 
point-input-value specifies an F-floating number. For MTH$DLOG, 
floating-point-input-value specifies a D-floating number. For MTH$GLOG, 
floating-point-input-value specifies a G-floating number. 

DESCRIPTION 

Computation of the natural logarithm routine is based on the following: 

1 ln(X * Y) = ln(AT) + ln(y) 

2 ln(l + X) = X — X 2 /2 + X 3 /3 — X*/4 ... 

for IXI < 1 

3 \n(X) = ln(A) + 2*{V + V 3 /3 + V 5 /5 + V 7 /7...) 

= ln(A) + V * p(V 2 ), where V = (X - A)/(X + A), 

A > 0, and p{y) = 2 * (1 + y/3 + y 2 /5...) 
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For x = 2 n * /, where n is an integer and f is in the interval of 0.5 to 1, 
define the following quantities: 

If n> 1, then N = n — 1 and F = 2/ 

If n< 0, then N = n and F = / 

From (1) above it follows that: 

4 ln(X)=iV*ln(2)+ln(F) 

Based on the above relationships, zLOG is computed as follows: 

1 If |F - 1| < 2 -5 , zLOG(X) = N * zLOG{2) + W + W * p{W), 

where W = F-l. 

2 Otherwise, zLOG(X) = N * zLOG{2) + zLOG(A) + V * p(V 2 ), 

where V = (F — A)/(F A) and A and zLOG(A) 
are obtained by table look up. 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HLOG. 

CONDITION 

VALUES 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTH$xLOG procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_LOGZERNEG Logarithm of zero or negative value. Argument 

floating-point-input-value is less than or equal 


to 0.0. LIBSSIGNAL copies the floating-point 
reserved operand to the mechanism argument 

vector CHF$I_MCH—SAVR0/R1. The result is the 

floating-point reserved operand unless you have 
written a condition handler to change CHF$I— 
MCH—SAVR0/R1. 
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MTH$xL0G2 


MTH$xL0G2 Base 2 Logarithm 



The Base 2 Logarithm routine returns the base 2 logarithm of the input 
value specified by floating-point-input-value. 

FORMAT 

MTH$ALOG2 floating-point-input-value 

MTH$DLOG2 floating-point-input-value 

MTH$GLOG2 floating-point-input-value 

Each of the above formats accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 

type: F_ floating, D_floating, G_ floating 

access: write only 

mechanism: by value 

The base 2 logarithm of floating-point-input-value. MTH$ALOG2 
returns an F-floating number. MTH$DLOG2 returns a D-floating number. 
MTH$GLOG2 returns a G-floating number. 

ARGUMENTS 

floating-point-input- value 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address of 
a floating-point number that is this input value. For MTH$ALOG2, floating¬ 
point-input-value specifies an F-floating number. For MTH$DLOG2, 
floating-point-input-value specifies a D-floating number. For MTH$GLOG2, 
floating-point-input-value specifies a G-floating number. 

DESCRIPTION 

The base 2 logarithm function is computed as follows: 

zLOG2{X) = zLOG2(E) * zLOG{X) 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HLOG2. 
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MTH$xL0G2 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_LOGZERNEG 


Reserved operand. The MTH$xLOG2 procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Logarithm of zero or negative value. Argument 
floating-point-input-value is less than or equal 
to 0.0. LIB$SIGNAL copies the floating-point 
reserved operand to the mechanism argument 

vector CHF$I_MCH_SAVR0/R1. The result is the 

floating-point reserved operand unless you have 
written a condition handler to change CHF$L_ 
MCH_SAVR0/R1. 
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MTH$xLOG10 Common Logarithm 



The Common Logarithm routine returns the common (base 10) logarithm 
of the input argument. 

FORMAT 

MTH$ALOG10 floating-point-input-value 

MTH$DLOG10 floating-point-input-value 

MTH$GLOG10 floating-point-input-value 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$ALOG10_R5 

MTH$DLOG10_R8 

MTH$GLOG10_R8 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 

type: F_floating, D—floating, G_floating 

access: write only 

mechanism: by value 

The common logarithm of floating-point-input-value. MTH$ALOG10 
returns an F-floating number. MTH$DLOG10 returns a D-floating number. 
MTH$GLOG10 returns a G-floating number. 

ARGUMENTS 

floating-point-input- value 

VMS usage: floating-point 

type: F_floating, D—floating, G_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address 
of a floating-point number. For MTH$ALOG10, floating-point-input-value 
specifies an F-floating number. For MTH$DLOG10, floating-point-input- 
value specifies a D-floating number. For MTH$GLOG10, floating-point- 
input-value specifies a G-floating number. 

DESCRIPTION 

The common logarithm function is computed as follows: 

zLOG10{X) = zLOG10{E) * zLOG[X) 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HLOG10. 
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MTHSxLOGIO 


CONDITION 

VALUES 

SIGNALED 


SSS—ROPRAND 


MTH$_LOGZERNEG 


Reserved operand. The MTH$xLOG10 procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Logarithm of zero or negative value. Argument 
floating-point-input-value is less than or equal 
to 0.0. LIB$SIGNAL copies the floating-point 
reserved operand to the mechanism argument 

vector CHF$I_MCH—SAVRO/R1. The result is the 

floating-point reserved operand unless you have 
written a condition handler to change CHF$I— 
MCH—SAVR0/R1. 
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MTH$RANDOM Random Number Generator, 



Uniformly Distributed 

The Random Number Generator, Uniformly Distributed routine is a general 
random number generator. 

FORMAT 

MTHSRANDOM seed 

RETURNS 

VMS usage: floating-point 
type: F_ floating 

access: write only 

mechanism: by value 

MTH$RANDOM returns an F-floating random number. 

ARGUMENT 

seed 

VMS usage: long word—unsigned 
type: longword (unsigned) 

access: modify 

mechanism: by reference 

The integer seed, a 32-bit number whose high-order 24 bits are converted by 
MTH$RANDOM to an F-floating random number. The seed argument is the 
address of an unsigned longword that contains this integer seed. The seed is 
modified by each call to MTH$RANDOM. 

DESCRIPTION 

This routine must be called again to obtain the next pseudorandom number. 
The seed is updated automatically. 

The result is a floating-point number that is uniformly distributed between 

0.0 inclusively and 1.0 exclusively. 

There are no restrictions on the seed, although it should be initialized 
to different values on separate runs in order to obtain different random 
sequences. MTH$RANDOM uses the following method to update the seed 
passed as the argument: 

SEED = (69069 * SEED + l)(modulo2 32 ) 
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MTH$RANDOM 


CONDITION 

VALUE 

SIGNALED 


SS$_ROPRAND 


Reserved operand. The MTH$RANDOM procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 


EXAMPLE 


RAND: PROCEDURE OPTIONS (MAIN); 

DECLARE FORSSECNDS ENTRY (FLOAT BINARY (24)) 

RETURNS (FLOAT BINARY (24)); 

DECLARE MTH$RANDOM ENTRY (FIXED BINARY (31)) 

RETURNS (FLOAT BINARY (24)); 

DECLARE TIME FLOAT BINARY (24); 

DECLARE SEED FIXED BINARY (31); 

DECLARE I FIXED BINARY (7); 

DECLARE RESULT FIXED DECIMAL (2); 

/* Get floating random time value */ 

TIME = FOR$SECNDS (0E0); 

/* Convert to fixed */ 

SEED = TIME; 

/* Generate 100 random numbers between 1 and 10 */ 

DO I = 1 TO 100; 

RESULT = 1 + FIXED ( (10E0 * MTH$RANDOM (SEED) ),31 ); 
PUT LIST (RESULT); 

END; 

END RAND; 


7 4 6 
4 4 2 
1 10 10 
6 2 3 

8 5 4 
1 2 2 
3 8 1 


5 9 

4 4 

6 7 

6 10 

9 8 

3 6 

5 


This PL/I program demonstrates the use of MTH$RANDOM. The value 
returned by FOR$SECNDS is used as the seed for the random-number 
generator to insure a different sequence each time the program is run. The 
random value returned is scaled so as to represent values between 1 and 10. 


Because this program generates random numbers, the output generated will 
be different each time the program is executed. One example of the outut 
generated by this program is as follows: 



3 2 

9 10 

5 8 

4 2 

5 7 

5 5 
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MTH$xREAL 


MTH$xREAL Real Part of a Complex Number 



The Real Part of a Complex Number routine returns the real part of a 
complex number. 

FORMAT 

MTH$REAL complex-number 

MTH$DREAL complex-number 

MTH$GREAL complex-number 

Each of the above three formats accepts as input one of the three floating¬ 
point complex types. 

RETURNS 

VMS usage: floating-point 
type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

Real part of the complex number. MTH$REAL returns an F-floating number. 
MTH$DREAL returns a D-floating number. MTH$GREAL returns a G- 
floating number. 

ARGUMENT 

complex-number 

VMS usage: complex—number 

type: F_floating complex, D_floating complex, G_floating 

complex 
access: read only 

mechanism: by reference 

The complex number whose real part is returned by MTH$REAL. The 
complex-number argument is the address of this floating-point complex 
number. For MTH$REAL, complex-number is an F-floating complex number. 
For MTH$DREAL, complex-number is a D-floating complex number. For 
MTH$GREAL, complex-number is a G-floating complex number. 

CONDITION 

VALUE 

SIGNALED 

SS$_ROPRAND Reserved operand. The MTHSxREAL procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 


of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$xREAL 


EXAMPLE 

c+ 

C This FORTRAN example forms the real 
C part of an F-floating complex number using 
C MTH$REAL and the FORTRAN random number 
C generator RAN. 

C 

C Declare Z as a complex value and MTH$REAL as a 

C REAL*4 value. MTH$REAL will return the real 

C part of Z: Z_NEW = MTH$REAL(Z). 

C- 

COMPLEX Z 
COMPLEX CMPLX 
REAL*4 MTH$REAL 
INTEGER M 
M = 1234567 

C+ 

C Generate a random complex number with the FORTRAN 
C generic CMPLX. 

C- 

Z = CMPLX(RAN(M),RAN(M)) 

C+ 

C Z is a complex number (r,i) with real part "r" and imaginary 
C part "i". 

C- 

TYPE *, * The complex number z is’,z 
TYPE *, * It has real part *,MTH$REAL(Z) 

END 

This FORTRAN example demonstrates the use of MTH$REAL. The output of 
this program is as follows: 

The complex number z is (0.8535407,0.2043402) 

It has real part 0.8535407 
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MTH$xSIN Sine of Angle Expressed in Radians 



The Sine of Angle Expressed in Radians routine returns the sine of a given 
angle (in radians). 

FORMAT 

MTH$SIN angle-in-radians 

MTH$DSIN angle-in-radians 

MTH$GSIN angle-in-radians 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$SIN_R4 

MTH$DSIN_R7 

MTH$GSIN_R7 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F_floating, D_floating, G-floating 

access: write only 

mechanism: by value 

Sine of the angle specified by angle-in-radians. MTH$SIN returns an F- 
floating number. MTH$DSIN returns a D-floating number. MTH$GSIN 
returns a G-floating number. 

ARGUMENTS 

angle-in-radians 

VMS usage: floating-point 
type: F_floating, D_floating, G-floating 

access: read only 

mechanism: by reference 

Angle (in radians). The angle-in-radians argument is the address of a 
floating-point number that is this angle. For MTH$SIN, angle-in-radians 
specifies an F-floating number. For MTH$DSIN, angle-in-radians specifies a 
D-floating number. For MTH$GSIN, angle-in-radians specifies a G-floating 
number. 

DESCRIPTION 

See the MTH$SINCOS routine for the algorithm used to compute this sine. 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HSIN. 
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MTH$xSIN 


CONDITION 

VALUE 

SIGNALED 


SS$_ROPRAND 


Reserved operand. The MTH$xSIN procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 


MTH-125 







MTH$xSINCOS 


MTH$xSINCOS Sine and Cosine of Angle 



Expressed in Radians 

The Sine and Cosine of Angle Expressed in Radians routine returns the 
sine and cosine of a given angle (in radians). 

FORMAT 

MTH$SINCOS angle-in-radians ,sine ,cosine 
MTH$DSINCOS angle-in-radians ,sine,cosine 
MTH$GSINCOS angle-in-radians ,sine,cosine 
MTH$HSINCOS angle-in-radians ,sine,cosine 

Each of the above four formats accepts as input one of the four floating-point 
types. 

jsb entries 

MTH$SINCOS_R5 

MTH$DSINCOS_R7 

MTH$GSINCOS_R7 

MTH$HSINCOS_R7 

Each of the above four JSB entries accepts as input one of the four floating¬ 
point types. 

RETURNS 

MTH$SINCOS, MTH$DSINCOS, MTH$GSINCOS, and MTH$HSINCOS 
return the sine and cosine of the input angle by reference in the sine and 
cosine arguments. 

ARGUMENTS 

angle-in-radians 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating, H-floating 

access: read only 

mechanism: by reference 

Angle (in radians) whose sine and cosine are to be returned. The angle- 
in-radians argument is the address of a floating-point number that 
is this angle. For MTH$SINCOS, angle-in-radians is an F-floating 
number. For MTHSDSINCOS, angle-in-radians is a D-floating number. 

For MTH$GSINCOS, angle-in-radians is a G-floating number. For 
MTH$HSINCOS, angle-in-radians is an H-floating number. 

sine 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating, H_floating 

access: write only 

mechanism: by reference 

Sine of the angle specified by angle-in-radians. The sine argument is the 
address of a floating-point number. MTH$SINCOS writes an F-floating 
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MTH$xSINCOS 


DESCRIPTION 


number into sine. MTH$DSINCOS writes a D-floating number into sine. 
MTH$GSINCOS writes a G-floating number into sine. MTH$HSINCOS 
writes an H-floating number into sine. 

cosine 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating, H-floating 

access: write only 

mechanism: by reference 

Cosine of the angle specified by angle-in-radians. The cosine argument is 
the address of a floating-point number. MTH$SINCOS writes an F-floating 
number into cosine. MTH$DSINCOS writes a D-floating number into cosine, 
MTH$GSINCOS writes a G-floating number into cosine. MTH$HSINCOS 
writes an H-floating number into cosine. 


All routines with JSB entry points accept a single argument in R0:Rm, where 

m, which is defined below, is dependent on the data type. 


Data Type 

m 

F_floating 

0 

D_floating 

1 

G_floating 

1 

H-floating 

3 


In general, Run-Time Library routines with JSB entry points return one value 
in R0:Rm. The MTH$SINCOS routine returns two values, however. The sine 
of angle-in-radians is returned in R0:Rm and the cosine of angle-in-radians 
is returned in (R <m+l> :R <2*m+l> ). 

In radians, the computation of zSIN(X) and zCOS(X) is based on the following 
polynomial expansions: 

sin(X) = X - X 3 /(3!) + X 5 /(5!) - X 7 /(7!)... 

= X + X * P{X 2 ), where 
P{y) = 2//(3!) + J/ 2 /(5!) + 2/ 3 /(7!)— 

cos(X) = 1 - X 2 /(2!) + * 4 /(4!) - X 6 /(6!)... 

= Q(X 2 ), where 

Q(y) = (i - v/(2!) + y 2 /( 4!) + y 3 /(6!)-) 

1 If \X\ < 

then zSIN(X) = X and zCOS{X) = 1 
(see the section on MTH$zCOSH for 
the definition of f) 

2 If 2-7/ 2 <|X| < tt/4, 

then zSIN(X) = X + P{X 2 ) 
and zCOS(X) = Q{X 2 ) 

3 If 7 t/4<|X| and X > 0, 

a. Let J = INT(X/(tt/4)) 
and I = Jmodulo 8 
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b. If J is even, let Y = X - J * (n/4) 
otherwise, 

let Y = (J + 1) * (tt/4) - X 

With the above definitions, the following table relates zSIN(X) 
and zCOS(X) to zSIN(Y) and zCOS(Y): 


Value of / 

zSIN(X) 

zCOS(X) 

0 

zSIN(Y) 

zCOS(Y) 

1 

zCOS(Y) 

zSIN(Y) 

2 

zCOS(Y) 

-zSIN(Y) 

3 

zSIN(Y) 

-zCOS(Y) 

4 

-zSIN(Y) 

-zCOS(Y) 

5 

-zCOS(Y) 

-zSIN(Y) 

6 

-zCOS(Y) 

zSIN(Y) 

7 

-zSIN(Y) 

zCOS(Y) 


c. zSIN(Y) and zCOS(Y) are computed as follows: 
zSIN(Y) = Y + P(Y 2 ), 
and zCOS{Y) = Q(Y 2 ) 

4 If */4<\X\ and X < 0, 

then zSIN(X) = -zS/AT(|X|) 
and zCOS{X ) = zCOS{\X\) 


CONDITION 

VALUE 

RETURNED 


SS$_ROPRAND 


Reserved operand. The MTH$xSINCOS procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$xSINCOSD Sine and Cosine of Angle 

Expressed in Degrees 

The Sine and Cosine of Angle Expressed in Degrees routine returns the 
sine and cosine of a given angle (in degrees). 

FORMAT 

MTH$SINCOSD angle-in-degrees ,sine,cosine 
MTH$DSINCOSD angle-in-degrees ,sine,cosine 
MTH$GSINCOSD angle-in-degrees ,sine,cosine 
MTH$HSINCOSD angle-in-degrees ,sine,cosine 

Each of the above four formats accepts as input one of the four floating-point 
types. 

jsb entries 

MTH$SINCOSD_R5 

MTH$DSINCOSD_R7 

MTH$GSINCOSD_R7 

MTH$HSINCOSD_R7 

Each of the above four JSB entries accepts as input one of the four floating¬ 
point types. 

RETURNS 

MTH$SINCOSD, MTH$DSINCOSD, MTH$GSINCOSD, and 
MTH$HSINCOSD return the sine and cosine of the input angle by reference 
in the sine and cosine arguments. 


ARGUMENTS angle-in-degrees 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating, H_floating 

access: read only 

mechanism: by reference 


Angle (in degrees) whose sine and cosine are returned by MTH$xSINCOSD. 
The angle-in-degrees argument is the address of a floating-point number 
that is this angle. For MTH$SINCOSD, angle-in-degrees is an F-floating 
number. For MTH$DSINCOSD, angle-in-degrees is a D-floating number. 
For MTH$GSINCOSD, angle-in-degrees is a G-floating number. For 
MTH$HSINCOSD, angle-in-degrees is an H-floating number. 


sine 

VMS usage: 
type: 
access: 
mechanism: 


floating-point 

F_floating, D_floating, G_floating, H_floating 
write only 
by reference 


Sine of the angle specified by angle-in-degrees. The sine argument is the 
address of a floating-point number. MTH$SINCOSD writes an F-floating 
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DESCRIPTION 


number into sine. MTH$DSINCOSD writes a D-floating number into sine. 
MTH$GSINCOSD writes a G-floating number into sine. MTH$HSINCOSD 
writes an H-floating number into sine. 

cosine 

VMS usage: floating-point 

type: F_floating, D_floating, G-floating, H_floating 

access: write only 

mechanism: by reference 

Cosine of the angle specified by angle-in-degrees. The cosine argument 
is the address of a floating-point number. MTH$SINCOSD writes an F- 
floating number into cosine. MTH$DSINCOSD writes a D-floating number 
into cosine. MTH$GSINCOSD writes a G-floating number into cosine. 
MTH$HSINCOSD writes an H-floating number into cosine. 


All routines with JSB entry points accept a single argument in R0:Rm, where 
m, which is defined below, is dependent on the data type. 


Data Type 

m 

F__floating 

0 

D__floating 

1 

G—floating 

1 

H floating 

3 


In general, Run-Time Library routines with JSB entry points return one value 
in R0:Rm. The MTH$SINCOSD routine returns two values, however. The 
sine of angle-in-degrees is returned in R0:Rm and the cosine of angle-in¬ 
degrees is returned in (R <m+l> :R <2*m+l> ). 

In degrees, the computation of zSIND(X) and zCOSD(X) is based on the 
following polynomial expansions: 

SIND(X) = (C*X)-(C* X) 3 /(3!)+ 

(C*X) 5 /(5!)-(C*X) 7 /(7!)... 

= X/2 6 + X * P(X 2 ), where 
P(y) = -0/(3!) + y 2 /( 5!) - y 3 /{ 7!)... 

COSD{X) = 1 - (C * X) 2 /{2\)+ 

(C*X) 4 /(4!) - (C*X) 6 /( 6!)... 

= Q(X 2 ), where 

Q(V) = 1 - 0/(2!) + y 2 /(4!) - y 3 /(6!)... 
and C = 7 t/180 

1 If |A'| < (180/7 t) *2~ 2 and underflow signaling is enabled, 
underflow is signaled for zSIND(X) and zSINCOSD(X). 

See MTH$zCOSH for the definition of e. 

OTHERWISE: 

2 If \X\ < (180/tt) * 2<-// 2 ), 

then zSIND(X) = (tt/180) * X and zCOSD{X) = 1. 

(See MTH$zCOSH for the definition of /.) 
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3 If (180 /tt) * 2<-// 2 > <1X1 <45 

then zSIND{X) = X/2 6 + P{X 2 ) 
and zCOSD(X) = Q{X 2 ) 

4 If 45<|X| and X > 0, 

a. Let J = I NT (X/{4b)) and 

I = J modulo 8 

b. If J is even, let Y = X - J * 45; 
otherwise, let Y = ( J + 1) * 45 — X. 

With the above definitions, the following table relates 
zSIND(X) and zCOSD(X) to zSIND(Y) and zCOSD(Y): 


Value of / 

zSIND(X) 

zCOSD(X) 

0 

zSIND(Y) 

zCOSD(Y) 

1 

zCOSD(Y) 

zSIND(Y) 

2 

zCOSD(Y) 

-zSIND(Y) 

3 

zSIND(Y) 

-zCOSD(Y) 

4 

-zSIND(Y) 

-zCOSD(Y) 

5 

-zCOSD(Y) 

-zSIND(Y) 

6 

-zCOSD(Y) 

zSIND(Y) 

7 

-zSIND(Y) 

zCOSD(Y) 


c. zSIND(Y) and zCOSD(Y) are computed as follows: 

zSIND(Y) = Y/2 6 + P{Y 2 ) 
zCOSD{Y) = Q{Y 2 ) 

d. If 45<|X| and X < 0, 

then zSIND(X) = -zSIND{\X\) 
and zCOSD(X) = zCOSD(\X \) 


SS$_ROPRAND Reserved operand. The MTHSxSINCOSD 

procedure encountered a floating-point reserved 
operand due to incorrect user input. A floating¬ 
point reserved operand is a floating-point datum 
with a sign bit of 1 and a biased exponent of zero 
Floating-point reserved operands are reserved for 
future use by DIGITAL. 

MTH$_FLOUNDMAT Floating-point underflow in math library. The 

absolute value of the input angle is less than 
180/tt * 2" m (where m = 128 for F-floating and 
D-floating, 1,024 for G-floating, and 16,384 for 
H-floating). 


CONDITION 

VALUES 

SIGNALED 
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MTH$xSIND Sine of Angle Expressed in Degrees 



The Sine of Angle Expressed in Degrees routine returns the sine of a given 
angle (in degrees). 

FORMAT 

MTH$SIND angle-in-degrees 

MTH$DSIND angle-in-degrees 

MTH$GSIND angle-in-degrees 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$SIND_R4 

MTH$DSIND_R7 

MTH$GSIND_R7 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

The sine of the angle. MTH$SIND returns an F-floating number. 
MTH$DSIND returns a D-floating number. MTH$GSIND returns a G-floating 
number. 

ARGUMENTS 

angle-in-degrees 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

Angle (in degrees). The angle-in-degrees argument is the address of a 
floating-point number that is this angle. For MTH$SIND, angle-in-degrees 
specifies an F-floating number. For MTH$DSIND, angle-in-degrees specifies 
a D-floating number. For MTH$GSIND, angle-in-degrees specifies a G- 
floating number. 

DESCRIPTION 

See MTH$SINCOSD for the algorithm that is used to compute the sine. 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HSIND. 
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CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_FLOUNDMAT 


Reserved operand. The MTHSSIND procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased ecponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point underflow in math library. The 
absolute value of the input angle is less than 
180/tt * 2" m (where m = 128 for F-floating and 
D-floating, and 1,024 for G-floating). 
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MTH$xSINH Hyperbolic Sine 

The Hyperbolic Sine routine returns the hyperbolic sine of the input value 
specified by floating-point-input-value. 

FORMAT 

MTH$SI N H floating-point-input-value 

MTH$DSI N H floating-point-input-value 

MTH$GSI N H floating-point-input-value 

Each of the above formats accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

The hyperbolic sine of floating-point-input-value. MTH$SINH returns an 
F-floating number. MTH$DSINH returns a D-floating number. MTH$GSINH 
returns a G-floating number. 

ARGUMENTS 

floating-point-input- value 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the address 
of a floating-point number that is this value. For MTH$SINH, floating¬ 
point-input-value specifies an F-floating number. For MTH$DSINH, 
floating-point-input-value specifies a D-floating number. For MTH$GSINH, 
floating-point-input-value specifies a G-floating number. 


DESCRIPTION 


Computation of the hyperbolic sine function depends on the magnitude of 
the input argument. The range of the function is partitioned using four data 
type dependent constants: a(z), b(z), and c(z). The subscript z indicates the 
data type. The constants depend on the number of exponent bits (e) and the 
number of fraction bits (/) associated with the data type (z). 

The values of e and / are: 


z 

e 

f 

F 

8 

24 

D 

8 

56 

G 

11 

53 
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The values of the constants in terms of e and f are: 


Variable 

Value 

a(z) 

2<-//2> 

b(z) 

CEILING[ (/ + l)/2 * ln(2)] 

c(z) 

(2 <e-1> * ln(2)) 

Based on the above definitions, zSINH(X) is computed as follows: 

Value of X 

Value Returned 

IXI < a(z) 

X 

a(z) < IXI < 1.0 

zSINH(X) is computed using a 
power series expansion in |X| 2 

1.0 < IXI < b(z) 

(zEXP(X) - zEXP(-X ))/2 

b(z) <1X1 < c(z) 

SIGN(X) * zEXP( |X|)/2 

c(z) < IXI 

Overflow occurs 


The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HSINH. 


CONDITION 

VALUES 

SIGNALED 


SS$_ROPRAND 


MTH$_FLOOVEMAT 


Reserved operand. The MTH$xSINH procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

Floating-point overflow in Math Library: the 
absolute value of floating-point-input-value 
is greater than yyy. LIB$SIGNAL copies the 
floating-point reserved operand to the mechanism 

argument vector CHF$I_MCH_SAVR0/R1. The 

result is the floating-point reserved operand unless 
you have written a condition handler to change 
CHF$L _MCH _S A VRO/R1. 

The values of yyy are approximately: 

MTHSSINH—88.722 

MTH$DSINH—88.722 

MTH$GSINH—709.782 
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MTH$xSQRT Square Root 



The Square Root routine returns the square root of the input value 

floating-point-input-value. 

FORMAT 

MTH$SQRT floating-point-input-value 

MTH$DSQRT floating-point-input-value 

MTH$GSQRT floating-point-input-value 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$SQRT_R3 

MTH$DSQRT_R5 

MTH$GSQRT_R5 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

The square root of floating-point-input-value. MTH$SQRT returns an F- 
floating number. MTH$DSQRT returns a D-floating number. MTH$GSQRT 
returns a G-floating number. 

ARGUMENTS 

floating-point-input- value 

VMS usage: floating-point 
type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

Input value. The floating-point-input-value argument is the address of 
a floating-point number that contains this input value. For MTH$SQRT, 
floating-point-input-value specifies an F-floating number. For 

MTH$DSQRT, floating-point-input-value specifies a D-floating number. 

For MTH$GSQRT, floating-point-input-value specifies a G-floating number. 

DESCRIPTION 

The square root of X is computed as follows: 

If X < 0, an error is signaled. 

Let X = 2 K *F 

where: 

K is the exponential part of the floating-point data 

F is the fractional part of the floating-point data 
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CONDITION 

VALUES 

SIGNALED 


If K is even: 

X = 2 (2 * p| * F, 

zSQRT(X) = 2 P * zSQRT(F), 

1/2 <F < 1, where P = K/2 

If K is odd: 

X = 2 (2 * p+1) * F = 2 (2 * p+2) * (F/2), 
zSQRT(X) = 2 (P+1) * zSQRT(F/ 2), 
l/4<F/2 < 1/2, where p = (K-1)/2 

Let F' = A * F + B, when K is even: 

A = 0.95F6198 (hex) 

B = 0.6BA5918 (hex) 

Let F' = A* (F/2) + B, when K is odd: 

A = 0.D413CCC (hex) 

B = 0.4C1E248 (hex) 

Let K' = P, when K is even 
Let K' - P+1, when K is odd 

Let Y[0] = 2 k> * F' be a straight line approximation within the given interval 
using coefficients A and B which minimize the absolute error at the midpoint 
and endpoint. 

Starting with Y[0], n Newton-Raphson iterations are performed: 

Y[rH-l] = l/2*(Y[n] + X/K[n]) 

where n = 2, 3, or 3 for z = F-floating, D-floating, or G-floating, respectively. 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HSQRT. 


SS$_ROPRAND Reserved operand. The MTHSxSQRT procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_SQUROONEG Square root of negative number. Argument 

floating-point-input-value is less than 0.0. 
LIBSSIGNAL copies the floating-point reserved 
operand to the mechanism argument vector 

CHF$I_MCH_SAVR0/R1. The result is the 

floating-point reserved operand unless you have 
written a condition handler to change CHF$L_ 
MCH—SAVR0/R1. 
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MTH$xTAN Tangent of Angle Expressed in 



Radians 

The Tangent of Angle Expressed in Radians routine returns the tangent of 
a given angle (in radians). 

FORMAT 

MTH$TAN angle-in-radians 

MTH$DTAN angle-in-radians 

MTH$GTAN angle-in-radians 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$TAN_R4 

MTH$DTAN_R7 

MTH$GTAN_R7 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

The tangent of the angle specified by angle-in-radians. MTH$TAN 
returns an F-floating number. MTH$DTAN returns a D-floating number. 
MTH$GTAN returns a G-floating number. 

ARGUMENTS 

angle-in-radians 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

The input angle (in radians). The angle-in-radians argument is the address of 
a floating-point number that is this angle. For MTH$TAN, angle-in-radians 
specifies an F-floating number. For MTH$DTAN, angle-in-radians specifies a 
D-floating number. For MTH$GTAN, angle-in-radians specifies a G-floating 
number. 
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DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


When the input argument is expressed in radians, the tangent function is 
computed as follows: 

1 If |X| < then zTAN(X) = X (see the section on MTH$zCOSH 

for the definition of /) 

2 Otherwise, call MTH$zSINCOS to obtain zSIN(X) and zCOS(X); then 

a. If zCOS(X) = 0, signal overflow 

b. Otherwise, zTAN{X) = zSIN{X)/zCOS{X) 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HTAN. 


SSS—ROPRAND Reserved operand. The MTH$xTAN procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_FLOOVEMAT Floating-point overflow in Math Library. 
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MTH$xTAND Tangent of Angle Expressed in 



Degrees 

The Tangent of Angle Expressed in Degrees routine returns the tangent of 
a given angle (in degrees). 

FORMAT 

MTH$TAND angle-in-degrees 

MTH$DTAND angle-in-degrees 

MTH$GTAND angle-in-degrees 

Each of the above formats accepts as input one of the floating-point types. 

jsb entries 

MTH$TAND_R4 

MTH$DTAND_R7 

MTH$GTAND_R7 

Each of the above JSB entries accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 
type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

Tangent of the angle specified by angle-in-degrees. MTH$TAND returns 
an F-floating number. MTH$DTAND returns a D-floating number. 
MTH$GTAND returns a G-floating number. 

ARGUMENTS 

angle-in-degrees 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

The input angle (in degrees). The angle-in-degrees argument is the address 
of a floating-point number which is this angle. For MTH$TAND, angle-in- 
degrees specifies an F-floating number. For MTH$DTAND, angle-in-degrees 
specifies a D-floating number. For MTH$GTAND, angle-in-degrees specifies 
a G-floating number. 
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DESCRIPTION 


CONDITION 

VALUES 

SIGNALED 


When the input argument is expressed in degrees, the tangent function is 
computed as follows: 

1 If |X| < (180/ir)*2^ -2 ^ e_1 ^ and underflow signaling is enabled, underflow 
is signaled (see the section on MTH$zCOSH for the definition of e). 

2 Otherwise, if \X\ < (180/tt) * 2then zTAND(X) = (tt/180) * X. See 
the description of MTH$zCOSH for the definition of f. 

3 Otherwise, call MTH$zSINCOSD to obtain zSIND(X) and zCOSD(X). 

a. Then, if zCOSD(X) = 0, signal overflow 

b. Else, zTAND(X) = zSIND{X)/zCOSD{X ) 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HTAND. 


SS$_ROPRAND Reserved operand. The MTH$xTAND procedure 

encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 

MTH$_FLOOVEMAT Floating-point overflow in Math Library. 
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MTH$xTANH Compute the Hyperbolic Tangent 



The Compute the Hyperbolic Tangent routine returns the hyperbolic 
tangent of the input value. 

FORMAT 

MTH$T AN H floating-point-input-value 

MTH$DTAN H floating-point-input-value 

MTH$GTAN H floating-point-input-value 

Each of the above formats accepts as input one of the floating-point types. 

RETURNS 

VMS usage: floating-point 

type: F_floating, D_floating, G_floating 

access: write only 

mechanism: by value 

The hyperbolic tangent of floating-point-input-value. MTH$TANH 
returns an F-floating number. MTH$DTANH returns a D-floating number. 
MTH$GTANH returns a G-floating number. Unlike the other three routines, 
MTH$HTANH returns the hyperbolic tangent by reference in the h-tanh 
argument. 

ARGUMENTS 

floating-point-input- value 

VMS usage: floating-point 
type: F_floating, D_floating, G_floating 

access: read only 

mechanism: by reference 

The input value. The floating-point-input-value argument is the 
address of a floating-point number that contains this input value. For 
MTH$TANH, floating-point-input-value specifies an F-floating number. 

For MTH$DTANH, floating-point-input-value specifies a D-floating number. 
For MTH$GTANH, floating-point-input-value specifies a G-floating number. 
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DESCRIPTION 


CONDITION 

VALUE 

SIGNALED 


In calculating the hyperbolic tangent of x, the values of g and h are: 


z 

9 

h 

F 

12 

10 

D 

28 

21 

G 

26 

20 

For MTH$TANH, MTH$DTANH, and MTH$GTANH the hyperbolic tangent 

of x is then computed as follows: 

Value of x 

Hyperbolic Tangent Returned 


|x|<2 -9 

X 

2~° < |X| < 0.5 

xTANH(X) = X + X 3 * R{X 2 ), where R(X 2 ) is a 
rational function of X 2 . 

0.5<|X| < 1.0 

xTANH(X) = xTANH(xHI)+xTANH(xLO)*C/B 
where C = 1 - xTANH(xHI) * xTANH(xHI), 
B=l + xTANH(xHI) * xTANH(xLO), 
xHI = 1/2 + A/16 + 1/32 for N=0,1.7, 

and xLO = X — xHI. 

1.0 <\X\<h 

xTANH(X) = (xEXP(2*X)-l)/{xEXP(2*X) + l) 

h<\X\ 

xTANH(X) = sign(X) * 1 

The routine description for the H-floating point version of this routine is listed 
alphabetically under MTH$HTANH. 

SS$_ROPRAND 

Reserved operand. The MTHSxTANH procedure 
encountered a floating-point reserved operand due 
to incorrect user input. A floating-point reserved 
operand is a floating-point datum with a sign bit 
of 1 and a biased exponent of zero. Floating-point 
reserved operands are reserved for future use by 
DIGITAL. 
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MTH$UMAX 


MTH$UMAX Compute Unsigned Maximum 



The Compute Unsigned Maximum routine computes the unsigned 
longword maximum of n unsigned longword arguments, where n is greater 
than or equal to 1. 

FORMAT 

MTH$UMAX argument [argument,...] 

RETURNS 

VMS usage: longword_unsigned 
type: longword (unsigned) 

access: write only 

mechanism: by value 

Maximum value returned by MTH$UMAX. 

ARGUMENTS 

argument 

VMS usage: longword_unsigned 
type: longword (unsigned) 

access: read only 

mechanism: by reference 

Argument whose maximum MTH$UMAX computes. Each argument 
argument is an unsigned longword that contains one of these values. 

argument 

VMS usage: longword_unsigned 
type: longword (unsigned) 

access: read only 

mechanism: by reference 

Additional arguments whose maximum MTH$UMAX computes. Each 
argument argument is an unsigned longword that contains one of these 
values. 

DESCRIPTION 

MTH$UMAX is the unsigned version of MTH$JMAX0. 

CONDITION 

VALUES 

RETURNED 

None. 
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MTHSUMIN 


MTH$UMIN Compute Unsigned Minimum 



The Compute Unsigned Minimum routine computes the unsigned longword 
minimum of n unsigned longword arguments, where n is greater than or 
equal to 1. 

FORMAT 

MTH$UMIN argument [argument,...] 

RETURNS 

VMS usage: longword—unsigned 
type: longword (unsigned) 

access: write only 

mechanism: by value 

Minimum value returned by MTH$UMIN. 

ARGUMENTS 

argument 

VMS usage: longword—unsigned 
type: longword (unsigned) 

access: read only 

mechanism: by reference 

Argument whose minimum MTH$UMIN computes. Each argument 
argument is an unsigned longword that contains one of these values. 

argument 

VMS usage: longword—unsigned 
type: longword (unsigned) 

access: read only 

mechanism: by reference 

Additional arguments whose minimum MTH$UMIN computes. Each 
argument argument is an unsigned longword that contains one of these 
values. 

DESCRIPTION 

MTH$UMIN is the unsigned version of MTH$JMIN0. 

CONDITION 

VALUES 

RETURNED 

None. 
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y\ Undocumented MTH$ Routines 


This appendix lists all of the entry point and argument information for the 
MTH$ routines not documented in the MTH$ Reference Section of this 
manual. 


Table A—1 Undocumented MTH$ Routines 


Routine Name 


Entry Point Information 

MTH$ABS 

Format: 

Returns: 

f-floating: 

F-floating Absolute Value Routine 

MTHSABS f-floating 

floating-point, F-floating, write only, by value 
floating-point, F_floating, read only, by reference 

MTHSDABS 

Format: 

Returns: 

d-floating: 

D-floating Absolute Value Routine 

MTH$DABS d-floating 

floating-point, D_floating, write only, by value 
floating-point, D_floating, read only, by reference 

MTH$GABS 

Format: 

Returns: 

g-floating: 

G-floating Absolute Value Routine 

MTHSGABS g-floating 

floating-point, G_floating, write only, by value 
floating-point, G_floating, read only, by reference 

MTHSHABS 

Format: 

Returns: 

h-abs-val: 

h-floating: 

H-floating Absolute Value Routine 

MTH$ABS h-abs-val, h-floating 

None 

floating-point, H_floating, write only, by reference 
floating-point, H_floating, read only, by reference 

MTHSIIABS 

Format: 

Returns: 

word: 

Word Absolute Value Routine 

MTHSIIABS word 

word-signed, word (signed), write only, by value 
word—signed, word (signed), read only, by reference 

MTHSJIABS 

Format: 

Returns: 

long word: 

Long word Absolute Value Routine 

MTHSJIABS longword 

longword—signed, longword (signed), write only, by value 
longword—signed, longword (signed), read only, by reference 
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Table A-1 (Cont.) Undocumented MTH$ Routines 


Routine Name 

Entry Point Information 

MTH$IIAND 

Format: 

Returns: 

wordl: 

word 2: 

Bitwise AND of Two Word Parameters Routine 

MTH$IIAND wordl, word2 

word—unsigned, word (unsigned), write only, by value 
word-unsigned, word (unsigned), read only, by reference 
word—unsigned, word (unsigned), read only, by reference 

MTH$JIAND 

Format: 

Returns: 

longwordl: 
longword2: 

Bitwise AND of Two Long word Parameters Routine 

MTH$JIAND longwordl, longword2 

longword—unsigned, longword (unsigned), write only, by value 
longword—unsigned, longword (unsigned), read only, by reference 
longword—unsigned, longword (unsigned), read only, by reference 

MTHSDBLE 

Format: 

Returns: 

f-floating: 

Convert F-floating to D-floating (Exact) Routine 

MTHSDBLE f-floating 

floating-point, D_floating, write only, by value 
floating-point, F_floating, read only, by reference 

MTH$GDBLE 

Format: 

Returns: 

f-floating: 

Convert F-floating to G-floating (Exact) Routine 

MTHSGDBLE f-floating 

floating-point, G_floating, write only, by value 

floating-point, F_floating, read only, by reference 

MTH$DIM 

Format: 

Returns: 

f-floating 1: 
f-floating2: 

Positive Difference of Two F-floating Parameters Routine 

MTHSDIM f-floating 1, f-floating2 
floating-point, F_floating, write only, by value 
floating-point, F_floating, read only, by reference 
floating-point, F_floating, read only, by reference 

MTHSDDIM 

Format: 

Returns: 

d-floatingl: 

d-floating2: 

Positive Difference of Two D-floating Parameters Routine 

MTHSDDIM d-floatingl, d-floating2 
floating-point, D_floating, write only, by value 
floating-point, D_floating, read only, by reference 
floating-point, D_floating, read only, by reference 
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Table A-1 (Cont.) Undocumented MTH$ Routines 


Routine Name 

Entry Point Information 

MTH$GDIM 

Format: 

Returns: 

g-floatingl: 
g-floating2: 

Positive Difference of Two G-floafing Parameters Routine 

MTHSGDIM g-floatingl, g-floating2 
floating-point, G_floating, write only, by value 
floating-point, G_floating, read only, by reference 
floating-point, G_floating, read only, by reference 

MTH$HDIM 

Format: 

Returns: 

h-floating: 
h-floatingl: 
h-floating2: 

Positive Difference of Two H-floating Parameters Routine 

MTH$HDIM h-floating, h-floatingl, h-floating2 

None 

floating-point, H_floating, write only, by reference 
floating-point, H_floating, read only, by reference 
floating-point, H_floating, read only, by reference 

MTH$IIDIM 

Format: 

Returns: 

wordl: 

word 2: 

Positive Difference of Two Word Parameters Routine 

MTH$IIDIM wordl, word2 
word—signed, word (signed), write only, by value 
word—signed, word (signed), read only, by reference 
word—signed, word (signed), read only, by reference 

MTHSJIDIM 

Format: 

Returns: 

long wordl: 
longword2: 

Positive Difference of Two Long word Parameters Routine 
MTHSJIDIM longwordl, longword2 
longword—signed, longword (signed), write only, by value 
longword—signed, longword (signed), read only, by reference 
longword—signed, longword (signed), read only, by reference 

MTH$IIEOR 

Format: 

Returns: 

wordl: 

word 2: 

Bitwise Exclusive OR of Two Word Parameters Routine 

MTHSIIEOR wordl, word2 

word—unsigned, word (unsigned), write only, by value 
word-unsigned, word (unsigned), read only, by reference 
word—unsigned, word (unsigned), read only, by reference 

MTH$JIEOR 

Format: 

Returns: 

longwordl: 
longword2: 

Bitwise Exclusive OR of Two Longword Parameters Routine 
MTHSJIEOR longwordl, longword2 

longword—unsigned, longword (unsigned), write only, by value 
longword—unsigned, longword (unsigned), read only, by reference 
longword—unsigned, longword (unsigned), read only, by reference 
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Routine Name 

Entry Point Information 

MTHSIIFIX 

Format: 

Returns: 

f-floating: 

Convert F-floating to Word (Truncated) Routine 

MTHSIIFIX f-floating 

word—signed, word (signed), write only, by value 
floating-point, F_floating, read only, by reference 

MTH$JIFIX 

Format: 

Returns: 

f-floating: 

Convert F-floating to Longword (Truncated) Routine 

MTHSJIFIX f-floating 

longword_signed, longword (signed), write only, by value 
floating-point, F_floating, read only, by reference 

MTHSFLOATI 

Format: 

Returns: 

word: 

Convert Word to F-floating (Exact) Routine 

MTHSFLOATI word 

floating-point, F_floating, write only, by value 
word—signed, word (signed), read only, by reference 

MTH$DFLOTI 

Format: 

Returns: 

word: 

Convert Word to D-floating (Exact) Routine 

MTHSDFLOTI word 

floating-point, D_floating, write only, by value 
word-signed, word (signed), read only, by reference 

MTH$GFLOTI 

Format: 

Returns: 

word: 

Convert Word to G-floating (Exact) Routine 

MTHSGFLOTI word 

floating-point, G_floating, write only, by value 
word-signed, word (signed), read only, by reference 

MTH$FLOATJ 

Format: 

Returns: 

long word: 

Convert Longword to F-floating (Exact) Routine 

MTHSFLOATJ longword 

floating-point, F-floating, write only, by value 

longword_signed, longword (signed), read only, by reference 

MTH$DFLOTJ 

Format: 

Returns: 

long word: 

Convert Longword to D-floating (Exact) Routine 

MTHSDFLOTJ longword 

floating-point, D_floating, write only, by value 

longword—signed, longword (signed), read only, by reference 
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Routine Name 


Entry Point Information 

MTH$GFLOTJ 

Format: 

Returns: 

long word: 

Convert Longword to G-floating (Exact) Routine 
MTHSGFLOTJ longword 
floating-point, G-floating, write only, by value 
longword—signed, longword (signed), read only, by reference 

MTHSFLOOR 

Format: 

JSB: 

Returns: 

f-floating: 

Convert F-floating to Greatest F-floating Integer Routine 
MTHSFLOOR f-floating 

MTHSFLOOR—R1 f-floating 

floating-point, F_floating, write only, by value 

floating-point, F_floating, read only, by reference 

MTHSDFLOOR 

Format: 

JSB: 

Returns: 

d-floating: 

Convert D-floating to Greatest D-floating Integer Routine 
MTHSDFLOOR d-floating 

MTHSDFLOOR—R3 d-floating 
floating-point, D_floating, write only, by value 
floating-point, D_floating, read only, by reference 

MTH$GFLOOR 

Format: 

JSB: 

Returns: 

g-floating: 

Convert G-floating to Greatest G-floating Integer Routine 
MTHSGFLOOR g-floating 

MTHSGFLOOR-R3 g-floating 
floating-point, G_floating, write only, by value 
floating-point, G_floating, read only, by reference 

MTH$HFLOOR 

Format: 

JSB: 

Returns: 

max-h-float: 

h-floating: 

Convert H-floating to Greatest H-floating Integer Routine 
MTHSHFLOOR max-h-float, h-floating 

MTHSHFLOOR-R7 h-floating 

None 

floating-point, H_floating, write only, by reference 
floating-point, H_floating, read only, by reference 

MTH$AINT 

Format: 

JSB: 

Returns: 

f-floating: 

Convert F-floating to Truncated F-floating Routine 
MTHSAINT f-floating 

MTHSAINT-R2 f-floating 

floating-point, F_floating, write only, by value 

floating-point, F_floating, read only, by reference 
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Routine Name 


Entry Point Information 

MTHSDINT 

Format: 

JSB: 

Returns: 

d-floating: 

Convert D-floating to Truncated D-floating Routine 
MTHSDINT d-floating 

MTH$DINT_R4 d-floating 

floating-point, D_floating, write only, by value 

floating-point, D_floating, read only, by reference 

MTH$IIDINT 

Format: 

Returns: 

d-floating: 

Convert D-floating to Word (Truncated) Routine 

MTHSIIDINT d-floating 

word—signed, word (signed), write only, by value 
floating-point, D_floating, read only, by reference 

MTH$JIDINT 

Format: 

Returns: 

d-floating: 

Convert D-floating to Longword (Truncated) Routine 

MTHSJIDINT d-floating 

longword—signed, longword (signed), write only, by value 
floating-point, D_floating, read only, by reference 

MTHSGINT 

Format: 

JSB: 

Returns: 

g-floating: 

Convert G-floating to G-floating (Truncated) Routine 
MTHSGINT g-floating 

MTH$GINT_R4 g-floating 

floating-point, G_floating, write only, by value 

floating-point, G_floating, read only, by reference 

MTH$IIGINT 

Format: 

Returns: 

g-floating: 

Convert G-floating to Word (Truncated) Routine 

MTHSIIGINT g-floating 

word-signed, word (signed), write only, by value 
floating-point, G_floating, read only, by reference 

MTH$JIGINT 

Format: 

Returns: 

g-floating: 

Convert G-floating to Longword (Truncated) Routine 

MTHSJIGINT g-floating 

longword—signed, longword (signed), write only, by value 
floating-point, G_floating, read only, by reference 

MTH$HINT 

Format: 

JSB: 

Returns: 

trunc-h-flt: 

h-floating: 

Convert H-floating to H-floating (Truncated) Routine 
MTHSHINT trunc-h-flt, h-floating 

MTH$HINT_R8 h-floating 

None 

floating-point, H_floating, write only, by reference 
floating-point, H_floating, read only, by reference 
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Routine Name 

Entry Point Information 

MTHSIIHINT 

Format: 

Returns: 

h-floating: 

Convert H-floating to Truncated Word Routine 

MTHSIIHINT h-floating 

word—signed, word (signed), write only, by value 
floating-point, H_floating, read only, by reference 

MTH$JIHINT 

Format: 

Returns: 

h-floating: 

Convert H-floating to Truncated Longword Routine 

MTHSJIHINT h-floating 

longword—signed, longword (signed), write only, by value 
floating-point, H_floating, read only, by reference 

MTH$IINT 

Format: 

Returns: 

f-floating: 

Convert F-floating to Word (Truncated) Routine 

MTHSIINT f-floating 

word—signed, word (signed), write only, by value 
floating-point, F_floating, read only, by reference 

MTH$JINT 

Format: 

Returns: 

f-floating: 

Convert F-floating to Longword (Truncated) Routine 

MTHSJINT f-floating 

longword—signed, longword (signed), write only, by value 
floating-point, F_floating, read only, by reference 

MTH$IIOR 

Format: 

Returns: 

wordl: 

word 2: 

Bitwise Inclusive OR of Two Word Parameters Routine 

MTHSIIOR wordl, word2 

word—unsigned, word (unsigned), write only, by value 
word—unsigned, word (unsigned), read only, by reference 
word—unsigned, word (unsigned), read only, by reference 

MTHSJIOR 

Format: 

Bitwise Inclusive OR of Two Longword Parameters Routine 
MTHSJIOR longword 1, longword2 

Returns: 

long wordl: 
longword2: 

longword—unsigned, longword (unsigned), write only, by value 
longword—unsigned, longword (unsigned), read only, by reference 
longword—unsigned, longword (unsigned), read only, by reference 

MTH$AIMAXO 

Format: 

Returns: 

word: 

F-floating Maximum of N Word Parameters Routine 

MTHSAIMAXO word, . . . 

floating-point, F_floating, write only, by value 
word—signed, word (signed), read only, by reference 
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Routine Name 


MTH$AJMAXO 


MTHSIMAXO 


MTHSJMAXO 


MTH$AMAX1 


MTH$DMAX1 


MTH$GMAX1 


MTH$HMAX1 


Entry Point Information 


Format: 
Returns: 
long word: 


Format: 

Returns: 

word: 


Format: 
Returns: 
long word: 


Format: 

Returns: 

f-floating: 


Format: 

Returns: 

d-floating: 


Format: 

Returns: 

g-floating: 


Format: 

Returns: 

h-float-max: 

h-floating: 


F-floating Maximum of N Longword Parameters Routine 
MTH$AJMAXO longword, . . . 
floating-point, F_floating, write only, by value 
longword—signed, longword (signed), read only, by reference 

Word Maximum of N Word Parameters Routine 
MTHSIMAXO word, . . . 

word_signed, word (signed), write only, by value 
word—signed, word (signed), read only, by reference 

Longword Maximum of N Longword Parameters Routine 
MTHSJMAXO longword, . . . 

longword—signed, longword (signed), write only, by value 
longword_signed, longword (signed), read only, by reference 

F-floating Maximum of N F-floating Parameters Routine 

MTHSAMAX1 f-floating_ 

floating-point, F_floating, write only, by value 
floating-point, F_floating, read only, by reference 

D-floating Maximum of N D-floating Parameters Routine 
MTHSDMAX1 d-floating, . . . 
floating-point, D_floating, write only, by value 
floating-point, D_floating, read only, by reference 

G-floating Maximum of N G-floating Parameters Routine 

MTHSGMAX1 g-floating_ 

floating-point, G_floating, write only, by value 
floating-point, G_floating, read only, by reference 

FI-floating Maximum of N H-floating Parameters Routine 
MTHSHMAX1 h-float-max, h-floating, . . . 

None 

floating-point, H_floating, write only, by reference 
floating-point, H_floating, read only, by reference 
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Routine Name 

Entry Point Information 

MTH$IMAX1 

Format: 

Returns: 

f-floating: 

Word Maximum of N F-floating Parameters Routine 

MTHSIMAX1 f-floating, . . . 

word—signed, word (signed), write only, by value 

floating-point, F_floating, read only, by reference 

MTH$JMAX1 

Format: 

Returns: 

f-floating: 

Longword Maximum of N F-floating Parameters Routine 

MTH$JMAX1 f-floating_ 

longword_signed, longword (signed), write only, by value 
floating-point, F_floating, read only, by reference 

MTHSAIMINO 

Format: 

Returns: 

word: 

F-floating Minimum of N Word Parameters Routine 

MTHSAIMINO word, . . . 

floating-point, F_floating, write only, by value 
word-signed, word (signed), read only, by reference 

MTH$AJMINO 

Format: 

Returns: 

long word: 

F-floating Minimum of N Longword Parameters Routine 
MTHSAJMINO longword, . . . 
floating-point, F_floating, write only, by value 
longword—signed, longword (signed), read only, by reference 

MTH$IMINO 

Format: 

Returns: 

word: 

Word Minimum of N Word Parameters Routine 

MTHSIMINO word, . . . 

word—signed, word (signed), write only, by value 
word—signed, word (signed), read only, by reference 

MTH$JMINO 

Format: 

Returns: 

long word: 

Longword Minimum of N Longword Parameters Routine 

MTHSJMINO longword, . . . 

longword—signed, longword (signed), write only, by value 
longword_signed, longword (signed), read only, by reference 

MTH$AMIN1 

Format: 

Returns: 

f-floating: 

F-floating Minimum of N F-floating Parameters Routine 

MTHSAMIN1 f-floating, . . . 

floating-point, F_floating, write only, by value 

floating-point, F_floating, read only, by reference 
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Routine Name 


Entry Point Information 

MTH$DMIN1 

Format: 

Returns: 

d-floating: 

D-floating Minimum of N D-floating Parameters Routine 
MTH$DMIN1 d-floating, . . . 
floating-point, D_floating, write only, by value 
floating-point, D_floating, read only, by reference 

MTH$GMIN1 

Format: 

Returns: 

g-floating: 

G-floating Minimum of N G-floating Parameters Routine 

MTH$GMIN1 g-floating_ 

floating-point, G_floating, write only, by value 
floating-point, G_floating, read only, by reference 

MTH$HMIN1 

Format: 

Returns: 

h-float-max: 

h-floating: 

H-floating Minimum of N H-floating Parameters Routine 
MTH$HMIN1 h-float-max, h-floating, . . . 

None 

floating-point, H_floating, write only, by reference 
floating-point, H_floating, read only, by reference 

MTH$IMIN1 

Format: 

Returns: 

f-fioating: 

Word Minimum of N F-floating Parameters Routine 

MTHSIMINI f-floating_ 

word—signed, word (signed), write only, by value 
floating-point, F_floating, read only, by reference 

MTHSJMIN1 

Format: 

Returns: 

f-floating: 

Longword Minimum of N F-floating Parameters Routine 
MTH$JMIN1 f-floating, . . . 

longword—signed, longword (signed), write only, by value 
floating-point, F_floating, read only, by reference 

MTH$AM0D 

Format: 

Returns: 

dividend: 

divisor: 

Remainder of Two F-floating Parameters Routine 
MTH$AMOD dividend, divisor 
floating-point, F_floating, write only, by value 
floating-point, F_floating, read only, by reference 
floating-point, F_floating, read only, by reference 

MTHSDMOD 

Format: 

Returns: 

dividend: 

divisor: 

Remainder of Two D-floating Parameters Routine 
MTH$DMOD dividend, divisor 
floating-point, D_floating, write only, by value 
floating-point, D_floating, read only, by reference 
floating-point, D_floating, read only, by reference 
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Routine Name 

Entry Point Information 

MTHSGMOD 

Format: 

Returns: 

dividend: 

divisor: 

Remainder of Two G-floating Parameters Routine 

MTHSGMOD dividend, divisor 
floating-point, G_floating, write only, by value 
floating-point, G_floating, read only, by reference 
floating-point, G_floating, read only, by reference 

MTH$HMOD 

Format: 

Returns: 

h-mod: 

dividend: 

divisor: 

Remainder of Two H-fbating Parameters Routine 

MTHSHMOD h-mod, dividend, divisor 

None 

floating-point, H_floating, write only, by reference 
floating-point, H-floating, read only, by reference 
floating-point, H_floating, read only, by reference 

MTH$IMOD 

Format: 

Returns: 

dividend: 

divisor: 

Remainder of Two Word Parameters Routine 

MTHSIMOD dividend, divisor 
word—signed, word (signed), write only, by value 
word—signed, word (signed), read only, by reference 
word—signed, word (signed), read only, by reference 

MTH$JMOD 

Format: 

Returns: 

dividend: 

divisor: 

Remainder of Two Longword Parameters Routine 

MTHSJMOD dividend, divisor 

longword—signed, longword (signed), write only, by value 
longword—signed, longword (signed), read only, by reference 
longword_signed, longword (signed), read only, by reference 

MTHSANINT 

Format: 

Returns: 

f-floating: 

Convert F-floating to Nearest F-floating Integer Routine 

MTHSANINT f-floating 

floating-point, F_floating, write only, by value 

floating-point, F_floating, read only, by reference 

MTH$DNINT 

Format: 

Returns: 

d-floating: 

Convert D-floating to Nearest D-floating Integer Routine 
MTHSDNINT d-floating 
floating-point, D_floating, write only, by value 
floating-point, D_floating, read only, by reference 
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Undocumented MTH$ Routines 


Table A-1 (Cont.) Undocumented MTH$ Routines 


Routine Name 


Entry Point Information 

MTH$IIDNNT 

Format: 

Returns: 

d-floating: 

Convert D-floating to Word Integer Routine 

MTHSIIDNNT d-floating 

word—signed, word (signed), write only, by value 
floating-point, D_floating, read only, by reference 

MTH$JIDNNT 

Format: 

Returns: 

d-floating: 

Convert D-floating to Nearest Longword Integer Routine 

MTH$JIDNNT d-floating 

longword_signed, longword (signed), write only, by value 
floating-point, D_floating, read only, by reference 

MTH$GNINT 

Format: 

Returns: 

g-floating: 

Convert G-floating to Nearest G-floating Integer Routine 
MTHSGNINT g-floating 
floating-point, G_floating, write only, by value 
floating-point, G_floating, read only, by reference 

MTHSIIGNNT 

Format: 

Returns: 

g-floating: 

Convert G-floating to Nearest Word Integer Routine 

MTHSIIGNNT g-floating 

word—signed, word (signed), write only, by value 
floating-point, G_floating, read only, by reference 

MTH$JIGNNT 

Format: 

Returns: 

g-floating: 

Convert G-floating to Nearest Longword Integer Routine 

MTHSJIGNNT g-floating 

longword—signed, longword (signed), write only, by value 
floating-point, G_floating, read only, by reference 

MTH$HNINT 

Format: 

Returns: 

nearst-h-flt: 

h-floating: 

Convert H-floating to Nearest H-floating Integer Routine 
MTH$HNINT nearst-h-flt, h-floating 

None 

floating-point, H_floating, write only, by reference 
floating-point, H_floating, read only, by reference 

MTH$IIHNNT 

Format: 

Returns: 

h-floating: 

Convert H-floating to Nearest Word Integer Routine 

MTHSIIHNNT h-floating 

word—signed, word (signed), write only, by value 
floating-point, H_floating, read only, by reference 
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Table A-1 (Cont.) Undocumented MTH$ Routines 


Routine Name 


Entry Point Information 

MTH$JIHNNT 

Format: 

Returns: 

h-floating: 

Convert H-floating to Nearest Longword Integer Routine 

MTH$JIHNNT h-floating 

longword—signed, longword (signed), write only, by value 
floating-point, H_floating, read only, by reference 

MTH$ININT 

Format: 

Returns: 

f-floating: 

Convert F-floating to Nearest Word Integer Routine 

MTH$ININT f-floating 

word-signed, word (signed), write only, by value 
floating-point, F_floating, read only, by reference 

MTH$JNINT 

Format: 

Returns: 

f-floating: 

Convert F-floating to Nearest Longword Integer Routine 
MTHSJNINT f-floating 

longword—signed, longword (signed), write only, by value 
floating-point, F_floating, read only, by reference 

MTH$INOT 

Format: 

Returns: 

word: 

Bitwise Complement of Word Parameter Routine 

MTHSINOT word 

word—unsigned, word (unsigned), write only, by value 
word—unsigned, word (unsigned), read only, by reference 

MTH$JNOT 

Format: 

Returns: 

long word: 

Bitwise Complement of Longword Parameter Routine 

MTHSJNOT longword 

longword—unsigned, longword (unsigned), write only, by value 
longword—unsigned, longword (unsigned), read only, by reference 

MTHSDPROD 

Format: 

Returns: 

f-floating 1: 
f-floating2: 

D-floating Product of Two F-floating Parameters Routine 
MTHSDPROD f-floating 1, f-floating2 
floating-point, D_floating, write only, by value 
floating-point, F_floating, read only, by reference 
floating-point, F_floating, read only, by reference 

MTH$GPROD 

Format: 

Returns: 

f-floating 1: 
f-floating2: 

G-floating Product of Two F-floating Parameters Routine 
MTHSGPROD f-floating 1, f-floating2 
floating-point, G_floating, write only, by value 
floating-point, F_floating, read only, by reference 
floating-point, F_floating, read only, by reference 
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Table A-1 (Cont.) Undocumented MTH$ Routines 

Routine Name Entry Point Information 


MTH$SGN 

Format: 

Returns: 

f-floating: 

F-floating Sign Function 

MTHSSGN f-floating 

longword—signed, longword (signed), write only, by reference 
floating-point, F_floating, read only, by reference 

MTH$SGN 

Format: 

Returns: 

d-floating: 

D-floating Sign Function 

MTH$SGN d-floating 

longword—signed, longword (signed), write only, by reference 
floating-point, D_floating, read only, by reference 

MTH$IISHFT 

Format: 

Returns: 

word: 

shift-cnt: 

Bitwise Shift of Word Routine 

MTH$IISHFT word, shift-cnt 

word—unsigned, word (unsigned), write only, by value 
word—unsigned, word (unsigned), read only, by reference 
word—signed, word (signed), read only, by reference 

MTHSJISHFT 

Format: 

Returns: 

longword: 

shift-cnt: 

Bitwise Shift of Longword Routine 

MTHSJISHFT longword, shift-cnt 

longword—unsigned, longword (unsigned), write only, by value 
longword—unsigned, longword (unsigned), read only, by reference 
longword—signed, longword (signed), read only, by reference 

MTH$SIGN 

Format: 

Returns: 

f-float-x: 

f-float-y: 

F-floating Transfer of Sign of Y to Sign of X Routine 

MTHSSIGN f-float-x, f-float-y 
floating-point, F_floating, write only, by value 
floating-point, F_floating, read only, by reference 
floating-point, F_floating, read only, by reference 

MTH$DSIGN 

Format: 

Returns: 

d-float-x: 

d-float-y: 

D-floating Transfer of Sign of Y to Sign of X Routine 
MTHSDSIGN d-float-x, d-float-y 
floating-point, D_floating, write only, by value 
floating-point, D_floating, read only, by reference 
floating-point, D_floating, read only, by reference 
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Table A-1 (Cont.) Undocumented MTH$ Routines 


Routine Name 


Entry Point Information 

MTH$GSIGN 

Format: 

Returns: 

g-float-x: 

g-float-y: 

G-floating Transfer of Sign of Y to Sign of X Routine 

MTHSGSIGN g-float-x, g-float-y 
floating-point, G_floating, write only, by value 
floating-point, G_floating, read only, by reference 
floating-point, G_floating, read only, by reference 

MTH$HSIGN 

Format: 

Returns: 

h-result: 

h-float-x: 

h-float-y: 

H-floating Transfer of Sign of Y to Sign of X Routine 

MTHSHSIGN h-result, h-float-x, h-float-y 

None 

floating-point, H_floating, write only, by reference 
floating-point, H_floating, read only, by reference 
floating-point, H_floating, read only, by reference 

MTH$IISIGN 

Format: 

Returns: 

word-x: 

word-y: 

Word Transfer of Sign of Y to Sign of X Routine 

MTHSHSIGN word-x, word-y 
word—signed, word (signed), write only, by value 
word—signed, word (signed), read only, by reference 
word—signed, word (signed), read only, by reference 

MTH$JISIGN 

Format: 

Returns: 

longwrd-x: 

longwrd-y: 

Longword Transfer of Sign of Y to Sign of X Routine 

MTHSJISIGN longwrd-x, longwrd-y 

longword—signed, longword (signed), write only, by reference 
longword—signed, longword (signed), read only, by reference 
longword—signed, longword (signed), read only, by reference 

MTH$SNGL 

Format: 

Returns: 

d-floating: 

Convert D-floating to F-floating (Rounded) Routine 

MTHSSNGL d-floating 

floating-point, F_floating, write only, by value 
floating-point, D_floating, read only, by reference 

MTH$SNGLG 

Format: 

Returns: 

g-floating: 

Convert G-floating to F-floating (Rounded) Routine 

MTHSSNGLG g-floating 

floating-point, F_floating, write only, by value 

floating-point, G_floating, read only, by reference 
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A D 


Absolute value 

See also Mathematics routine 
of complex number •MTH-23 
Arc cosine 

in degrees* MTH-6, MTH-71 
in radians *MTH-3, MTH-69 
Arc sine 

in degrees*MTH-11, MTH-75 
in radians* MTH-9, MTH-73 
Arc tangent 

hyperbolic •MTH-21, MTH-85 
in degrees*MTH-15, MTH-19, MTH-79, 
MTH-83 

in radians*MTH-13, MTH-17, MTH-77, 
MTH-81 

Arrays 

conversion of*MTH-64 


c 


Complex number* 1-3 ; MTH-57, MTH-59, 
MTH-112, MTH-122 
absolute value of*MTH-23 
complex exponential of*MTH-31, MTH-33 
conjugate of*MTH-44, MTH-45 
cosine of*MTH-26, MTH-28 
made from floating-point*MTH-40, MTH-42 
natural logarithm of*MTH-36, MTH-38 
sine of • MTH-53, MTH-54 
Conjugate of complex number*MTH-44, 
MTH-45 
Cosine 

in radians*MTH-126 
Cosine 

hyperbolic*MTH-51, MTH-89 
in degrees • MTH-49, MTH—88, MTH-129 
in radians*MTH-47, MTH-87 
of complex number *MTH-26, MTH-28 


Double-precision value 
converting • MTH-62 
converting an array of*MTH-64 


E 


Exponential *MTH-66, MTH-91 

of complex number • MTH-31, MTH-33 


H 


Hyperbolic arc tangent • MTH-21, MTH-85 
Hyperbolic cosine • MTH-51, MTH-89 
Hyperbolic sine*MTH-102, MTH-134 
Hyperbolic tangent*MTH-110, MTH-142 


L 


Logarithm 

base 2*MTH-95, MTH-116 
common*MTH-97, MTH-118 
natural *MTH-93, MTH-114 
natural complex *MTH-36, MTH-38 


M 


Mathematics routine* 1-1 
absolute value* 1-4 
algorithm* 1-2 
bitwise AND operator* 1-4 
bitwise complement operator* 1-8 
bitwise exclusive OR operator* 1-5 
bitwise inclusive OR operator* 1-6 
bitwise shift* 1-8 
calling convention* 1-2 
complex number* 1-3 
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Mathematics routine (cont’d.) 
condition handling* 1-3 
conversion of double to single floating-point 
value* 1-9 

conversion to greatest floating-point integer • 
1-5 

entry point name* 1-1 

F-floating conversion • 1 -4 

floating-point conversion to nearest value* 1-7 

floating-point multiplication* 1-8 

floating-point positive difference* 1-5 

floating-point sign function* 1-8 

integer to floating-point conversion* 1-5 

JSB entry point* 1-2 

maximum value* 1-6 

minimum value* 1-7 

remainder* 1-7 

truncation of floating-point value • 1 -6 
undocumented routines* A-1 to A-15 
list of* 1-4 to 1-9 
MTH$ ACOS • MTH-3 
MTH$ ACOSD • MTH-6 
MTH$ AIM AG • MTH-112 
MTH$ ALOG • MTH-114 
MTH$ ALOG10 • MTH-118 
MTH$ AL0G2 • MTH-116 
MTH$ ASIN • MTH-9 
MTH$ ASIND • MTH-11 
MTH$ AT AN • MTH-13 
MTH$ AT AN2 • MTH-17 
MTH$ AT AND • MTH-15 
MTH$ AT AND2 • MTH-19 
MTH$ AT ANH • MTH-21 
MTH$CABS • MTH-23 
MTH$CCOS • MTH-26 
MTH$CDABS* MTH-23 
MTH$CDCOS • MTH-28 
MTH$CDEXP • MTH-33 
MTH$CDLOG • MTH-38 
MTH$CDSIN • MTH-54 
MTH$CDSQRT • MTH-59 
MTHSCEXP* MTH-31 
MTH$CGABS* MTH-23 
MTH$CGCOS* MTH-28 
MTH$CGEXP* MTH-33 
MTH$CGLOG* MTH-38 
MTH$CGSIN* MTH-54 
MTHSCGSQRT • MTH-59 
MTH$CLOG • MTH-36 
MTHSCMPLX • MTH-40 
MTH$CONJG • MTH-44 
MTH$COS • MTH-47 


MTH$COSD • MTH-49 
MTHSCOSH • MTH-51 
MTH$CSIN • MTH-53 
MTH$CSQRT • MTH-57 
MTH$CVT_DA_GA • MTH-64 
MTH$CVT_D_G • MTH-62 
MTH$CVT_GA_DA • MTH-64 
MTH$CVT_G_D • MTH-62 
MTHSDACOS* MTH-3 
MTH$D ACOSD* MTH-6 
MTHSDASIN* MTH-9 
MTH$D ASIND* MTH-11 
MTH$D AT AN • MTH-13 
MTH$D AT AN2 • MTH-17 
MTHSDAT AND • MTH-15 
MTH$D AT AND2 • MTH-19 
MTH$D ATANH • MTH-21 
MTHSDCMPLX • MTH-42 
MTHSDCONJG • MTH-45 
MTHSDCOS* MTH-47 
MTHSDCOSD* MTH-49 
MTHSDCOSH* MTH-51 
MTHSDEXP • MTH-66 
MTH$DIM AG • MTH-112 
MTH$DLOG • MTH-114 
MTH$DLOG 10 • MTH-118 
MTHSDL0G2 • MTH-116 
MTHSDRE AL • MTH-122 
MTH$DSIN • MTH-124 
MTHSDSINCOS • MTH-126 
MTHSDSINCOSD • MTH-129 
MTHSDSIND • MTH-132 
MTH$DSINH • MTH-134 
MTHSDSQRT • MTH-136 
MTH$DT AN • MTH-138 
MTHSDTAND • MTH-140 
MTH$DT ANH • MTH-142 
MTHSEXP* MTH-66 
MTHSGACOS* MTH-3 
MTH$G ACOSD* MTH-6 
MTHSGASIN* MTH-9 
MTH$G ASIND* MTH-11 
MTH$G AT AN • MTH-13 
MTH$G AT AN2 • MTH-17 
MTHSGATAND • MTH-15 
MTH$G ATAND2 • MTH-19 
MTH$G AT ANH • MTH-21 
MTHSGCMPLX • MTH-42 
MTHSGCONJG • MTH-45 
MTHSGCOS* MTH-47 
MTHSGCOSD • MTH-49 
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MTH$GCOSH • MTH-51 
MTHSGEXP • MTH-66 
MTH$GIM AG • MTH-112 
MTH$GLOG • MTH-114 
MTH$GLOG 10 • MTH-118 
MTHSGL0G2 • MTH-116 
MTH$GRE AL • MTH-122 
MTH$GSIN • MTH-124 
MTHSGSINCOS • MTH-126 
MTHSGSINCOSD • MTH-129 
MTH$GSIND • MTH-132 
MTH$GSINH • MTH-134 
MTHSGSQRT • MTH-136 
MTH$GT AN • MTH-138 
MTH$GT AND • MTH-140 
MTH$GT ANH • MTH-142 
MTH$H ACOS • MTH-69 
MTH$H ACOSD • MTH-71 
MTH$HASIN • MTH-73 
MTH$HASIND • MTH-75 
MTH$H AT AN • MTH-77 
MTH$H ATAN2 • MTH-81 
MTH$H AT AND • MTH-79 
MTH$H ATAND2 • MTH-83 
MTH$H AT ANH • MTH-85 
MTHSHCOS • MTH-87 
MTH$HCOSD • MTH-88 
MTH$HCOSH • MTH-89 
MTH$HEXP • MTH-91 
MTH$HLOG • MTH-93 
MTH$HLOG 10 • MTH-97 
MTH$HLOG2 • MTH-95 
MTHSHSIN • MTH-99 
MTH$HSINCOS • MTH-126 
MTH$HSINCOSD • MTH-129 
MTHSHSIND • MTH-100 
MTH$HSINH • MTH-102 
MTH$HSQRT • MTH-104 
MTH$HT AN • MTH-106 
MTHSHT AND • MTH-108 
MTH$HT ANH • MTH-110 
MTHSRANDOM • MTH-120 
MTH$RE AL • MTH-122 
MTH$SIN • MTH-124 
MTHSSINCOS • MTH-126 
MTH$SINCOSD • MTH-129 
MTH$SIND • MTH-132 
MTHSSINH • MTH-134 
MTH$SQRT • MTH-136 
MTH$T AN • MTH-138 
MTH$T AND • MTH-140 


MTH$T ANH • MTH-142 
MTH$UM AX • MTH-144 
MTH$UMIN • MTH-145 



Random number generator •MTH-120 



Sine 

hyperbolic • MTH-102, MTH-134 
in degrees • MTH-100, MTH-129, MTH-132 
in radians • MTH-99, MTH-124, MTH-126 
of complex number *MTH-53, MTH-54 
Square root*MTH-104, MTH-136 


T 


Tangent*MTH-106, MTH-108, MTH-138, 
MTH-140 

hyperbolic*MTH-110, MTH-142 
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